geo::Geometry Class Reference

The geometry of one entire detector (near, far, ipnd) More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-14/Geometry/Geometry.h"

Inheritance diagram for geo::Geometry:

Public Types
using	Parameters = art::ServiceTable< GeometryBase::Params >
template<class T >
using	Atom = fhicl::Atom< T >
template<class T >
using	Table = fhicl::Table< T >
using	Comment = fhicl::Comment
using	Name = fhicl::Name

Public Member Functions
	Geometry (const Parameters &params, art::ActivityRegistry &reg)
	~Geometry ()
void	preBeginRun (art::Run const &run)
TVector3	getDirectionInBeamCoordinates (TVector3 const &startPos, TVector3 const &dir) const
	Get direction of the Prong in the beam coordinates. More...
const PlaneGeo *	Plane (unsigned int i) const
const std::set< unsigned int > &	GetPlanesByView (View_t v=kXorY) const
const unsigned int	NextPlaneInView (unsigned int p, int d=+1) const
const unsigned int	NextPlaneOtherView (unsigned int p, int d=+1) const
const unsigned int	FirstPlaneInMuonCatcher () const
	Returns the index of the first plane contained in the muon catcher. More...
double	CellTpos (unsigned int ip, unsigned int ic, double w=0.0) const
void	CellInfo (unsigned int ip, unsigned int ic, View_t *view=0, double *pos=0, double *dpos=0) const
const CellGeo *	IdToCell (const CellUniqueId &id, int *iplane, int *icell) const
const PlaneGeo *	IdToPlane (const CellUniqueId &id, int *iplane) const
int	getPlaneID (const CellUniqueId &id) const
int	getPlaneID (const double &x, const double &y, const double &z) const
int	getCellID (const CellUniqueId &id) const
int	getCellID (const double &x, const double &y, const double &z) const
bool	getPlaneAndCellID (const CellUniqueId &id, int &plane, int &cell) const
bool	getPlaneAndCellID (const double &x, const double &y, const double &z, int &plane, int &cell) const
void	CountCellsOnLine (double X1, double Y1, double Z1, double X2, double Y2, double Z2, std::vector< OfflineChan > &Xhitsonline, std::vector< OfflineChan > &Yhitsonline)
void	CountCellsOnLine (TVector3 r1, TVector3 r2, std::vector< OfflineChan > &Xhitsonline, std::vector< OfflineChan > &Yhitsonline)
void	CountCellsOnLine (double X1, double Y1, double Z1, double X2, double Y2, double Z2, std::vector< CellOnLine > &Xhitsonline, std::vector< CellOnLine > &Yhitsonline)
	Make list of cells in each view traversed by a line segment and the entry and exit points into each cell. More...
void	CountCellsOnLine (TVector3 r1, TVector3 r2, std::vector< CellOnLine > &Xhitsonline, std::vector< CellOnLine > &Yhitsonline)
void	CountCellsOnLineFast (double x1, double y1, double z1, double x2, double y2, double z2, std::vector< OfflineChan > &Xhitsonline, std::vector< OfflineChan > &Yhitsonline)
	Make list of cells in each view traversed by a line segment. More...
void	CountCellsOnLineFast (TVector3 r1, TVector3 r2, std::vector< OfflineChan > &Xhitsonline, std::vector< OfflineChan > &Yhitsonline)
void	CountCellsOnLineFast (geo::View_t view, double v1, double z1, double v2, double z2, std::vector< OfflineChan > &hitsonline)
	2D version More...
unsigned int	NPlanes () const
double	DetHalfWidth () const
double	DetHalfHeight () const
double	DetLength () const
double	SurfaceY () const
TVector3	NuMIBeamDirection () const
	Direction of neutrinos from the NuMI beam (unit vector) More...
TVector3	BoosterBeamDirection () const
	Direction of neutrinos from the Booster beam (unit vector) More...
double	DistToElectronics (double localz, const CellGeo &cell) const
	get distance from local z position in cell to apd in cm, including pigtail More...
double	GetPigtail (const CellUniqueId &id) const
	Return length of fiber in cm from end of cell to apd. More...
void	WorldBox (double *xlo_cm, double *xhi_cm, double *ylo_cm, double *yhi_cm, double *zlo_cm, double *zhi_cm) const
void	DetectorEnclosureBox (double *xlo_cm, double *xhi_cm, double *ylo_cm, double *yhi_cm, double *zlo_cm, double *zhi_cm) const
void	DetectorBigBox (double *xlo_cm, double *xhi_cm, double *ylo_cm, double *yhi_cm, double *zlo_cm, double *zhi_cm) const
bool	IntersectsDetector (double *xyz_cm, double *dxyz) const
bool	IntersectsBigBox (double *xyz_cm, double *dxyz) const
bool	isDetectorBigBoxUsed () const
	Do we use the Detector Big Box. More...
bool	isInsideDetectorBigBox (const double x_cm, const double y_cm, const double z_cm) const
	Is the particle inside the detector Big Box? More...
bool	isInsideFiducialVolume (const double x_cm, const double y_cm, const double z_cm) const
	Is the particle inside the detector Fiducial Volume? More...
bool	isInsideFiducialVolume (const double *xyz_dimensionless) const
void	FiducialBox (TVector3 &r0, TVector3 &r1) const
int	CurrentCell (int *ip, int *ic) const
const CellUniqueId	CurrentCellId () const
const CellUniqueId	CellId (const double &x, const double &y, const double &z, double dxds=0., double dyds=0., double dzds=1., double step=0.01) const
TGeoMaterial *	Material (double x, double y, double z) const
bool	IsActive (const TGeoMaterial *m) const
double	TotalMass (const char *volume="vDetEnclosure") const
double	MassBetweenPoints (double *p1, double *p2) const
void	MaterialsBetweenPoints (const double *p1, const double *p2, std::vector< double > &ds, std::vector< const TGeoMaterial * > &mat) const
void	MaterialsBetweenPoints (TVector3 v1, TVector3 v2, std::vector< double > &ds, std::vector< const TGeoMaterial * > &mat) const
bool	calculateMassesLong (const unsigned int number_of_points, CLHEP::HepRandomEngine &engine) const
double	DEdge (const double *x0, const double *dx, double *exitp=0) const
double	ClosestApproach (unsigned int ip, unsigned int ic, const TVector3 &vtx, const TVector3 &dir, double dx=0.0, double dy=0.0, double *w=0, double *s=0, TVector3 *point_in_cell=0, TVector3 *point_on_track=0) const
	Find the distance of closest approach between a cell and a track. More...
std::string	GDMLFile () const
std::string	ROOTFile () const
std::string	GDMLFromFCL () const
std::string	FileBaseName () const
int	Run () const
TGeoManager *	ROOTGeoManager () const
novadaq::cnv::DetId	DetId () const
	Prefer ds::DetectorService::DetId() instead. More...
CoordinateTransformation *	getCoordinateTransformation () const
std::string	ExtractGDML () const
	Extract contents from fGDMLFile and return as a string. More...

Static Public Member Functions
static std::string	ExtractGDML (std::string fname, bool fullpath)
	Extract contents from fname and return as a string. More...

Protected Types
typedef std::unordered_map< CellUniqueId, OfflineChan >	UIDMap

Protected Member Functions
virtual void	setDetectorID (novadaq::cnv::DetId) override
	Method to set DetectorID. More...
geo::OfflineChan	getPlaneCellMap (const CellUniqueId &id) const
void	setDetectorBigBox (double detector_bigbox_range)
void	FindPlanes (std::vector< const TGeoNode * > &n, unsigned int depth, bool inMuonCatcher=false)
void	MakePlane (std::vector< const TGeoNode * > &n, unsigned int depth, bool inMuonCatcher)
void	SetDrawOptions ()
void	BuildMaps ()
void	SetDetectorSize ()
bool	LoadGeometryFile (std::string gdmlfile, novadaq::cnv::DetId det_id=novadaq::cnv::kUNKNOWN_DET)
	In case of unknown ID, guesses from the filename. More...
int	MakeTmpFile (std::string gdmlInfo)
void	RemoveTmpFile (std::string fileName)
	Method to remove temporary file that holds detector geometry file. More...

Static Protected Member Functions
static std::string	FindGDMLFile (std::string fname)
	Search for Geometry/gdml/fname in FW_SEARCH_PATH, return full path. More...

Protected Attributes
std::string	fLoadedGDMLInfo
Params	fParams
std::string	fGDMLFile
	gdml file holding the geometry More...
std::string	fROOTFile
	root file holding the geometry More...
std::string	fGDMLFromFCL
	keep track of original fcl parameter More...
novadaq::cnv::DetId	fDetId
	id for the detector being used More...
CoordinateTransformation	fCoordinateTransformation
	Coordinate Transformation class. More...
std::vector< PlaneGeo * >	fPlanes
	The detector planes. More...
std::vector< PlaneGeo * >	fPlanesInMuonCatcher
	Same pointers as fPlanes. More...
UIDMap	fIdMap
	Unique ID -> Plane,Cell. More...
std::set< unsigned int >	fAllPlanes
	List of all planes. More...
std::set< unsigned int >	fXplanes
	List of X measuring planes. More...
std::set< unsigned int >	fYplanes
	List of Y measuring planes. More...
double	fDetLength
	Detector length (cm) More...
double	fDetHalfHeight
	Detector 1/2 height (cm) More...
double	fDetHalfWidth
	Detector 1/2 width (cm) More...
bool	fIsDetectorBigBoxUsed
	Do we need to use the BigBox cut? More...
double	fBigBoxRange
	Range of big box. More...
double	fDetectorBigBoxXHi
	Dimensions of the DetectorBigBox in cm. More...
double	fDetectorBigBoxXLo
double	fDetectorBigBoxYHi
double	fDetectorBigBoxYLo
double	fDetectorBigBoxZHi
double	fDetectorBigBoxZLo
fhicl::ParameterSet	fConfigNDOS
fhicl::ParameterSet	fConfigND
fhicl::ParameterSet	fConfigFD
fhicl::ParameterSet	fConfigTB
double	fFiducialVolumeXHi
double	fFiducialVolumeXLo
double	fFiducialVolumeYHi
double	fFiducialVolumeYLo
double	fFiducialVolumeZHi
double	fFiducialVolumeZLo
std::string	fStoreTempGeo
	folder to store temp geo file More...
unsigned int	fRunNumber
	Run number of configuration. More...

Detailed Description

The geometry of one entire detector (near, far, ipnd)

Definition at line 23 of file Geometry.h.

Member Typedef Documentation

template<class T >

using geo::GeometryBase::Atom = fhicl::Atom<T>

inherited

Definition at line 51 of file GeometryBase.h.

using geo::GeometryBase::Comment = fhicl::Comment

inherited

Definition at line 53 of file GeometryBase.h.

using geo::GeometryBase::Name = fhicl::Name

inherited

Definition at line 54 of file GeometryBase.h.

using geo::Geometry::Parameters = art::ServiceTable<GeometryBase::Params>

Definition at line 27 of file Geometry.h.

template<class T >

using geo::GeometryBase::Table = fhicl::Table<T>

inherited

Definition at line 52 of file GeometryBase.h.

typedef std::unordered_map<CellUniqueId, OfflineChan> geo::GeometryBase::UIDMap

protectedinherited

Definition at line 255 of file GeometryBase.h.

Constructor & Destructor Documentation

geo::Geometry::Geometry	(	const Parameters &	params,
		art::ActivityRegistry &	reg
	)

Definition at line 30 of file Geometry_service.cc.

References geo::GeometryBase::fDetId, preBeginRun(), setDetectorID(), and art::ActivityRegistry::sPreBeginRun.

    : GeometryBase(params())
    , fLoadedGDMLInfo("(none)")
  {
    // This line is not useless as it might first appear. The call to
    // setDetectorID in the GeometryBase constructor does not call the override
    // in this class and so communicate with DetectorService. Now that the
    // derived class is more fully constructed call it again so it will.
    setDetectorID(fDetId);

    // tell the activity registry to check in after runs have been loaded
    // into the event processing so that we can update configurations if necessary
    reg.sPreBeginRun.watch(this, &Geometry::preBeginRun);
  }

geo::Geometry::~Geometry

(

)

Definition at line 46 of file Geometry_service.cc.

References geo::GeometryBase::fGDMLFile, and geo::GeometryBase::RemoveTmpFile().

  {
    // Clear up temporary files at shutdown, if necessary
    if(fGDMLFile.find("TEMP_GEO_GDML") != std::string::npos)
      RemoveTmpFile(fGDMLFile);
  }

Member Function Documentation

TVector3 geo::GeometryBase::BoosterBeamDirection ( ) const

inherited

Direction of neutrinos from the Booster beam (unit vector)

Definition at line 557 of file GeometryBase.cxx.

References e, geo::GeometryBase::fDetId, novadaq::cnv::kNDOS, and novadaq::cnv::kNEARDET.

  {
    // No simulation available for NDOS Booster. These values copied from the Cana
    // package, where the comment read:
    //
    // NDOS BNB direction from Zelimir Djurcic <zdjurcic@hep.anl.gov> 2/26/2011,
    // private communication.
    //
    // ND BNB direction from Ryan Murphy <rwmurphy@indiana.edu> 10/14/2016
    // Used same method as described in NuMIBeamDirection()
    // CAF File used:
    //   NearDet:
    //   /nova/ana/users/rwmurphy/bnb_mc_all_hadd_caf.caf.root

    switch(fDetId){
    case novadaq::cnv::kNDOS: return TVector3(-0.367, 0.012, 0.930).Unit();

    case novadaq::cnv::kNEARDET: return TVector3(-2.961e-01, -1.206e-01,+9.475e-01).Unit();

    default:
      throw cet::exception("BadGeoConfig")
        << "Bad geometry configuration: " << fDetId << "\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }
  }

void geo::GeometryBase::BuildMaps ( )

protectedinherited

Build maps used for quick access to the geometry

Definition at line 996 of file GeometryBase.cxx.

References plot_validation_datamc::c, geo::PlaneGeo::Cell(), geo::GeometryBase::fAllPlanes, geo::GeometryBase::fIdMap, geo::GeometryBase::fPlanes, geo::GeometryBase::fXplanes, geo::GeometryBase::fYplanes, MECModelEnuComparisons::i, geo::CellGeo::Id(), calib::j, geo::kX, geo::kY, geo::PlaneGeo::Ncells(), and geo::GeometryBase::Plane().

Referenced by geo::GeometryBase::LoadGeometryFile().

  {
    fAllPlanes.clear();
    fXplanes.clear();
    fYplanes.clear();

    // Sets which list planes by view
    for (unsigned int i=0; i<fPlanes.size(); ++i) {
      fAllPlanes.insert(i);
      if      (fPlanes[i]->View()==kX) fXplanes.insert(i);
      else if (fPlanes[i]->View()==kY) fYplanes.insert(i);
      else {
        throw cet::exception("BAD_GEO_CONFIG")
          << "Bad plane view" << fPlanes[i]->View() << ", plane " << i << "\n"
          << __FILE__ << ":" << __LINE__ << "\n";
      }
    }

    // Unique Id -> cell map
    for (unsigned int i=0; i<fPlanes.size(); ++i) {
      const PlaneGeo* p = fPlanes[i];
      for (unsigned int j=0; j<p->Ncells(); ++j) {
        // Assert that the cell ids are in fact unique
        if (fIdMap.find(this->Plane(i)->Cell(j)->Id())!=fIdMap.end()) {
          throw cet::exception("BAD_CELL_UNIQUE_ID")
            << "Duplicate Cell UIDs" << this->Plane(i)->Cell(j)->Id()
            << ", plane " << i << ", cell " << j << "\n"
            << __FILE__ << ":" << __LINE__ << "\n";
        }
        OfflineChan c(i,j);
        fIdMap[this->Plane(i)->Cell(j)->Id()] = c;
      }
    }
  }// end of BuildMaps

bool geo::GeometryBase::calculateMassesLong ( const unsigned int  number_of_points, CLHEP::HepRandomEngine &  engine  ) const

inherited

Method to calculate masses of the detector

Definition at line 1068 of file GeometryBase.cxx.

References CLHEP::cm, CLHEP::cm3, om::cout, barite::density, geo::GeometryBase::fDetHalfHeight, geo::GeometryBase::fDetHalfWidth, geo::GeometryBase::fDetLength, CLHEP::RandFlat::fire(), CLHEP::g, geo::GeometryBase::isInsideFiducialVolume(), CLHEP::kg, CLHEP::m3, geo::GeometryBase::Material(), geo::GeometryBase::ROOTGeoManager(), string, submit_syst::x, submit_syst::y, and test::z.

                                                                             {
    
    std::map <std::string, double> density;         // Map of densities of all materials in the detector. Dimensionless. Mapped by material names.
    std::map <std::string, double> volume;          // Map of volumes of all materials in the detector. Dimensionless. Mapped by material names.
    std::map <std::string, double> mass;            // Map of total masses of all materials the detector. Dimensionless. Mapped by material names.
    std::map <std::string, double> volume_fiducial; // Map of volumes of all materials in the fiducial region of the detector. Dimensionless. Mapped by material names.
    std::map <std::string, double> mass_fiducial;   // Map of total masses of all materials in the fiducial region of the detector. Dimensionless. Mapped by material names.
    
    CLHEP::RandFlat   flat(engine);
    
    TStopwatch timer;
    timer.Start();
    
    // Effective cell volume for a point
    const double cell_volume = (2.0*fDetHalfWidth*CLHEP::cm)*(2.0*fDetHalfHeight*CLHEP::cm)*(fDetLength*CLHEP::cm) / double(number_of_points);

    // Loop over points
    for(unsigned int ipoint = 0; ipoint < number_of_points; ++ipoint){

      // Set the progress bar
      ROOTGeoManager()->GetGeomPainter()->OpProgress("Calculating masses", ipoint, number_of_points-1, &timer, kFALSE);

      // X, Y, Z of the point in cm
      const double x = flat.fire(-fDetHalfWidth , fDetHalfWidth ) * CLHEP::cm;
      const double y = flat.fire(-fDetHalfHeight, fDetHalfHeight) * CLHEP::cm;
      const double z = flat.fire(0.0            , fDetLength    ) * CLHEP::cm;
      const double current_position[3] = {x,y,z};
      
      // Get the material of the point
      const TGeoMaterial* current_material = Material(x/CLHEP::cm,
                                                      y/CLHEP::cm,
                                                      z/CLHEP::cm
                                                      );
      // Get the material name
      const std::string current_material_name = current_material->GetName();

      // See if we have not defined the material
      if(density.find(current_material_name) == density.end()){
        density        [current_material_name] = current_material->GetDensity() * (CLHEP::g/CLHEP::cm3); //http://indra.in2p3.fr/KaliVedaDoc/1.8.1/src/KVedaLossMaterial.h.html
        // Initialize these to zero
        volume         [current_material_name] = 0.0;
        volume_fiducial[current_material_name] = 0.0;
        mass           [current_material_name] = 0.0;
        mass_fiducial  [current_material_name] = 0.0;
      }// end of searching current_material_name

      volume[current_material_name] += cell_volume;
      mass  [current_material_name] += cell_volume * density[current_material_name];
      
      // See if the point is inside fiducial volume
      if(isInsideFiducialVolume(current_position)){
        volume_fiducial[current_material_name] += cell_volume;
        mass_fiducial  [current_material_name] += cell_volume * density[current_material_name];
      }// end of checking if the point is inside fiducial volume
      
    }// end of looping over points
    
    std::cout<<"Calculated masses, volumes etc.:\n";
    
    // Printing the mass table for each detector material
    for(std::map<std::string, double>::iterator iter = density.begin(); iter != density.end(); ++iter)  {
      
      const std::string key_str = iter->first;
      
      std::cout<<std::setw(9)<<"\t\tMaterial="<<std::setw(20)<<key_str
               <<std::setw(10)<<" density="<<std::setw(10)<<density[key_str]/CLHEP::kg*CLHEP::m3<<std::setw(6)<<"kg/m3"
               <<std::setw(9)<<" volume="<<std::setw(10)<<volume[key_str]/CLHEP::m3<<std::setw(2)<<"m3"
               <<std::setw(7)<<" mass="<<std::setw(10)<<mass[key_str] / CLHEP::kg<<std::setw(2)<<"kg"
               <<std::setw(17)<<" volume_fiducial="<<std::setw(10)<<volume_fiducial[key_str] / CLHEP::m3<<std::setw(1)<<"m3"
               <<std::setw(15)<<" mass_fiducial="<<std::setw(10)<<mass_fiducial[key_str] / CLHEP::kg<<std::setw(3)<<"kg\n";
      
    }// end of printing the mass table for each detector material

    return true;
  }

const CellUniqueId geo::GeometryBase::CellId ( const double &  x, const double &  y, const double &  z, double  dxds = 0., double  dyds = 0., double  dzds = 1., double  step = 0.01  ) const

inherited

Return the unique cel id for a specified location and direction in the detector

Parameters

x	- x position coordinate (cm, world coordinates)
y	- y position coordinate (cm, world coordinates)
z	- z position coordinate (cm, world coordinates)
dxds	- dxds direction in which to explore a step
dyds	- dyds direction in which to explore a step
dzds	- dzds direction in which to explore a step
step	- step size to take

For the sake of efficiency, this method does no checking to ensure that the location is in fact inside a cell. The user is expected to perform that check before using the position (x,y,z).

Note further, that this uses the gGeoManager->InitTrack method and will disrupt TGeo-based tracking and should be called once any TGeo-based tracking is completed.

Definition at line 396 of file GeometryBase.cxx.

References geo::GeometryBase::CurrentCellId(), E, e, and std::max().

Referenced by fuzz::FuzzyKVertex::AddProng(), calib::MuondEdx::filter(), calib::MuonTrackHits::filter(), calib::MuondEdxAna::filter(), lem::FindCellRobust(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), trk::KalmanGeoHelper::MatchTrajectoryToPlane(), slid::Recluster::produce(), jmshower::RecoJMShower::RecoShowers(), g4n::FLSHitListAction::SteppingAction(), geo::GeometryTest::testCellId(), geo::GeometryTest::testCellIdFromPos(), and slid::NuEEnergyAlg::VertexEnergy().

                                                             {
    // The strategy for finding the cell ID is to try the exact point
    // we've been handed. If that does not return a valid cell ID, then
    // try shifting the point around slightly to account for tracking
    // and edge finding uncertainties. If after trying that we still
    // don't have a valid ID, return 0.

    gGeoManager->InitTrack(x, y, z, dxds, dyds, dzds);
    gGeoManager->FindNextBoundary(1.E-6);
    CellUniqueId ID  = this->CurrentCellId();
    if (ID != 0) return ID;

    double       eps = 1.0E-1;
    const double s   = std::max(step, 1.e-3);

    // walk the step points backwards and forwards to see if we can find the
    // cell
    while(eps < 1.){

      const double eps_s = eps*s;

      //     std::cout << "step back pos " << x-dxds*eps*s
      //               << "," << y-dyds*eps*s << "," << z-dzds*eps*s
      //               << " " << eps << " " << step << std::endl;
      gGeoManager->InitTrack(x-dxds*eps_s, y-dyds*eps_s, z-dzds*eps_s, dxds, dyds, dzds);
      gGeoManager->FindNextBoundary(1.E-6);
      ID = this->CurrentCellId();
      if (ID != 0) return ID;

      //     std::cout << "step forward pos " << x+dxds*eps*s
      //               << "," << y+dyds*eps*s << "," << z+dzds*eps*s << std::endl;
      gGeoManager->InitTrack(x+dxds*eps_s, y+dyds*eps_s, z+dzds*eps_s, dxds, dyds, dzds);
      ID = this->CurrentCellId();
      if (ID != 0) return ID;

      eps += 1.e-1;
    }

    // No valid ID found
    return 0;
  }// end of CellId

void geo::GeometryBase::CellInfo ( unsigned int  ip, unsigned int  ic, View_t *  view = 0, double *  pos = 0, double *  dpos = 0  ) const

inherited

Given a plane and cell number, return which view the cell measures, its nominal center position, and nominal size about that center position.

Parameters

ip	- plane number
ic	- cell number
view	- on return, the view this cell measures in
pos	- the nominal (xyz) center position of this cell (cm)
dpos	- half-dimensions of the cell in xyz (cm)

Exceptions

geo::Exception

Definition at line 813 of file GeometryBase.cxx.

References plot_validation_datamc::c, geo::PlaneGeo::Cell(), geo::CellGeo::GetCenter(), geo::CellGeo::HalfD(), geo::CellGeo::HalfL(), geo::CellGeo::HalfW(), geo::kX, geo::kY, geo::PlaneGeo::Ncells(), geo::GeometryBase::NPlanes(), geo::GeometryBase::Plane(), and geo::PlaneGeo::View().

Referenced by tbana::SpillStreamAna::analyze(), tbana::DetectorRateShutOff::analyze(), htk::HoughTrack::analyze(), comi::Leana::AnaNHit(), noe::build_cell_lookup_table(), noe::cart_to_cp(), evd::HitTower::Draw(), evd::RawDataDrawer::Draw3DHitBoxes(), skim::CosmicBeamComparison::FillNumuHistograms(), tut::TutProducer::FitPhoton(), evd::RecoBaseDrawer::GetClusterOutlines(), evd::GeometryDrawer::GetDCMBoxes(), evd::SimulationDrawer::GetLimits(), mono::MonopoleTrack::hit_is_on_track(), mono::operator<<(), earms::ElasticArmsHS::produce(), htk::Track3D::time_fit(), mono::Track3D::time_fit(), and dt::ViewMerger::ZipVerticalTracks().

                                                                              {
    if (ip>=this->NPlanes()) {
      throw cet::exception("BAD_PLANE_NUMBER")
        << " iplane " << ip << " [0," << this->NPlanes() << ")\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }
    const PlaneGeo* p = this->Plane(ip);
    if (p) {
      if (view) *view = p->View();

      if (ic>=p->Ncells()) {
        throw cet::exception("BAD_CELL_NUMBER")
          << " iplane " << ic << " [0," << p->Ncells() << ")\n"
          << __FILE__ << ":" << __LINE__ << "\n";
      }
      const CellGeo* c = p->Cell(ic);
      if (pos) c->GetCenter(pos);

      if (dpos) {
        if (p->View() == kX) {
          dpos[0] = c->HalfW();
          dpos[1] = c->HalfL();
          dpos[2] = c->HalfD();
        }
        else if (p->View() == kY) {
          dpos[0] = c->HalfL();
          dpos[1] = c->HalfW();
          dpos[2] = c->HalfD();
        }
      }
    }
  }// end of CellInfo

double geo::GeometryBase::CellTpos ( unsigned int  ip, unsigned int  ic, double  w = 0.0  ) const

inherited

Return the transverse position of the cell in which ever view its in

Parameters

ip	- plane number
ic	- cell number
w	- optional position along the local cell z-axis

Returns

- transverse position in units of cm

Exceptions

-	geo::Exception is plane, cell, or view are out of range

Definition at line 772 of file GeometryBase.cxx.

References plot_validation_datamc::c, geo::PlaneGeo::Cell(), geo::CellGeo::GetCenter(), geo::kX, geo::kY, geo::PlaneGeo::Ncells(), geo::GeometryBase::NPlanes(), geo::GeometryBase::Plane(), elec2geo::pos, and geo::PlaneGeo::View().

Referenced by mono::operator<<(), htk::Track3D::time_fit(), and mono::Track3D::time_fit().

                                                                                {
    if (ip>=this->NPlanes()) {
      throw cet::exception("BAD_PLANE_NUMBER")
        << " iplane " << ip << " [0," << this->NPlanes() << ")\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }
    const PlaneGeo* p = this->Plane(ip);
    if (p) {
      if (ic>=p->Ncells()) {
        throw cet::exception("BAD_CELL_NUMBER")
          << " iplane " << ic << " [0," << p->Ncells() << ")\n"
          << __FILE__ << ":" << __LINE__ << "\n";
      }
      const CellGeo* c = p->Cell(ic);
      if (c) {
        double pos[3];
        c->GetCenter(pos,w);
        if (p->View() == kX) { return pos[0]; }
        if (p->View() == kY) { return pos[1]; }
        throw cet::exception("BAD_PLANE_VIEW")
          << "view " << p->View() << "\n"
          << __FILE__ << ":" << __LINE__ << "\n";
        return 0.0;
      }
    }
    return 0.0;
  }// end of CellTpos

double geo::GeometryBase::ClosestApproach ( unsigned int  ip, unsigned int  ic, const TVector3 &  vtx, const TVector3 &  dir, double  dx = 0.0, double  dy = 0.0, double *  w = 0, double *  s = 0, TVector3 *  point_in_cell = 0, TVector3 *  point_on_track = 0  ) const

inherited

Find the distance of closest approach between a cell and a track.

Parameters

ip	- plane number
ic	- cell number
vtx	- track start point
dir	- track direction
dx	- offset within the cell (in cm, dx=0 is center of cell)
dy	- offset within the cell (in cm, dy=0 is center of cell)
w	- Distance from read-out end where the closest approach occurs (cm)
s	- Distance down the track where the closest approach occurs (cm)
point_in_cell	- Point inside the cell where the closest approach occurs
point_on_track	- Point on track where the closest approach occurs

Returns

The distance^2 of closest approach

Definition at line 1462 of file GeometryBase.cxx.

References getBrightness::cell, geo::PlaneGeo::Cell(), geo::ClosestApproach(), dir, dx, dy, dz, geo::CellGeo::GetCenter(), geo::CellGeo::HalfL(), geo::kX, geo::kY, PandAna.Demos.pi0_spectra::p0, geo::GeometryBase::Plane(), NDAPDHVSetting::plane, POTSpillRate::view, and geo::PlaneGeo::View().

Referenced by fuzz::ViewMatchAlg::CalcEnergyProfile(), and nerd::ViewMatchAlg::CalcEnergyProfile().

  {
    // ==
    // Find the line that runs down the cell at the request position
    //
    const PlaneGeo* plane = this->Plane(ip);
    View_t view = plane->View();;

    const CellGeo* cell = plane->Cell(ic);

    double p0[3], p1[3];
    double dz = cell->HalfL();
    cell->GetCenter(p0, +dz);
    cell->GetCenter(p1, -dz);
    
    // Shift to the ray requested
    p0[2] += dy;
    p1[2] += dy;
    if      (view==kX) { p0[0] += dx; p1[0] += dx; }
    else if (view==kY) { p0[1] += dx; p1[1] += dx; }

    TVector3 P0(p0[0], p0[1], p0[2]);
    TVector3 P1(p1[0], p1[1], p1[2]);
    // ==
    
    // Get two points on the track, the vertex and a point at vertex+dir
    TVector3 Q1(vtx); Q1 += dir;
    
    // Use the closest approach function to compute the rest
    return geo::ClosestApproach(P0,  P1, 
                                vtx, Q1, 
                                w,   s, 
                                point_in_cell,
                                point_on_track);
  }

void geo::GeometryBase::CountCellsOnLine ( double  x1, double  y1, double  z1, double  x2, double  y2, double  z2, std::vector< OfflineChan > &  xhitsonline, std::vector< OfflineChan > &  yhitsonline  )

inherited

Method to count the number of cells passed through by a line segment between two points in the detector.

Definition at line 1556 of file GeometryBase.cxx.

References plot_validation_datamc::c.

Referenced by calib::HitEfficiency::analyze(), bpfit::dEdxCalculator::computeDEDX(), and geo::GeometryBase::CountCellsOnLine().

  {
    std::vector<CellOnLine> xhs, yhs;
    CountCellsOnLine(x1, y1, z1, x2, y2, z2, xhs, yhs);

    for(const CellOnLine& c: xhs) xhitsonline.push_back(c.chan);
    for(const CellOnLine& c: yhs) yhitsonline.push_back(c.chan);
  }  // end CountCellsOnLine

void geo::GeometryBase::CountCellsOnLine ( TVector3  r1, TVector3  r2, std::vector< OfflineChan > &  xhitsonline, std::vector< OfflineChan > &  yhitsonline  )

inherited

Method to count the number of cells passed through by a line segment between two points in the detector.

Definition at line 1572 of file GeometryBase.cxx.

References geo::GeometryBase::CountCellsOnLine().

  {
    CountCellsOnLine(r1.X(), r1.Y(), r1.Z(), r2.X(), r2.Y(), r2.Z(),
                     xhitsonline, yhitsonline);
  }

void geo::GeometryBase::CountCellsOnLine ( double  x1, double  y1, double  z1, double  x2, double  y2, double  z2, std::vector< CellOnLine > &  xhitsonline, std::vector< CellOnLine > &  yhitsonline  )

inherited

Make list of cells in each view traversed by a line segment and the entry and exit points into each cell.

Method to count the number of cells passed through by a line segment between two points in the detector and keep track of the entry and exit points of the line on the cells

Definition at line 1586 of file GeometryBase.cxx.

References getBrightness::cell, geo::CellOnLine::chan, col, geo::GeometryBase::CurrentCellId(), d, geo::CellOnLine::entry, geo::CellOnLine::exit, geo::GeometryBase::IdToCell(), geo::kX, geo::kY, geo::OfflineChan::Plane(), geo::GeometryBase::Plane(), NDAPDHVSetting::plane, elec2geo::pos, util::pythag(), POTSpillRate::view, x1, and y1.

  {

    // This function returns a vector of cells crossed by the line and the
    // entry and exit points of the line through the cells for each view. If
    // the end point is within a cell, the entry or exit point is set to be
    // the end of the line.
    //
    // WARNING: This function makes use of gGeoManager and is therefore
    //            accessing global variables!  Take caution to properly
    //            initialize these global variables when accessing them after
    //            using this function!

    // define initial position and normalized direction vector
    double xyz[3]  = {x1, y1, z1};
    double d       = util::pythag(x2-x1, y2-y1, z2-z1);
    double dxyz[3] = {(x2-x1)/d, (y2-y1)/d, (z2-z1)/d};

    gGeoManager->InitTrack(xyz, dxyz);
    const double* pos = gGeoManager->GetCurrentPoint();

    double dtemp = util::pythag(pos[0]-x1, pos[1]-y1, pos[2]-z1);

    bool enter = false;

    // We update the properties of this cell intersection as we go, and store
    // it one every exit.
    CellOnLine col;

    // create a list of positions in scintillator
    while(dtemp <= d) {
      // if the current material is scintillator, count that as a cell passed through
      if(strncmp("Scintillator", gGeoManager->GetCurrentNode()->GetVolume()->GetMaterial()->GetName(), 12) == 0) {
        enter = true;

        col.entry = TVector3(pos[0],pos[1],pos[2]);

        int plane, cell;
        IdToCell(CurrentCellId(), &plane, &cell);

        col.chan = OfflineChan(plane, cell);
      }

      // step forward to the next material boundary
      gGeoManager->FindNextBoundary();
      gGeoManager->Step(kTRUE,kTRUE);

      dtemp = util::pythag(pos[0]-x1, pos[1]-y1, pos[2]-z1);

      if(enter){
        enter = false;

        if(dtemp <= d){
          col.exit = TVector3(pos[0], pos[1], pos[2]);
        }
        else{
          col.exit = TVector3(x2, y2, z2);
        }

        const View_t view = Plane(col.chan.Plane())->View();
        if(view == geo::kX){ xhitsonline.push_back(col); }
        if(view == geo::kY){ yhitsonline.push_back(col); }

      }
    }

  }  // end CountCellsOnLine

void geo::GeometryBase::CountCellsOnLine ( TVector3  r1, TVector3  r2, std::vector< CellOnLine > &  Xhitsonline, std::vector< CellOnLine > &  Yhitsonline  )

inherited

Method to count the number of cells passed through by a line segment between two points in the detector and keep track of the entry and exit points of the line on the cells

Definition at line 1661 of file GeometryBase.cxx.

References geo::GeometryBase::CountCellsOnLine().

  {
    CountCellsOnLine(r1.X(), r1.Y(), r1.Z(), r2.X(), r2.Y(), r2.Z(),
                     Xhitsonline, Yhitsonline);
  }

void geo::GeometryBase::CountCellsOnLineFast ( double  x1, double  y1, double  z1, double  x2, double  y2, double  z2, std::vector< OfflineChan > &  xhitsonline, std::vector< OfflineChan > &  yhitsonline  )

inherited

Make list of cells in each view traversed by a line segment.

Unlike CountCellsOnLine, which uses the full geometry, here we make some assumptions about cuboidal cells and the effects of alignment. For this loss in accuracy we gain a lot of speed.

Fast method to count the number of cells passed through by a line segment between two points in the detector.

Definition at line 1693 of file GeometryBase.cxx.

Referenced by chaninfo::CosmicEff::analyze(), hough::HoughValidate::analyze(), geo::LiveGeometry::CellsToEdgeInMC(), geo::GeometryBase::CountCellsOnLineFast(), dt::ViewMerger::CountUnexplainedOnLine(), chaninfo::FindSwappedChannels::filter(), geo::LiveGeometry::ProjectedCellsInMC(), geo::LiveGeometry::ProjectedCellsToEdge(), geo::LiveGeometry::ProjectedCellsToNextBadRegion(), and geo::LiveGeometry::ProjectedLiveCellsToEdge().

  {
    CountCellsOnLineFast(TVector3(x1, y1, z1), TVector3(x2, y2, z2),
                         xhitsonline, yhitsonline);
  }

void geo::GeometryBase::CountCellsOnLineFast ( TVector3  r1, TVector3  r2, std::vector< OfflineChan > &  xhitsonline, std::vector< OfflineChan > &  yhitsonline  )

inherited

Fast method to count the number of cells passed through by a line segment between two points in the detector.

Definition at line 1707 of file GeometryBase.cxx.

References getBrightness::cell, geo::PlaneGeo::Cell(), geo::ClampRayToDetectorVolume(), dz, geo::GeometryBase::fPlanes, geo::CellGeo::GetCenter(), geo::CellGeo::HalfD(), geo::CellGeo::HalfW(), geo::kX, geo::kY, std::max(), std::min(), geo::PlaneGeo::Ncells(), NDAPDHVSetting::plane, plot_timing_data::swap, POTSpillRate::view, and geo::PlaneGeo::View().

  {
    // Otherwise we would count when the ray was outside the detector in the
    // other view.
    ClampRayToDetectorVolume(&r1, &r2, this);

    // We assume downstream-going below
    if(r1.Z() > r2.Z()) std::swap(r1, r2);

    const double z1 = r1.Z();
    const double z2 = r2.Z();

    const unsigned int planeMax = fPlanes.size();
    for(unsigned int plane = 0; plane < planeMax; ++plane){
      const geo::PlaneGeo* geoplane = fPlanes[plane];
      const geo::CellGeo* geocell0 = geoplane->Cell(0);
      // Assuming here perfect alignment and identical cells across plane
      double xyzp[3];
      geocell0->GetCenter(xyzp);
      const double zc = xyzp[2];
      const double dz = geocell0->HalfD();

      if(z1 > zc+dz) continue;
      // Plane positions increase monotonically
      if(z2 < zc-dz) break;

      const geo::View_t view = geoplane->View();

      const double v1 = r1[view];
      const double v2 = r2[view];

      const double denom = z2-z1;
      const double dvdz = denom ? (v2-v1)/denom : 0;

      // Position entering the cell or starting
      const double zin = std::max(zc-dz, z1);
      double vlo = v1+(zin-z1)*dvdz;
      // Position exiting or ending
      const double zout = std::min(zc+dz, z2);
      double vhi = v1+(zout-z1)*dvdz;
      // Ensure they're from low to high instead
      if(vlo > vhi) std::swap(vlo, vhi);

      for(unsigned int cell = 0; cell < geoplane->Ncells(); ++cell){
        double xyz[3];
        geoplane->Cell(cell)->GetCenter(xyz);
        const double vc = xyz[view];
        const double dv = geoplane->Cell(cell)->HalfW();

        if(vlo > vc+dv) continue;
        // Cell positions increase monotonically
        if(vhi < vc-dv) break;

        if(view == geo::kX) xhitsonline.push_back(geo::OfflineChan(plane, cell));
        if(view == geo::kY) yhitsonline.push_back(geo::OfflineChan(plane, cell));
      } // end for cell
    } // end for plane
  }

void geo::GeometryBase::CountCellsOnLineFast ( geo::View_t  view, double  v1, double  z1, double  v2, double  z2, std::vector< OfflineChan > &  hitsonline  )

inherited

2D version

Fast method to count the number of cells passed through by a line segment between two points in the detector.

Definition at line 1674 of file GeometryBase.cxx.

References ana::assert(), geo::GeometryBase::CountCellsOnLineFast(), geo::kX, and geo::kY.

  {
    assert(view == geo::kX || view == geo::kY);

    std::vector<OfflineChan> junk;
    if(view == geo::kX)
      CountCellsOnLineFast(v1, 0, z1, v2, 0, z2, hitsonline, junk);
    else
      CountCellsOnLineFast(0, v1, z1, 0, v2, z2, junk, hitsonline);
  }

int geo::GeometryBase::CurrentCell ( int *  planeid, int *  cellid  ) const

inherited

Return the plane and cell number for the current tracking position held by the gGeoManager

Parameters

planeid	: pointer to plane number, >=0 on success, -1 on failure
cellid	: pointer to cell number, >=0 on success, -1 on failure

Returns

>0 if cell is found, <0 if not found.

Definition at line 333 of file GeometryBase.cxx.

References geo::GeometryBase::CurrentCellId(), and geo::GeometryBase::IdToCell().

Referenced by geo::GeometryTest::testFindCell().

                                                               {
    if (! this->IdToCell(this->CurrentCellId(), planeid, cellid)) {
      *planeid = -1;
      *cellid  = -1;
      throw cet::exception("UNKNOWN_ID")
        << "unknown id " << this->CurrentCellId() << "\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }
    return 1;
  }// end of

const CellUniqueId geo::GeometryBase::CurrentCellId ( ) const

inherited

Return the unique cell id for the current cell – zero if not in any cell

Definition at line 348 of file GeometryBase.cxx.

References allTimeWatchdog::endl, MECModelEnuComparisons::i, MF_LOG_DEBUG, geo::NodesToUniqueId(), ns, path, and string.

Referenced by geo::GeometryBase::CellId(), geo::GeometryBase::CountCellsOnLine(), and geo::GeometryBase::CurrentCell().

                                                       {
    const std::string path = gGeoManager->GetPath();
    //std::cout << "path to cell is " << path << std::endl;

    // check that the path actually contains a cell
    if(path.find("Cell") == std::string::npos){
      MF_LOG_DEBUG("Geometry") << "no Cell in the path to the current volume:" << path << std::endl;
      return 0;
    }
    if(path.find("HAirBubble") != std::string::npos){
      MF_LOG_DEBUG("Geometry") << "This an air bubble, not a cell:" << path << std::endl;
      return 0;
    }

    // We used to return PathToUniqueId(path.c_str()) at this point, but as of
    // the root6 upgrade, the paths are all like
    // vPlaneV_0/vModuleV_0/vExtruV_0/vCellV_0 rather than
    // vPlaneV_3/vModuleV_1/vExtruV_4/vCellV_1, so use the list of parent
    // pointers to determine the unique ID instead. The code here is derived
    // from what TGeoManager::GetPath() does.
    TGeoNodeCache* cache = gGeoManager->GetCurrentNavigator()->GetCache();
    std::vector<const TGeoNode*> ns;
    for(int i = 0; i <= cache->GetLevel(); ++i){
      ns.push_back(((TGeoNode**)cache->GetBranch())[i]);
    }
    return NodesToUniqueId(ns, ns.size()-1);
  }// end of CurrentCellId

double geo::GeometryBase::DEdge ( const double *  x0, const double *  dx, double *  exitp = 0  ) const

inherited

Compute the distance from a point to the detector wall along a particular direction.

Parameters

x0	- Input, [x,y,z] start point (cm)
dx	- Input, direction vector
exitp	- Ouput, Calculated [x,y,z] of exit point (cm)

Returns

distance to exit point in cm

Definition at line 1427 of file GeometryBase.cxx.

References geo::GeometryBase::fDetHalfHeight, geo::GeometryBase::fDetHalfWidth, geo::GeometryBase::fDetLength, geo::ProjectToBoxEdge(), and util::pythag().

Referenced by bpfit::Path::FindLayers(), comi::Cana::produce(), and crvtx::CosmicRayVertex::produce().

  {
    double xyzout[3];
    geo::ProjectToBoxEdge(x0, dx,
                          -fDetHalfWidth,  fDetHalfWidth,
                          -fDetHalfHeight, fDetHalfHeight,
                          0,               fDetLength,
                          xyzout);
    if (exitp!=0) {
      exitp[0] = xyzout[0];
      exitp[1] = xyzout[1];
      exitp[2] = xyzout[2];
    }
    return util::pythag(x0[0]-xyzout[0],
                        x0[1]-xyzout[1],
                        x0[2]-xyzout[2]);
  }// end of GeometryBase::DEdge

void geo::GeometryBase::DetectorBigBox ( double *  xlo, double *  xhi, double *  ylo, double *  yhi, double *  zlo, double *  zhi  ) const

inherited

---

Return the ranges of x,y and z for the "Detector volume" that the entire geometry lives in. If any pointers are 0, then those coordinates are ignored.

Parameters

xlo	: On return, lower bound on x positions
xhi	: On return, upper bound on x positions
ylo	: On return, lower bound on y positions
yhi	: On return, upper bound on y positions
zlo	: On return, lower bound on z positions
zhi	: On return, upper bound on z positions

Definition at line 658 of file GeometryBase.cxx.

References geo::GeometryBase::fDetectorBigBoxXHi, geo::GeometryBase::fDetectorBigBoxXLo, geo::GeometryBase::fDetectorBigBoxYHi, geo::GeometryBase::fDetectorBigBoxYLo, geo::GeometryBase::fDetectorBigBoxZHi, and geo::GeometryBase::fDetectorBigBoxZLo.

Referenced by evgen::CosmicsGen::DetectorBigBoxCut(), evgen::CORSIKAGen::DetectorBigBoxCut(), evgen::CORSIKAGen::populateNShowers(), and evgen::CosmicsGen::ProjectMuonsToDetectorBigBox().

                                                                    {
    //we don't want to return without passing along the big box coordinates
    //note the big box is set by default to just be the detector size
    //if 'isDetectorBigBox' is false
    //if(!this->isDetectorBigBoxUsed())return;
    *xlo = fDetectorBigBoxXLo;
    *xhi = fDetectorBigBoxXHi;
    *ylo = fDetectorBigBoxYLo;
    *yhi = fDetectorBigBoxYHi;
    *zlo = fDetectorBigBoxZLo;
    *zhi = fDetectorBigBoxZHi;
  }// end of DetectorBigBox

void geo::GeometryBase::DetectorEnclosureBox ( double *  xlo, double *  xhi, double *  ylo, double *  yhi, double *  zlo, double *  zhi  ) const

inherited

Return the ranges of x,y and z for the "Detector volume" that the entire geometry lives in. If any pointers are 0, then those coordinates are ignored.

Parameters

xlo	: On return, lower bound on x positions
xhi	: On return, upper bound on x positions
ylo	: On return, lower bound on y positions
yhi	: On return, upper bound on y positions
zlo	: On return, lower bound on z positions
zhi	: On return, upper bound on z positions

Definition at line 627 of file GeometryBase.cxx.

References ana::assert(), x1, submit_syst::x2, y1, and submit_syst::y2.

Referenced by g4n::FLSHitListAction::ParticleProjection(), and geo::GeometryBase::setDetectorBigBox().

                                                                           {
    const TGeoShape* s = gGeoManager->GetVolume("vDetEnclosure")->GetShape();
    assert(s);

    if (xlo || xhi) {
      double x1, x2;
      s->GetAxisRange(1,x1,x2); if (xlo) *xlo = x1; if (xhi) *xhi = x2;
    }
    if (ylo || yhi) {
      double y1, y2;
      s->GetAxisRange(2,y1,y2); if (ylo) *ylo = y1; if (yhi) *yhi = y2;
    }
    if (zlo || zhi) {
      double z1, z2;
      s->GetAxisRange(3,z1,z2); if (zlo) *zlo = z1; if (zhi) *zhi = z2;
    }
  }// end ofDetectorEnclosureBox

double geo::GeometryBase::DetHalfHeight ( ) const

inherited

Definition at line 474 of file GeometryBase.cxx.

References ana::assert(), geo::GeometryBase::fDetHalfHeight, geo::GeometryBase::fDetId, and novadaq::cnv::kUNKNOWN_DET.

Referenced by caf::AddCosmicTruthToVec(), caf::CAFMaker::AddMCTruthToVec(), fuzz::FuzzyKVertex::AddProng(), geo::GeometryTest::analyze(), calib::HitEfficiency::analyze(), mcchk::LeptonAna::analyze(), earms::ElasticArmsValidate::analyze(), mcchk::DetSimAna::analyze(), air::AirKalmanAna::analyze(), slid::LIDTraining::analyze(), trk::CosmicTrackAna::analyze(), mcchk::NeutrinoAna::AnalyzeNeutrinoInteraction(), calib::CosmicCalib::beginRun(), mcchk::CosmicAna::beginRun(), mcchk::NeutrinoAna::beginRun(), calib::MakeAttenuationProfiles::beginRun(), numue::NumuEAna::beginRun(), mcchk::LeptonAna::beginRun(), calib::AttenFit::beginRun(), calib::StopperThreshold::beginRun(), calib::AttenuationFit::beginRun(), trk::KalmanTrackMerge::CanJoinTracks(), calib::PCHitsList::CellHitToPCHit(), geo::ClampRayToDetectorVolume(), bpfit::BreakPoint::ClampToDetector(), trk::KalmanTrackAna::Containment(), trk::KalmanGeoHelper::CountMissedCellsOnLine(), evd::Display3DPad::Draw(), evd::HitTower::Draw(), evd::RawDataDrawer::Draw3DHitBoxes(), evd::RecoBaseDrawer::DrawHoughResults2D(), calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::StopperThreshold::FillHist(), calib::CosmicTrends::FillHistograms(), calib::StopperThreshold::FillTree(), filter::Filter::FillTruthVariables(), trk::KalmanTrack::FilterTracker(), vdt::VertexDT::FindMinMax(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), evgen::GENIENeutronOscGen::GENIENeutronOscGen(), evd::SliceNavigator::GetBounds(), evd::SliceNavigator::GetBox(), nuesand::FillNueSandbox::GetECF(), jmshower::RecoJMShower::GetInterCellTransDedxProb(), jmshower::RecoJMShower::GetInterPlaneDedxProb(), hough::HoughValidate::HoughLineToDetector(), bsf::BremShowerFilter::inFiducial(), dif::DiFShowerFinder::inFiducial(), g4n::NeutronSubstitutionProcess::InsideDetector(), calib::RockMuonStopperSelection::IsContainedStopper(), comi::CosmicMetrics::IsEntering(), align::AlignInspect::IsGoodTrack(), align::Alignment::IsGoodTrack(), comi::Leana::IsPointOnEdge(), calib::RockMuonStopperSelection::IsTrulyContainedStopper(), calib::RockMuonStopperSelection::IsTrulyInFiducial(), calib::ThroughgoingSelection::IsTrulyUncontained(), calib::ThroughgoingSelection::IsUncontained(), evd::TZProjPad::LayoutHisto(), calib::AttenFit::lowessFit(), calib::AttenuationFit::lowessFit(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), trk::KalmanGeoHelper::MatchTrajectoryToPlane(), numue::NumuEAlg::NumuEAlg(), cheat::TrackCheater::produce(), jmshower::JMTrackMerge::produce(), calib::FiberCalibration::produce(), earms::ElasticArmsHS::produce(), comi::Cana::produce(), evgen::GENIERockGen::produce(), crvtx::CosmicRayVertex::produce(), cheat::MCCheater::produce(), remid::RecoMuon::produce(), evd::RawDataDrawer::RawDigit3D(), trk::RLFit::RLFit(), calib::AttenFit::rolloffFit(), calib::AttenuationFit::rolloffFit(), evd::TZProjPad::ShowFull(), evd::TZProjPad::ShowPartial(), numue::NumuEAna::TotalLengthInDetector(), cheat::RecoCheckAna::TotalLengthInDetector(), and trk::KalmanTrackAna::TotalLengthInDetector().

  {
    assert(fDetId != novadaq::cnv::kUNKNOWN_DET);
    return fDetHalfHeight;
  }

double geo::GeometryBase::DetHalfWidth ( ) const

inherited

Definition at line 467 of file GeometryBase.cxx.

References ana::assert(), geo::GeometryBase::fDetHalfWidth, geo::GeometryBase::fDetId, and novadaq::cnv::kUNKNOWN_DET.

Referenced by caf::AddCosmicTruthToVec(), caf::CAFMaker::AddMCTruthToVec(), fuzz::FuzzyKVertex::AddProng(), geo::GeometryTest::analyze(), calib::HitEfficiency::analyze(), mcchk::LeptonAna::analyze(), earms::ElasticArmsValidate::analyze(), mcchk::DetSimAna::analyze(), air::AirKalmanAna::analyze(), slid::LIDTraining::analyze(), trk::CosmicTrackAna::analyze(), mcchk::NeutrinoAna::AnalyzeNeutrinoInteraction(), calib::CosmicCalib::beginRun(), mcchk::CosmicAna::beginRun(), mcchk::NeutrinoAna::beginRun(), numue::NumuEAna::beginRun(), calib::MakeAttenuationProfiles::beginRun(), mcchk::LeptonAna::beginRun(), calib::AttenFit::beginRun(), calib::StopperThreshold::beginRun(), calib::AttenuationFit::beginRun(), trk::KalmanTrackMerge::CanJoinTracks(), calib::PCHitsList::CellHitToPCHit(), geo::ClampRayToDetectorVolume(), bpfit::BreakPoint::ClampToDetector(), trk::KalmanTrackAna::Containment(), trk::KalmanGeoHelper::CountMissedCellsOnLine(), evd::Display3DPad::Draw(), evd::HitTower::Draw(), evd::RawDataDrawer::Draw3DHitBoxes(), evd::RecoBaseDrawer::DrawHoughResults2D(), calib::SystematicUncertaintyScaleAna::endJob(), calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), calib::StopperThreshold::FillHist(), calib::CosmicTrends::FillHistograms(), calib::StopperThreshold::FillTree(), filter::Filter::FillTruthVariables(), trk::KalmanTrack::FilterTracker(), vdt::VertexDT::FindMinMax(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), evgen::GENIENeutronOscGen::GENIENeutronOscGen(), evd::SliceNavigator::GetBounds(), evd::SliceNavigator::GetBox(), nuesand::FillNueSandbox::GetECF(), airshower::AirSlicer::GetHoughAngles(), airshower::AirSlicer::GetHoughMap(), airshower::AirSlicer::GetHoughRhos(), jmshower::RecoJMShower::GetInterCellTransDedxProb(), jmshower::RecoJMShower::GetInterPlaneDedxProb(), hough::HoughValidate::HoughLineToDetector(), bsf::BremShowerFilter::inFiducial(), dif::DiFShowerFinder::inFiducial(), g4n::NeutronSubstitutionProcess::InsideDetector(), calib::RockMuonStopperSelection::IsContainedStopper(), comi::CosmicMetrics::IsEntering(), align::AlignInspect::IsGoodTrack(), align::Alignment::IsGoodTrack(), comi::Leana::IsPointOnEdge(), calib::RockMuonStopperSelection::IsTrulyContainedStopper(), calib::RockMuonStopperSelection::IsTrulyInFiducial(), calib::ThroughgoingSelection::IsTrulyUncontained(), calib::ThroughgoingSelection::IsUncontained(), evd::TZProjPad::LayoutHisto(), calib::AttenFit::lowessFit(), calib::AttenuationFit::lowessFit(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), trk::KalmanGeoHelper::MatchTrajectoryToPlane(), numue::NumuEAlg::NumuEAlg(), cheat::TrackCheater::produce(), jmshower::JMTrackMerge::produce(), calib::FiberCalibration::produce(), earms::ElasticArmsHS::produce(), comi::Cana::produce(), evgen::GENIERockGen::produce(), cheat::MCCheater::produce(), crvtx::CosmicRayVertex::produce(), remid::RecoMuon::produce(), evd::RawDataDrawer::RawDigit3D(), trk::RLFit::RLFit(), calib::AttenFit::rolloffFit(), calib::AttenuationFit::rolloffFit(), evd::TZProjPad::ShowFull(), evd::TZProjPad::ShowPartial(), numue::NumuEAna::TotalLengthInDetector(), cheat::RecoCheckAna::TotalLengthInDetector(), and trk::KalmanTrackAna::TotalLengthInDetector().

  {
    assert(fDetId != novadaq::cnv::kUNKNOWN_DET);
    return fDetHalfWidth;
  }

novadaq::cnv::DetId geo::GeometryBase::DetId ( ) const

inlineinherited

Prefer ds::DetectorService::DetId() instead.

Definition at line 243 of file GeometryBase.h.

Referenced by calib::ADCShapeFitAna::analyze(), geo::GeometryTest::analyze(), align::SplitTracks::analyze(), comi::DataCheck::analyze(), bpfit::BPFTmvaTrainer::analyze(), LightLevels::analyze(), CerenkovResponse::analyze(), util::StoreGDML::beginRun(), calib::ThresholdAna::beginRun(), slicer::Slicer4D::beginRun(), presel::RockPresel::beginRun(), presel::CosmicVeto::beginRun(), wsnumu::WSNumu::beginRun(), presel::CosmicCVNVeto::beginRun(), presel::NuePresel::beginRun(), mcchk::CloseInTime::beginRun(), dprf::ChannelPlots::beginRun(), comi::NumiFiltering::beginRun(), mcchk::DetAna::beginRun(), calib::CosmicCalib::beginRun(), dprf::TrackPlots::beginRun(), calib::FiberCalibration::beginRun(), evgen::CosmicsGen::beginRun(), evgen::CosmicPionGen::beginRun(), evgen::IsoGen::beginRun(), numue::NumuEAna::beginRun(), calib::MakeAttenuationProfiles::beginRun(), supernova::SnovaGen::beginRun(), evgen::SupernovaGen::beginRun(), comi::NumiFilteringAna::beginRun(), evgen::GENIERockGen::beginRun(), calib::TimingCalibration::beginRun(), evgen::SingleGen::beginRun(), remid::ReMIdDedx::beginRun(), chaninfo::BadChanValidate::beginRun(), dif::DiFShowerFinder::beginRun(), calib::DriftResponseCalc::beginRun(), calib::PEResponse::beginRun(), evgen::CORSIKAGen::beginRun(), evgen::GENIEGen::beginRun(), calib::AttenFit::beginRun(), jmshower::NueSel::beginRun(), jmshower::NueSelLID::beginRun(), slid::SPIDBuilder::beginRun(), comi::NearlineAna::beginRun(), remid::ReMIdDedxRock::beginRun(), murem::MuonRemove::beginRun(), photrans::ImprovedTransport::beginRun(), calib::AttenuationFit::beginRun(), remid::ReMIdDedxFD::beginRun(), remid::ReMIdDedxStudies::beginRun(), rsim::ReadoutSim::beginRun(), calib::DCMTimingOffset::beginRun(), evgen::GENIENeutronOscGen::beginRun(), slid::LIDTraining::beginRun(), jmshower::RecoJMShower::beginRun(), slicemergeviews::SliceMergeViews::beginRun(), evgen::TextFileGen::beginRun(), tdslicer::TDSlicer::beginRun(), caf::CAFMaker::beginRun(), DataCheck::HardwareStatusCheck::beginSubRun(), comi::DataCheck::beginSubRun(), slicer::Slicer::BlockSlice(), slid::NuEEnergyAlg::CellADCToGeV(), slid::NuEEnergyAlg::CellEDataMC(), slid::NuEEnergyAlg::CellEnergy(), calib::PCHitsList::CellHitToPCHit(), geo::LiveGeometry::CellsToEdgeInMC(), calib::AttenFit::ConsolidateViews(), calib::AttenuationFit::ConsolidateViews(), remid::ReMIdTrain::ContainedEvent(), chaninfo::BadChanList::CountDCMHits(), rsim::ReadoutSim::CreateRawDigits(), evgen::CosmicsGen::DetectorBigBoxCut(), evgen::CORSIKAGen::DetectorBigBoxCut(), geo::LiveGeometry::DistanceToEdgeInMC(), geo::LiveGeometry::DistToClosestBadRegion(), geo::LiveGeometry::DoWeEnterMC(), evd::HardwarePad::Draw(), evd::GeometryDrawer::DrawCells2D(), calib::SystematicUncertaintyScaleAna::endJob(), comi::CosmicMetrics::endJob(), comi::DataCheck::endJob(), calib::AttenFit::endJob(), calib::DCMTimingOffset::endJob(), calib::Calibrator::EnsureTimingCacheVldTime(), calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), geo::LiveGeometry::FillBadBoxes(), dqsf::DQSpillFlags::FillCountVars(), meta::Metadata::FillMetadata(), calib::CosmicCalib::FillPathLengthInfo(), calib::MakeAttenuationProfiles::FillPathLengthInfo(), util::RemoveBeamSpills::filter(), util::RemoveBeamSpillsPrecise::filter(), comi::NumiFiltering::filter(), calib::MuondEdx::filter(), calib::MuonTrackHits::filter(), calib::MuondEdxAna::filter(), calib::AttenFit::fit_channel_prof(), calib::AttenuationFit::fit_channel_prof(), geo::LiveGeometry::FullNDProjectedCells(), calib::Calibrator::GetAdcPerPE(), calib::Calibrator::GetAttenCurve(), calib::Calibrator::GetAttenScale(), geo::LiveGeometry::GetBCInfo(), LSTME::SliceVars::getCoarseTiming(), evd::GeometryDrawer::GetDCMBoxes(), geo::LiveGeometry::GetDetectorEdges(), skim::ParametersNumu::GetDetectorEdges(), LSTME::SliceVars::getLowGain(), numusand::NumuSandFxs::getMissingE(), calib::Calibrator::GetPECorrToGeVScale(), geo::GeometryBase::GetPigtail(), photrans::ImprovedTransport::GetPosCorr(), geo::LiveGeometry::GetRHInfo(), calib::Calibrator::GetTimeRes(), calib::Calibrator::GetTimingOffset(), calib::Calibrator::GetTNS(), calib::ThresholdAna::GetView(), calib::CosmicCalib::GetView(), calib::MakeAttenuationProfiles::GetView(), calib::AttenFit::GetView(), calib::AttenuationFit::GetView(), calib::AttenCache::GetView(), chaninfo::BadChanList::InDroppedDCM(), mcdatadaq::Configuration::init(), qeef::QeFinder::Init(), remid::RecoMuon::Init(), slid::SPIDAlg::Initialize(), slid::LIDAlg::Initialize(), geo::LiveGeometry::InMuonCatcher(), geo::LiveGeometry::InstrumentedDetBack(), geo::LiveGeometry::InstrumentedDetEnds(), geo::LiveGeometry::InstrumentedDetFront(), geo::LiveGeometry::InstrumentedDetLength(), chaninfo::BadChanList::IsBadUncached(), calib::RockMuonStopperSelection::IsContainedStopper(), calib::StopperSelection::IsContainedStopper(), remid::RecoMuon::IsTrackContained(), filter::Filter::isVertexInFiducial(), evd::TZProjPad::LayoutHisto(), calib::AttenCache::LoadFromDatabase(), calib::AttenCache::LoadFromDatabaseUncalibratedMask(), chaninfo::BadChanList::LoadNewData(), calib::AttenFit::lowessFit(), calib::AttenuationFit::lowessFit(), calib::Calibrator::MakeCellHit(), evtsum::EventSummary::MakeOutput(), chaninfo::BadChanList::NBadInSubRun(), slicer::Slicer::NeighborSlice(), numue::NumuEAlg::NumuEAlg(), trk::CosmicTrackUtilities::PointInDetector(), calib::Calibrator::postBeginRun(), preBeginRun(), SliceLID::FillSliceLID::predict(), LSTME::FillLSTME::predictEnergy(), wsnumu::WSNumu::produce(), calib::ResizeDetector::produce(), calib::FiberCalibration::produce(), cosrej::CosRej::produce(), numue::NumuEnergy::produce(), calib::TimingCalibration::produce(), rwgt::MakeGENIEReweightTable::produce(), calib::DriftResponseCalc::produce(), slid::LIDBuilder::produce(), dqsf::DQSpillFlags::produce(), bpfit::BPFEnergyEstimator::produce(), remid::RecoMuon::produce(), bpfit::BPFPIdMaker::produce(), g4rwgt::MakeG4WeightTable::produce(), rsim::ReadoutSim::produce(), evgen::CosmicsGen::ProjectCosmicsToSurface(), geo::LiveGeometry::ProjectedAirDist(), geo::LiveGeometry::ProjectedCellsInMC(), geo::LiveGeometry::ProjectedDistanceInMC(), geo::LiveGeometry::ProjectedSteelDist(), evgen::CosmicsGen::ProjectMuonsToDetectorBigBox(), evgen::CosmicsGen::ProjectToSurface(), evd::RawDataDrawer::RawDigitHardware(), lem::GenFromLib::readNext(), cheat::BackTracker::Rebuild(), calib::AttenCache::RepresentativePlane(), calib::AttenFit::rolloffFit(), calib::AttenuationFit::rolloffFit(), LSTME::FillLSTME::selectModel(), geo::LiveGeometry::SetInstrumentation(), trk::KalmanGeoHelper::SetupHelper(), geo::LiveGeometry::SetupLiveGeo(), photrans::ImprovedTransport::StepAlongHit(), calib::Calibrator::SystematicUncertaintyScale(), calib::AttenCache::TranslateChannelData(), calib::AttenCache::TranslateChannelMC(), calib::PCHitsList::TrueW(), evd::TZProjPad::TZProjPad(), slid::NuEEnergyAlg::VertexEnergy(), rb::Cluster::W(), rb::Prong::W(), rb::Track::W(), calib::AttenuationFit::WriteDummyCSVRow(), calib::AttenuationFit::writeResults(), geo::LiveGeometry::XPositionAtTransition(), and geo::LiveGeometry::YPositionAtTransition().

{ return fDetId; }

double geo::GeometryBase::DetLength ( ) const

inherited

Definition at line 481 of file GeometryBase.cxx.

References ana::assert(), geo::GeometryBase::fDetId, geo::GeometryBase::fDetLength, and novadaq::cnv::kUNKNOWN_DET.

Referenced by caf::AddCosmicTruthToVec(), caf::CAFMaker::AddMCTruthToVec(), fuzz::FuzzyKVertex::AddProng(), mcchk::CosmicAna::analyze(), mcchk::LeptonAna::analyze(), earms::ElasticArmsValidate::analyze(), mcchk::DetSimAna::analyze(), air::AirKalmanAna::analyze(), slid::LIDTraining::analyze(), trk::CosmicTrackAna::analyze(), mcchk::NeutrinoAna::AnalyzeNeutrinoInteraction(), mcchk::CosmicAna::beginRun(), mcchk::NeutrinoAna::beginRun(), numue::NumuEAna::beginRun(), mcchk::LeptonAna::beginRun(), calib::PCHitsList::CellHitToPCHit(), geo::ClampRayToDetectorVolume(), bpfit::BreakPoint::ClampToDetector(), trk::KalmanTrackAna::Containment(), trk::KalmanGeoHelper::CountMissedCellsOnLine(), evd::Display3DPad::Draw(), evd::RecoBaseDrawer::DrawHoughResults2D(), calib::CosmicTrends::FillHistograms(), filter::Filter::FillTruthVariables(), trk::KalmanTrack::FilterTracker(), vdt::VertexDT::FindMinMax(), evgen::GENIENeutronOscGen::GENIENeutronOscGen(), evd::SliceNavigator::GetBounds(), evd::SliceNavigator::GetBox(), nuesand::FillNueSandbox::GetECF(), airshower::AirSlicer::GetHoughMap(), airshower::AirSlicer::GetHoughRhos(), hough::HoughValidate::HoughLineToDetector(), g4n::NeutronSubstitutionProcess::InsideDetector(), calib::RockMuonStopperSelection::IsContainedStopper(), comi::CosmicMetrics::IsEntering(), align::AlignInspect::IsGoodTrack(), align::Alignment::IsGoodTrack(), calib::RockMuonStopperSelection::IsTrulyContainedStopper(), calib::ThroughgoingSelection::IsTrulyUncontained(), calib::ThroughgoingSelection::IsUncontained(), evd::TZProjPad::LayoutHisto(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), numue::NumuEAlg::NumuEAlg(), cheat::TrackCheater::produce(), evgen::CosmicsGen::produce(), earms::ElasticArmsHS::produce(), evgen::CosmicPionGen::produce(), comi::Cana::produce(), evgen::GENIERockGen::produce(), vdt::VertexDT::produce(), crvtx::CosmicRayVertex::produce(), cheat::MCCheater::produce(), remid::RecoMuon::produce(), geo::LiveGeometry::ProjectedWall(), trk::RLFit::RLFit(), evd::TZProjPad::ShowFull(), evd::TZProjPad::ShowPartial(), numue::NumuEAna::TotalLengthInDetector(), cheat::RecoCheckAna::TotalLengthInDetector(), trk::KalmanTrackAna::TotalLengthInDetector(), and dt::Chain::ToTrack().

  {
    assert(fDetId != novadaq::cnv::kUNKNOWN_DET);
    return fDetLength;
  }

double geo::GeometryBase::DistToElectronics ( double  localz, const CellGeo &  cell  ) const

inherited

get distance from local z position in cell to apd in cm, including pigtail

Return distance from a point in the cell to readout electronics, including pigtail

Parameters

localz

: Z position in cm in local cell frame

Definition at line 1771 of file GeometryBase.cxx.

References geo::GeometryBase::GetPigtail(), geo::CellGeo::HalfL(), and geo::CellGeo::Id().

Referenced by calib::HitEfficiency::analyze(), jmshower::RecoJMShower::GetTrkHitPos(), calib::GetXYZD(), and slid::NuEEnergyAlg::VertexEnergy().

  {
    return (cell.HalfL() - localz + this->GetPigtail(cell.Id()));
  }

std::string geo::GeometryBase::ExtractGDML ( ) const

inherited

Extract contents from fGDMLFile and return as a string.

Definition at line 180 of file GeometryBase.cxx.

References geo::GeometryBase::fGDMLFile.

Referenced by util::StoreGDML::beginRun(), evgen::CosmicsGen::beginRun(), evgen::CosmicPionGen::beginRun(), evgen::IsoGen::beginRun(), supernova::SnovaGen::beginRun(), evgen::SupernovaGen::beginRun(), evgen::SingleGen::beginRun(), evgen::CORSIKAGen::beginRun(), evgen::GENIEGen::beginRun(), evgen::GENIENeutronOscGen::beginRun(), daq2raw::NOvARawInputDriver::readNext(), and daq2raw::TestBeamRawInputDriver::readNext().

  {
    return ExtractGDML(fGDMLFile, true);
  }

std::string geo::GeometryBase::ExtractGDML ( std::string  fname, bool  fullpath  )

staticinherited

Extract contents from fname and return as a string.

Parameters

fullpath

If false, search for the file in FW_SEARCH_PATH

Definition at line 186 of file GeometryBase.cxx.

References geo::GeometryBase::FindGDMLFile(), make_syst_table_plots::line, and string.

  {
    if(!fullpath) fname = FindGDMLFile(fname);

    std::string gdmlInfo;

    std::ifstream geoFile(fname.c_str());
    std::string line;
    
    if (geoFile.is_open()){
      while (geoFile.good()){
        getline(geoFile, line);
        gdmlInfo += line;
        gdmlInfo += '\n';
      }
      geoFile.close();
    }
    
   return gdmlInfo;
  }

void geo::GeometryBase::FiducialBox ( TVector3 &  r0, TVector3 &  r1  ) const

inherited

Definition at line 1410 of file GeometryBase.cxx.

References geo::GeometryBase::fFiducialVolumeXHi, geo::GeometryBase::fFiducialVolumeXLo, geo::GeometryBase::fFiducialVolumeYHi, geo::GeometryBase::fFiducialVolumeYLo, geo::GeometryBase::fFiducialVolumeZHi, and geo::GeometryBase::fFiducialVolumeZLo.

  {
    r0 = TVector3(fFiducialVolumeXLo, fFiducialVolumeYLo, fFiducialVolumeZLo);
    r1 = TVector3(fFiducialVolumeXHi, fFiducialVolumeYHi, fFiducialVolumeZHi);
  }

std::string geo::GeometryBase::FileBaseName ( ) const

inherited

Method to extract gdml file basename and store as an std::string

Definition at line 307 of file GeometryBase.cxx.

References ana::assert(), geo::GeometryBase::fGDMLFile, and string.

Referenced by geo::GeometryTest::analyze(), util::StoreGDML::beginRun(), evgen::CosmicsGen::beginRun(), evgen::CosmicPionGen::beginRun(), evgen::IsoGen::beginRun(), supernova::SnovaGen::beginRun(), evgen::SupernovaGen::beginRun(), evgen::GENIERockGen::beginRun(), evgen::SingleGen::beginRun(), evgen::CORSIKAGen::beginRun(), evgen::GENIEGen::beginRun(), evgen::GENIENeutronOscGen::beginRun(), evgen::TextFileGen::beginRun(), and lem::GenFromLib::readNext().

                                             {
    assert(fGDMLFile.find(".gdml") != std::string::npos);

    // to find the base name of the file
    // use string operations to find the last
    // instance of "/" in the fROOTFile string
    // and the last instance of ".", the base
    // string should be between those
    size_t lastSlash   = fGDMLFile.find_last_of("/");
    std::string almost = fGDMLFile.substr(lastSlash+1);

    // resize almost to drop the file ".gdml"
    almost.resize(almost.length() - 5);

    return almost;
  }// end of FileBaseName

std::string geo::GeometryBase::FindGDMLFile ( std::string  fname )

staticprotectedinherited

Search for Geometry/gdml/fname in FW_SEARCH_PATH, return full path.

Definition at line 208 of file GeometryBase.cxx.

References cet::search_path::find_file(), runNovaSAM::ret, and string.

Referenced by geo::GeometryBase::ExtractGDML(), and preBeginRun().

  {
    cet::search_path sp("FW_SEARCH_PATH");
    size_t lastSlash   = fname.find_last_of("/");
    std::string basename = fname.substr(lastSlash+1);

    fname = "Geometry/gdml/"+basename;

    std::string ret;
    if(!sp.find_file(fname, ret)){
      throw cet::exception("FindGDMLFile")
        << "cannot find GDML file " << fname << " bail ungracefully\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }
    return ret;
  }

void geo::GeometryBase::FindPlanes ( std::vector< const TGeoNode * > &  n, unsigned int  depth, bool  inMuonCatcher = false  )

protectedinherited

Recursively search for planes in the geometry

Parameters

n	- Array holding the path of nodes traversed
depth	- Depth of the search to this point
inMuonCatcher	- internal flag, recursed into the muon catcher?

Exceptions

geo::Exception

is maximum depth is exceeded.

Definition at line 703 of file GeometryBase.cxx.

References MECModelEnuComparisons::i, geo::GeometryBase::MakePlane(), nd, Munits::nm, and registry_explorer::v.

Referenced by geo::GeometryBase::LoadGeometryFile().

  {
    const char* nm = n[depth]->GetName();
    if (strncmp(nm,"vPlane", 6)==0) {
      this->MakePlane(n, depth, inMuonCatcher);
      return;
    }

    if(strncmp(nm, "vBlockMuon", 10) == 0) inMuonCatcher = true;

    // Explore the next layer down
    unsigned int deeper = depth+1;
    if (deeper>=n.size()) {
      throw cet::exception("BAD_NODE")
        << "Exceeded maximum depth (" << n.size() << ") " << deeper << "\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }
    const TGeoVolume* v = n[depth]->GetVolume();
    int nd = v->GetNdaughters();
    for (int i=0; i<nd; ++i) {
      n[deeper] = v->GetNode(i);
      this->FindPlanes(n, deeper, inMuonCatcher);
    }
  }// end of FindPlanes

const unsigned int geo::GeometryBase::FirstPlaneInMuonCatcher ( ) const

inherited

Returns the index of the first plane contained in the muon catcher.

Or kPLANE_NOT_FOUND if there is no muon catcher

Definition at line 979 of file GeometryBase.cxx.

References ana::assert(), geo::GeometryBase::fPlanes, geo::GeometryBase::fPlanesInMuonCatcher, it, and geo::kPLANE_NOT_FOUND.

Referenced by LightLevels::analyze(), CerenkovResponse::analyze(), calib::SystematicUncertaintyScaleAna::endJob(), calib::AttenuationFit::expFit(), dqsf::DQSpillFlags::FillCountVars(), calib::MuondEdx::filter(), calib::MuonTrackHits::filter(), calib::MuondEdxAna::filter(), calib::AttenuationFit::fit_channel_prof(), calib::Calibrator::GetAttenCurve(), calib::ThresholdAna::GetView(), calib::CosmicCalib::GetView(), calib::MakeAttenuationProfiles::GetView(), calib::AttenFit::GetView(), calib::AttenuationFit::GetView(), calib::AttenCache::GetView(), remid::RecoMuon::IsTrackContained(), calib::AttenCache::LoadFromDatabase(), calib::AttenuationFit::lowessFit(), numue::NumuEAlg::NumuEAlg(), calib::AttenCache::RepresentativePlane(), calib::AttenuationFit::rolloffFit(), photrans::ImprovedTransport::StepAlongHit(), calib::Calibrator::SystematicUncertaintyScale(), calib::AttenCache::TranslateChannelMC(), calib::AttenuationFit::WriteDummyCSVRow(), and calib::AttenuationFit::writeResults().

  {
    if(fPlanesInMuonCatcher.empty()){
      mf::LogWarning("GeometryBase") << "Requesting muon catcher info for a detector without one";
      return kPLANE_NOT_FOUND;
    }

    auto it = std::find(fPlanes.begin(), fPlanes.end(), fPlanesInMuonCatcher[0]);

    assert(it != fPlanes.end());
    return it-fPlanes.begin();
  }

std::string geo::GeometryBase::GDMLFile ( ) const

inlineinherited

Definition at line 235 of file GeometryBase.h.

Referenced by geo::GeometryTest::analyze(), g4n::G4Alg::G4Alg(), and rawp::CheckDAQChannelMap::Reco().

{ return fGDMLFile;  }

std::string geo::GeometryBase::GDMLFromFCL ( ) const

inlineinherited

Definition at line 237 of file GeometryBase.h.

References string.

{ return fGDMLFromFCL; }

int geo::GeometryBase::getCellID ( const CellUniqueId &  id ) const

inherited

Definition at line 888 of file GeometryBase.cxx.

References geo::OfflineChan::Cell(), and geo::GeometryBase::getPlaneCellMap().

                                                           {
    return getPlaneCellMap(id).Cell();
  }

int geo::GeometryBase::getCellID ( const double &  x, const double &  y, const double &  z  ) const

inlineinherited

Definition at line 107 of file GeometryBase.h.

References geo::GeometryBase::getCellID().

Referenced by geo::GeometryBase::getCellID().

{return getCellID(CellId(x,y,z));}

CoordinateTransformation* geo::GeometryBase::getCoordinateTransformation ( ) const

inlineinherited

Definition at line 245 of file GeometryBase.h.

References plot_validation_datamc::fname, submit_nova_art::fullpath, and string.

Referenced by mcchk::NeutrinoAna::AnalyzeNeutrinoInteraction(), and getDirectionInBeamCoordinates().

{ return (CoordinateTransformation*)&fCoordinateTransformation;}

TVector3 geo::Geometry::getDirectionInBeamCoordinates	(	TVector3 const &	startPos,
		TVector3 const &	dir
	)		const

Get direction of the Prong in the beam coordinates.

x0 in Detector coordinates. All dimensions are in cm

Small displacement in the Prong's direction

x1 in Detector coordinates. All dimensions are in cm

x0, x1 in Detector coordinates. All dimensions are CLHEP

x0 and x1 in beam coordinates

Calculate beam coordinates for x0 and x1

Transform into a vector class

Direction in beam coordinates. Not normalized

return Unit vector

Definition at line 154 of file Geometry_service.cc.

References CLHEP::cm, geo::GeometryBase::getCoordinateTransformation(), and geo::CoordinateTransformation::transformDetectorToBeamCoordinates().

  {

    /// x0 in Detector coordinates. All dimensions are in cm
    const TVector3 x0_det_cm = startPos;

    /// Small displacement in the Prong's direction
    const TVector3 dr        = dir * 0.001;

    /// x1 in Detector coordinates. All dimensions are in cm
    const TVector3 x1_det_cm = x0_det_cm + dr;

    /// x0, x1 in Detector coordinates. All dimensions are CLHEP
    double x0_det[3] = {x0_det_cm.X()*CLHEP::cm
                       ,x0_det_cm.Y()*CLHEP::cm
                       ,x0_det_cm.Z()*CLHEP::cm
                       };
    double x1_det[3] = {x1_det_cm.X()*CLHEP::cm
                       ,x1_det_cm.Y()*CLHEP::cm
                       ,x1_det_cm.Z()*CLHEP::cm
                       };

    /// x0 and x1 in beam coordinates
    double x0_beam[3] = {0.0, 0.0, 0.0};
    double x1_beam[3] = {0.0, 0.0, 0.0};

    /// Calculate beam coordinates for x0 and x1
    const geo::CoordinateTransformation* coordinate_transformation = getCoordinateTransformation();
    coordinate_transformation->transformDetectorToBeamCoordinates(x0_det, x0_beam);
    coordinate_transformation->transformDetectorToBeamCoordinates(x1_det, x1_beam);

    /// Transform into a vector class
    const TVector3 x0_beam_vec(x0_beam);

    /// Direction in beam coordinates. Not normalized
    TVector3 dir_beam(x1_beam);
    dir_beam -= x0_beam_vec;

    /// return Unit vector
    return dir_beam.Unit();
  }

double geo::GeometryBase::GetPigtail ( const CellUniqueId &  id ) const

inherited

Return length of fiber in cm from end of cell to apd.

Definition at line 1777 of file GeometryBase.cxx.

References geo::GeometryBase::DetId(), geo::GeometryBase::IdToCell(), geo::GeometryBase::IdToPlane(), novadaq::cnv::kTESTBEAM, geo::kY, and POTSpillRate::view.

Referenced by geo::GeometryBase::DistToElectronics(), photrans::PhotonTransport::produce(), jmshower::RecoJMShower::RecoShowers(), and photrans::ImprovedTransport::StepAlongHit().

  {
    int planeid,cellid;
    this->IdToCell(id, &planeid, &cellid);
    View_t view = (this->IdToPlane(id,&planeid))->View();
    //NOTE: Do we think 32 cells per module can ever change
    //This number can be found in DAQChannelMap but did not want to add that
    //dependency to geometry unless needed.
    const int kCellsPerModule = 32;
    cellid = cellid % kCellsPerModule;
    // In vertical planes, high cell numbers have longer fibres.
    // In horizontal planes it's the opposite, so flip it round.
    // For Test Beam, vertical planes actually follow the same pattern as
    // the horizontals, so flip both.
    if( view == geo::kY || this->DetId() == novadaq::cnv::kTESTBEAM ) cellid = kCellsPerModule-cellid-1;
    // This should really never happen, but just to be safe...
    if(cellid < 0 || cellid >= kCellsPerModule) return 100;
    // Email from Tom Chase 2011-04-29
    // NB: this isn't just a ~3.8cm change per cell. Every 8 cells something
    // different happens.
    const double kPigtails[kCellsPerModule] = {
      34.5738, 38.4379, 42.3020, 46.1660, 50.0301, 53.8942, 57.7582, 61.6223,
      64.7504, 68.6144, 72.4785, 76.3426, 80.2067, 84.0707, 87.9348, 91.0790,
      95.3301, 99.1941, 103.058, 106.922, 110.786, 114.650, 118.514, 122.379,
      125.507, 129.371, 133.235, 137.099, 140.963, 144.827, 148.691, 150.751
    };
    return kPigtails[cellid];
  }

bool geo::GeometryBase::getPlaneAndCellID ( const CellUniqueId &  id, int &  plane, int &  cell  ) const

inherited

Definition at line 892 of file GeometryBase.cxx.

References geo::OfflineChan::Cell(), geo::GeometryBase::getPlaneCellMap(), makeTrainCVSamples::int, and geo::OfflineChan::Plane().

                                                                                         {
    const OfflineChan offline_chan = getPlaneCellMap(id);
    plane = (int)offline_chan.Plane();
    cell  = (int)offline_chan.Cell();
    return true;
  }

bool geo::GeometryBase::getPlaneAndCellID ( const double &  x, const double &  y, const double &  z, int &  plane, int &  cell  ) const

inlineinherited

Definition at line 108 of file GeometryBase.h.

References geo::GeometryBase::getPlaneAndCellID(), ana::NuMIBeamDirection(), POTSpillRate::view, x1, X1, submit_syst::x2, X2, y1, Y1, submit_syst::y2, Y2, Z1, and Z2.

Referenced by geo::GeometryBase::getPlaneAndCellID().

{return getPlaneAndCellID(CellId(x,y,z), plane, cell);}

OfflineChan geo::GeometryBase::getPlaneCellMap ( const CellUniqueId &  id ) const

protectedinherited

Definition at line 899 of file GeometryBase.cxx.

References geo::GeometryBase::fIdMap.

Referenced by geo::GeometryBase::getCellID(), geo::GeometryBase::getPlaneAndCellID(), geo::GeometryBase::getPlaneID(), geo::GeometryBase::IdToCell(), and geo::GeometryBase::IdToPlane().

                                                                       {

    UIDMap::const_iterator itr = fIdMap.find(id);
    if (itr == fIdMap.end()) {
      throw cet::exception("BAD_CELL_UNIQUE_ID")
        << "CellUniqueId " << id << "\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }

    return itr->second;
  }

int geo::GeometryBase::getPlaneID ( const CellUniqueId &  id ) const

inherited

Same as IdToCell and IdToPlane, but doesn't calculate CellGeo (PlaneGeo), just the global plane number or cell number

Definition at line 884 of file GeometryBase.cxx.

References geo::GeometryBase::getPlaneCellMap(), and geo::OfflineChan::Plane().

Referenced by fuzz::FuzzyKVertex::AddProng(), calib::MuondEdx::filter(), calib::MuonTrackHits::filter(), calib::MuondEdxAna::filter(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), slid::Recluster::produce(), jmshower::RecoJMShower::RecoShowers(), and slid::NuEEnergyAlg::VertexEnergy().

                                                           {
    return getPlaneCellMap(id).Plane();
  }

int geo::GeometryBase::getPlaneID ( const double &  x, const double &  y, const double &  z  ) const

inlineinherited

Definition at line 106 of file GeometryBase.h.

References geo::GeometryBase::getPlaneID().

Referenced by geo::GeometryBase::getPlaneID().

{return getPlaneID(CellId(x,y,z));}

const std::set< unsigned int > & geo::GeometryBase::GetPlanesByView ( View_t  v = kXorY ) const

inherited

Return list of planes which measure the requested projection

Parameters

v	: X/Y or All (see View_t typedef in PlaneGeo.h)

Returns

: a set of plane indices satisfying the request

Exceptions

geo::Exception if initialization fails

Definition at line 918 of file GeometryBase.cxx.

References geo::GeometryBase::fAllPlanes, geo::GeometryBase::fXplanes, geo::GeometryBase::fYplanes, geo::kX, and geo::kY.

Referenced by mcchk::CosmicAna::analyze(), htk::HoughTrack::analyze(), dprf::ChannelPlots::beginRun(), dprf::TrackPlots::beginRun(), mcchk::DetAna::beginRun(), calib::DriftResponseCalc::beginRun(), calib::CosmicTrends::beginRun(), calib::DriftResponseCalc::endRun(), rsim::RecordNoiseSpectrumFile::endSubRun(), LightLevels::endSubRun(), calib::AttenCache::RepresentativePlane(), chaninfo::BadChanList::SetRandomBadChannels(), geo::GeometryTest::testFindPlanes(), and dt::View::View().

                                                                          {
    // Return the user's choice, fall back on all planes
    if      (v==kX) return fXplanes;
    else if (v==kY) return fYplanes;
    return fAllPlanes;
  }// end of GetPlanesByView

const CellGeo * geo::GeometryBase::IdToCell ( const CellUniqueId &  id, int *  iplane, int *  icell  ) const

inherited

Given a unique cell identifier, look up the plane and cell numbers.

Parameters

id	: Cell identifier
iplane	: pointer to return plane number
icell	: pointer to return cell number

Returns

cell geometry description

Exceptions

geo::Exception

if id is not found

Definition at line 859 of file GeometryBase.cxx.

References geo::OfflineChan::Cell(), geo::PlaneGeo::Cell(), geo::GeometryBase::getPlaneCellMap(), makeTrainCVSamples::int, geo::OfflineChan::Plane(), and geo::GeometryBase::Plane().

Referenced by geo::GeometryBase::CountCellsOnLine(), geo::GeometryBase::CurrentCell(), geo::GeometryBase::GetPigtail(), photrans::ImprovedTransport::GetPosCorr(), g4n::FLSHitListAction::SteppingAction(), geo::GeometryTest::testUniqueId(), and lem::VertexToPlaneAndCell().

                                                          {
    const OfflineChan offline_chan = getPlaneCellMap(id);
    *iplane = (int)offline_chan.Plane();
    *icell  = (int)offline_chan.Cell();
    return this->Plane(*iplane)->Cell(*icell);
  }// end of IdToCell

const PlaneGeo * geo::GeometryBase::IdToPlane ( const CellUniqueId &  id, int *  iplane  ) const

inherited

Given a unique cell identifier, look up the plane number.

Parameters

id	: Cell identifier
iplane	: optional, pointer to return plane number

Returns

Plane geometry description

Exceptions

geo::Exception

if id is not found

Definition at line 879 of file GeometryBase.cxx.

References geo::GeometryBase::getPlaneCellMap(), geo::OfflineChan::Plane(), and geo::GeometryBase::Plane().

Referenced by geo::GeometryBase::GetPigtail(), and trk::KalmanGeoHelper::MatchTrajectoryToPlane().

                                                                                   {
    *iplane = getPlaneCellMap(id).Plane();
    return this->Plane(*iplane);
  }// end of IdToCell

bool geo::GeometryBase::IntersectsBigBox ( double *  xyz_cm, double *  dxyz  ) const

inherited

Definition at line 682 of file GeometryBase.cxx.

References geo::GeometryBase::fDetectorBigBoxXHi, geo::GeometryBase::fDetectorBigBoxXLo, geo::GeometryBase::fDetectorBigBoxYHi, geo::GeometryBase::fDetectorBigBoxYLo, geo::GeometryBase::fDetectorBigBoxZHi, geo::GeometryBase::fDetectorBigBoxZLo, and geo::IntersectsBox().

                                                                      {
    return geo::IntersectsBox(xyz, dxyz,
                              fDetectorBigBoxXLo, fDetectorBigBoxXHi,
                              fDetectorBigBoxYLo, fDetectorBigBoxYHi,
                              fDetectorBigBoxZLo, fDetectorBigBoxZHi);
  }// end of IntersectsBigBox

bool geo::GeometryBase::IntersectsDetector ( double *  xyz_cm, double *  dxyz  ) const

inherited

Definition at line 674 of file GeometryBase.cxx.

References geo::GeometryBase::fDetHalfHeight, geo::GeometryBase::fDetHalfWidth, geo::GeometryBase::fDetLength, and geo::IntersectsBox().

Referenced by crvtx::CosmicRayVertex::produce().

                                                                       {
    return geo::IntersectsBox(xyz, dxyz,
                              -fDetHalfWidth, fDetHalfWidth,
                              -fDetHalfHeight, fDetHalfHeight,
                              0., fDetLength);
  }// end of IntersectsDetector

bool geo::GeometryBase::IsActive ( const TGeoMaterial *  m ) const

inherited

Is the material m an active or passive material?

Definition at line 453 of file GeometryBase.cxx.

References string.

Referenced by bpfit::Path::FindLayers().

                                                         {
    const std::string scint = "Scintillator";
    if (scint == m->GetName()) return true;
    return false;
  }

bool geo::GeometryBase::isDetectorBigBoxUsed ( ) const

inlineinherited

Do we use the Detector Big Box.

Definition at line 190 of file GeometryBase.h.

References geo::ClosestApproach(), dir, dx, dy, CLHEP::m, mat, plot_validation_datamc::p1, plot_validation_datamc::p2, fillBadChanDBTables::step, w, submit_syst::x, submit_syst::y, and test::z.

Referenced by g4n::FLSHitListAction::SteppingAction().

{return fIsDetectorBigBoxUsed;}

bool geo::GeometryBase::isInsideDetectorBigBox ( const double  x_cm, const double  y_cm, const double  z_cm  ) const

inherited

Is the particle inside the detector Big Box?

Definition at line 1371 of file GeometryBase.cxx.

References geo::GeometryBase::fDetectorBigBoxXLo, geo::GeometryBase::fDetectorBigBoxYLo, and geo::GeometryBase::fDetectorBigBoxZLo.

Referenced by nuesand::FillNueSandbox::FillTruthVars(), cheat::BackTracker::InterceptsDetector(), g4n::ParticleListAction::PostTrackingAction(), rwgt::MakeGENIEReweightTable::produce(), RockFilter_A(), RockFilter_B(), RockFilter_C(), RockFilter_D(), g4n::ParticleListAction::SteppingAction(), and g4n::FLSHitListAction::SteppingAction().

                                                                                                {
    if( x < fDetectorBigBoxXHi && x > fDetectorBigBoxXLo &&
        y < fDetectorBigBoxYHi && y > fDetectorBigBoxYLo &&
        z < fDetectorBigBoxZHi && z > fDetectorBigBoxZLo)
      return true;
    return false;
  }// end of GeometryBase::isInsideDetectorBigBox

bool geo::GeometryBase::isInsideFiducialVolume ( const double  x_cm, const double  y_cm, const double  z_cm  ) const

inherited

Is the particle inside the detector Fiducial Volume?

Is the particle inside the fiducial volume?

Parameters

x_cm	: x coordinate in cm
y_cm	: y coordinate in cm
z_cm	: z coordinate in cm

Definition at line 1389 of file GeometryBase.cxx.

References geo::GeometryBase::fFiducialVolumeXLo, geo::GeometryBase::fFiducialVolumeYLo, and geo::GeometryBase::fFiducialVolumeZLo.

Referenced by geo::GeometryBase::calculateMassesLong(), filter::Filter::FillTruthVariables(), bsf::BremShowerFilter::findShowerByTruth(), geo::GeometryBase::isInsideFiducialVolume(), nuesand::FillNueSandbox::produce(), presel::NuePresel::produce(), and lem::Preselection::produce().

                                                                                                         {

    return (x_cm < fFiducialVolumeXHi && x_cm > fFiducialVolumeXLo &&
            y_cm < fFiducialVolumeYHi && y_cm > fFiducialVolumeYLo &&
            z_cm < fFiducialVolumeZHi && z_cm > fFiducialVolumeZLo);
  }// end of GeometryBase::isInsideFiducialVolume

bool geo::GeometryBase::isInsideFiducialVolume ( const double *  xyz ) const

inherited

Is the particle inside the fiducial volume?

Parameters

xyz	: pointer to array holding xyz in world coordinates

Definition at line 1402 of file GeometryBase.cxx.

References CLHEP::cm, and geo::GeometryBase::isInsideFiducialVolume().

                                                                   {

    return isInsideFiducialVolume(xyz[0]/CLHEP::cm,
                                  xyz[1]/CLHEP::cm,
                                  xyz[2]/CLHEP::cm);
  }// end of GeometryBase::isInsideFiducialVolume

bool geo::GeometryBase::LoadGeometryFile ( std::string  gdmlfile, novadaq::cnv::DetId  detId = novadaq::cnv::kUNKNOWN_DET  )

protectedinherited

In case of unknown ID, guesses from the filename.

Method to load detector geometry file

Definition at line 229 of file GeometryBase.cxx.

References geo::GeometryBase::BuildMaps(), geo::CellGeo::clear(), geo::GeometryBase::fBigBoxRange, geo::GeometryBase::fDetectorBigBoxXHi, geo::GeometryBase::fDetectorBigBoxXLo, geo::GeometryBase::fDetectorBigBoxYHi, geo::GeometryBase::fDetectorBigBoxYLo, geo::GeometryBase::fDetectorBigBoxZHi, geo::GeometryBase::fDetectorBigBoxZLo, geo::GeometryBase::fDetHalfHeight, geo::GeometryBase::fDetHalfWidth, geo::GeometryBase::fDetLength, geo::GeometryBase::fIdMap, geo::GeometryBase::FindPlanes(), geo::GeometryBase::fIsDetectorBigBoxUsed, geo::GeometryBase::fPlanes, geo::GeometryBase::fPlanesInMuonCatcher, MECModelEnuComparisons::i, novadaq::cnv::kFARDET, novadaq::cnv::kNDOS, novadaq::cnv::kNEARDET, novadaq::cnv::kTESTBEAM, novadaq::cnv::kUNKNOWN_DET, geo::plane_sort(), geo::GeometryBase::setDetectorBigBox(), geo::GeometryBase::setDetectorID(), geo::GeometryBase::SetDetectorSize(), and geo::GeometryBase::SetDrawOptions().

Referenced by geo::GeometryBase::GeometryBase(), and preBeginRun().

                                                             {

    int old_verbosity = gGeoManager->GetVerboseLevel();
    
    // TGeoManager is too verbose when loading geometry.
    // Make it quiet.
    gGeoManager->SetVerboseLevel(0);

    // Open the GDML file
    struct stat sb;
    if (gdmlfile.empty() || stat(gdmlfile.c_str(), &sb)!=0) {
      // failed to resolve the file name
      throw cet::exception("GeometryFileLoad")
        << "empty GDML string " << gdmlfile << " bail ungracefully\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }
    // clear the plane array and id map
    for (unsigned int i=0; i<fPlanes.size(); ++i) {
      if (fPlanes[i]) { delete fPlanes[i]; fPlanes[i] = 0; }
    }
    fPlanes.clear();
    fPlanesInMuonCatcher.clear();
    fIdMap.clear();
    CellGeo::clear();

    TGeoManager::Import(gdmlfile.c_str());

    this->SetDrawOptions();

    std::vector<const TGeoNode*> n(16);
    n[0] = gGeoManager->GetTopNode();
    this->FindPlanes(n, 0);
    sort(fPlanes.begin(), fPlanes.end(), plane_sort);
    sort(fPlanesInMuonCatcher.begin(), fPlanesInMuonCatcher.end(), plane_sort);

    this->SetDetectorSize();
    this->BuildMaps();

    if(detId == novadaq::cnv::kUNKNOWN_DET){
      // have to set the fDetId based on the input file then
      if(     gdmlfile.find("ndos") != std::string::npos) detId = novadaq::cnv::kNDOS;
      else if(gdmlfile.find("near") != std::string::npos) detId = novadaq::cnv::kNEARDET;
      else if(gdmlfile.find("far")  != std::string::npos) detId = novadaq::cnv::kFARDET;
      else if(gdmlfile.find("testbeam-1x1-2block") != std::string::npos) detId = novadaq::cnv::kTESTBEAM;
      else if(gdmlfile.find("testbeam-1x1-3block") != std::string::npos) detId = novadaq::cnv::kTESTBEAM;
      else if(gdmlfile.find("testbeam-2x2-2block") != std::string::npos) detId = novadaq::cnv::kTESTBEAM;
      else if(gdmlfile.find("testbeam-2x2-3block") != std::string::npos) detId = novadaq::cnv::kTESTBEAM;
    }
    setDetectorID(detId);

    if(fIsDetectorBigBoxUsed)
      this->setDetectorBigBox(fBigBoxRange);
    else{
      fDetectorBigBoxXHi =  fDetHalfWidth;
      fDetectorBigBoxXLo = -fDetHalfWidth;
      fDetectorBigBoxYHi =  fDetHalfHeight;
      fDetectorBigBoxYLo = -fDetHalfHeight;
      fDetectorBigBoxZHi =  fDetLength;
      fDetectorBigBoxZLo =  0.;
    }

    mf::LogInfo("LoadGeometryFile")
      << "Detector geometry loaded from " << gdmlfile;

    
    // Restore verbosity level before we leave because
    // we have to be polite and listen
    gGeoManager->SetVerboseLevel(old_verbosity);
    
    return true;
  }// end of LoadGeometryFile

void geo::GeometryBase::MakePlane ( std::vector< const TGeoNode * > &  n, unsigned int  depth, bool  inMuonCatcher  )

protectedinherited

Definition at line 731 of file GeometryBase.cxx.

References geo::GeometryBase::fPlanes, and geo::GeometryBase::fPlanesInMuonCatcher.

Referenced by geo::GeometryBase::FindPlanes().

  {
    fPlanes.push_back(new PlaneGeo(n, depth));

    if(inMuonCatcher) fPlanesInMuonCatcher.push_back(fPlanes.back());
  }

int geo::GeometryBase::MakeTmpFile ( std::string  gdmlInfo )

protectedinherited

Method to make a temporary file holding Detector geometry (gdml) file

Definition at line 102 of file GeometryBase.cxx.

References om::cerr, allTimeWatchdog::endl, fclose(), geo::GeometryBase::fGDMLFile, geo::GeometryBase::fParams, geo::GeometryBase::Params::StoreTempGeo, string, and tmp.

Referenced by preBeginRun().

  {
    int fStat = -1; //file status variable, -1 indicates error
    FILE *tmpFile; //file variable
    
    //hold temporary gdml file name
    std::string tmpFileName;
    
    //Use mkstemp to make a temporary file, mkstemp takes in a template
    //character array as an argument.  The last 6 characters must be XXXXXX
    //mkstemp will then attempt to open a file with that template and the 
    //last 6 characters will be changed to be unique.  If the file cannot
    //be successfully opened an error code will be returned.
    //If a file is properly opened it will be unique.  File is set to be
    //opened in /temp.
    char tmp[256];
    sprintf(tmp,"%s/TEMP_GEO_GDML_XXXXXX", fParams.StoreTempGeo().c_str());
    //check to see if file opens properly
    if ( ((fStat = mkstemp(tmp)) == -1) ||
           ((tmpFile = fopen(tmp, "w+")) == NULL) ) {
      
      if (fStat != -1){
        fclose(tmpFile);
        remove(tmp);
        std::cerr<<"Error making temporary file "
            <<tmp<<std::endl;
      }
      //if file could not be created properly exit and return -1
      std::cerr<<"Failed to produce temporary geometry file"<<std::endl;
      return -1;
    }
    
    //this means a unique file has been properly created and opened
    fclose(tmpFile); //properly close file
    
    //add temp file name to string with proper .gdml extension
    tmpFileName += tmp;
    tmpFileName += ".gdml";
    
    //now attempt to open new file with .gdml extension
    FILE *gdmlFile;
    gdmlFile = fopen(tmpFileName.c_str(), "w+");
    if (gdmlFile == NULL){
      //file was not opened properly, remove temp files
      //and exit with a status -1
      remove(tmp);
      remove(tmpFileName.c_str());
      std::cerr<<"Error making temporary file "
               <<tmpFileName<<std::endl;
      return -1;
    }
    
    //temporary gdml file was properly opened, now put our gdml info in file
    fputs (gdmlInfo.c_str(), gdmlFile);
    
    //now remove initial temp file without the .gdml extension
    remove(tmp);
    
    //close temporary gdml file
    fclose(gdmlFile);
    
    //now store name of temporary gdml file in GDMLFile
    //this way the temp file can also be accessed by GEANT4,
    //which loads the geometry directly
    fGDMLFile = tmpFileName;
    
    //return 0 indicating successful temporary file creation
    return 0;
  }

double geo::GeometryBase::MassBetweenPoints ( double *  p1, double *  p2  ) const

inherited

Return the column density between 2 points in gm/cm^2

Parameters

p1	: pointer to array holding xyz of first point in world coordinates
p2	: pointer to array holding xyz of second point in world coordinates

Definition at line 1155 of file GeometryBase.cxx.

References E, demo0::length, geo::GeometryBase::MaterialsBetweenPoints(), file_size_ana::node, and util::pythag().

  {
    // The purpose of this method is to determine the column density
    // between the two points given.  Do that by starting at p1 and
    // stepping until you get to the node of p2.  calculate the distance
    // between the point just inside that node and p2 to get the last
    // bit of column density
    double columnD = 0.;

    // first initialize a track - get the direction cosines
    double length = util::pythag(p2[0]-p1[0],
                                 p2[1]-p1[1],
                                 p2[2]-p1[2]);
    double dxyz[3] = {(p2[0]-p1[0])/length,
                      (p2[1]-p1[1])/length,
                      (p2[2]-p1[2])/length};

    gGeoManager->InitTrack(p1,dxyz);
    gGeoManager->FindNextBoundary(1.E-6);
    
    // might be helpful to have a point to a TGeoNode
    TGeoNode *node = gGeoManager->GetCurrentNode();

    // check that the points are not in the same volume already.
    // if they are in different volumes, keep stepping until you
    // are in the same volume as the second point
    while(!gGeoManager->IsSameLocation(p2[0], p2[1], p2[2])){
      gGeoManager->FindNextBoundary();

      columnD +=
        gGeoManager->GetStep() *
        node->GetMedium()->GetMaterial()->GetDensity();

      // the act of stepping puts you in the next node and returns that
      // node
      node = gGeoManager->Step();
    } //end loop to get to volume of second point

    // now you are in the same volume as the last point, but not at that
    // point.  get the distance between the current point and the last
    // one
    const double *current = gGeoManager->GetCurrentPoint();
    length = util::pythag(p2[0]-current[0],
                          p2[1]-current[1],
                          p2[2]-current[2]);
    columnD += length*node->GetMedium()->GetMaterial()->GetDensity();

    return columnD;
  }// end of MassBetweenPoints

TGeoMaterial * geo::GeometryBase::Material ( double  x, double  y, double  z  ) const

inherited

Method to obtain current detector material at given (x,y,z) coordinates

Definition at line 444 of file GeometryBase.cxx.

References file_size_ana::node.

Referenced by geo::GeometryBase::calculateMassesLong(), and calib::EnergyStandardCandles::FillTrueMichelInfo().

                                                                         {
    const TGeoNode* node = gGeoManager->FindNode(x,y,z);
    return node->GetMedium()->GetMaterial();
  }// end of Material

void geo::GeometryBase::MaterialsBetweenPoints ( const double *  p1, const double *  p2, std::vector< double > &  ds, std::vector< const TGeoMaterial * > &  mat  ) const

inherited

Definition at line 1208 of file GeometryBase.cxx.

References om::cerr, d, E, e, allTimeWatchdog::endl, file_size_ana::node, make_associated_cosmic_defs::p3, util::pythag(), and fillBadChanDBTables::step.

Referenced by calib::EnergyStandardCandles::FillTrueMichelInfo(), bpfit::Path::FindLayers(), bpfit::Path::IntegrateLeg(), geo::GeometryBase::MassBetweenPoints(), geo::GeometryBase::MaterialsBetweenPoints(), and geo::LiveGeometry::ProjectedSteelDist().

  {
    double d[3], dmag;
    dmag = util::pythag(p2[0]-p1[0],p2[1]-p1[1],p2[2]-p1[2]);
    d[0] = (p2[0]-p1[0])/dmag;
    d[1] = (p2[1]-p1[1])/dmag;
    d[2] = (p2[2]-p1[2])/dmag;
    
    gGeoManager->InitTrack(p1,d);
    gGeoManager->FindNextBoundary(1.E-6);
    TGeoNode *node = gGeoManager->GetCurrentNode();

    while (!gGeoManager->IsSameLocation(p2[0], p2[1], p2[2])) {
      gGeoManager->FindNextBoundary();
      
      // If we take a step and haven't gone anywhere, then we will be
      // stuck in a very long loop. So we will quit if the next step
      // takes us nowhere.
      double step = gGeoManager->GetStep();
      if(step < 1.0e-12) break;

      ds.push_back(step);
      m.push_back(node->GetMedium()->GetMaterial());
      
      node = gGeoManager->Step();
      if (node==0) {
        std::cerr << __FILE__ << ":" << __LINE__ 
                  << " Cannot continue tracking MaterialsBetweenPoints(). "
                  << "Null node pointer." << std::endl;
        return;
      }
    }
    //
    // Handle final step by hand
    //
    const double* p3 = gGeoManager->GetCurrentPoint();
    ds.push_back(util::pythag(p3[0]-p2[0],p3[1]-p2[1],p3[2]-p2[2]));
    m.push_back(node->GetMedium()->GetMaterial());
  }

void geo::GeometryBase::MaterialsBetweenPoints ( TVector3  v1, TVector3  v2, std::vector< double > &  ds, std::vector< const TGeoMaterial * > &  m  ) const

inherited

Overloaded MaterialsBetweenPoints to for using TVector3 instead of arrays because rb::Track trajectory points are TVector3 and there is a desire to track a 'track' through the geometry. This will make that easier. Function is really just a wrapper around the array version

Definition at line 1258 of file GeometryBase.cxx.

References geo::GeometryBase::MaterialsBetweenPoints().

  {
    double p1[3], p2[3];
    v1.GetXYZ(p1);
    v2.GetXYZ(p2);
    
    MaterialsBetweenPoints(p1, p2, ds, m);
  }

const unsigned int geo::GeometryBase::NextPlaneInView ( unsigned int  p1, int  d = +1  ) const

inherited

Return the index of the next plane in a particular view

Parameters

p1	: Index of current plane
d	: Direction to step in (>0 next down stream, <0 upstream)

Returns

Index of next plane in same view. If no plane can be found (eg. reach end of detector) return kPLANE_NOT_FOUND.

Exceptions

Asserts that p1 is in valid range

Definition at line 936 of file GeometryBase.cxx.

References ana::assert(), geo::GeometryBase::fPlanes, geo::kPLANE_NOT_FOUND, and plot_validation_datamc::p2.

Referenced by chaninfo::CosmicEff::analyze(), trk::KalmanTrack::FilterTracker(), calib::HasZAdjacents(), and geo::GeometryTest::testFindPlanes().

                                                                               {
    assert(p1<fPlanes.size());

    // March along from p1 in the direction d until we find a match
    int s = (d>0 ? 1 : -1);
    if (s<0 && p1==0)                return kPLANE_NOT_FOUND;
    if (s>0 && p1==fPlanes.size()-1) return kPLANE_NOT_FOUND;
    for (unsigned int p2=p1+s; 1; p2+=s) {
      if (fPlanes[p1]->View() == fPlanes[p2]->View()) return p2;
      if (p2==fPlanes.size()-1) return kPLANE_NOT_FOUND;
      if (p2==0)                return kPLANE_NOT_FOUND;
    }
  }// end of NextPlaneInView

const unsigned int geo::GeometryBase::NextPlaneOtherView ( unsigned int  p1, int  d = +1  ) const

inherited

Return the index of the next plane in a particular view

Parameters

p1	: Index of current plane
d	: Direction to step in (>0 next down stream, <0 upstream)

Returns

index of next plane in same view.

Exceptions

Asserts that p1 is in valid range

Definition at line 960 of file GeometryBase.cxx.

References ana::assert(), geo::GeometryBase::fPlanes, geo::kPLANE_NOT_FOUND, and plot_validation_datamc::p2.

Referenced by fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), lem::FindLEMMatches::MatchToVertex(), lem::LEM::MatchToVertex(), trk::KalmanGeoHelper::MatchTrajectoryToPlaneInView(), calib::AttenCache::RepresentativePlane(), geo::GeometryTest::testFindPlanes(), and lem::VertexToPlaneAndCell().

  {
    assert(p1<fPlanes.size());

    // March along from p1 in the direction d until we find a match
    int s = (d>0 ? 1 : -1);
    if (s<0 && p1==0)                return kPLANE_NOT_FOUND;
    if (s>0 && p1==fPlanes.size()-1) return kPLANE_NOT_FOUND;
    for (unsigned int p2=p1+s; 1; p2+=s) {
      if (fPlanes[p1]->View() != fPlanes[p2]->View()) return p2;
      if (p2==fPlanes.size()) return kPLANE_NOT_FOUND;
      if (p2==0)              return kPLANE_NOT_FOUND;
    }
  }// end of NextPlaneOtherView

unsigned int geo::GeometryBase::NPlanes ( ) const

inherited

Definition at line 460 of file GeometryBase.cxx.

References ana::assert(), geo::GeometryBase::fDetId, geo::GeometryBase::fPlanes, and novadaq::cnv::kUNKNOWN_DET.

Referenced by align::Alignment::analyze(), chaninfo::BadChanValidate::BadChansFromDB(), chaninfo::CosmicEff::beginRun(), dprf::ChannelPlots::beginRun(), dprf::TrackPlots::beginRun(), calib::HitEfficiency::beginRun(), rsim::ReadoutSim::beginRun(), calib::CosmicTrends::beginRun(), chaninfo::BadChanValidate::beginSubRun(), noe::build_cell_lookup_table(), chaninfo::BadChanList::CacheBadChannels(), geo::GeometryBase::CellInfo(), geo::GeometryBase::CellTpos(), trk::KalmanGeoHelper::CountMissedCellsOnLine(), LightLevels::CreateHistos(), rsim::ReadoutSim::CreateRawDigits(), lem::DefaultVertex(), evd::GeometryDrawer::DrawCells2D(), chaninfo::CosmicEff::endJob(), calib::AttenFit::endJob(), comi::ChanOcc::endSubRun(), rsim::IFPGAAlgorithm::FetchThresholds(), trk::KalmanGeoHelper::fillprivategeo(), trk::CosmicTrackUtilities::FindZBoundaries(), calib::FindZBoundaries(), evd::GeometryDrawer::GetBox(), evd::GeometryDrawer::GetDCMBoxes(), airshower::AirSlicer::GetHoughMap(), airshower::AirSlicer::GetHoughRhos(), calib::GetXYZD(), dt::DiscreteTracker::HighActivity(), calib::RockMuonStopperSelection::IsContainedStopper(), photrans::ImprovedTransport::LoadHistos(), chaninfo::BadChanList::LoadNewData(), lem::FindLEMMatches::MatchToVertex(), lem::LEM::MatchToVertex(), skim::ParametersNumu::ParametersNumu(), trk::KalmanTrack::produce(), lem::AlignLibToVtx::produce(), cvn::RegCVNMapper::produce(), rsim::ReadoutSim::produce(), evd::RawDataDrawer::RawDigit2D(), lem::GenFromLib::readNext(), trk::KalmanGeoHelper::SetupHelper(), geo::GeometryTest::testCellId(), geo::GeometryTest::testCellIdFromPos(), dt::View::View(), and calib::AttenuationFit::writeResults().

  {
    assert(fDetId != novadaq::cnv::kUNKNOWN_DET);
    return fPlanes.size();
  }

TVector3 geo::GeometryBase::NuMIBeamDirection ( ) const

inherited

Direction of neutrinos from the NuMI beam (unit vector)

Because of the varying location of pion decays in the decay pipe, this can only be an average value in the ND and NDOS.

Definition at line 523 of file GeometryBase.cxx.

References e, geo::GeometryBase::fDetId, novadaq::cnv::kFARDET, novadaq::cnv::kNDOS, novadaq::cnv::kNEARDET, and novadaq::cnv::kTESTBEAM.

Referenced by cvn::CVNEventProngDump::analyze(), comi::NearlineAna::analyze(), slid::LIDTraining::analyze(), cosrej::CosRejFxs::getAngle(), numusand::NumuSandFxs::getAngle(), nuesand::FillNueSandbox::GetTrackAngle(), numue::NumuEAlg::NumuEAlg(), nussand::FillNusSandbox::produce(), slid::SliceLIDBuilder::produce(), cosrej::CosRej::produce(), comi::Cana::produce(), slid::LIDBuilder::produce(), jmshower::NueSel::produce(), slid::SPIDBuilder::produce(), jmshower::NueSelLID::produce(), cvn::CVNAddTrainingData::produce(), cosrej::MakeNueCosRej::produce(), jmshower::RecoJMShower::RecoShowers(), ncid::MakeNCCosRej::TransMomFraction(), ncid::NCNNKeras::TransMomFraction(), and cosrej::MakeNueCosRej::TransMomFraction().

  {
    // These are obtained by plotting rec.mc.nu.p.fP.fZ/rec.mc.nu.p.fE etc from
    // CAF files and fitting the resulting distribution to a gaussian.
    // A true vertex inside of the detector is require as well as the NuE group FA reconstruction
    // and containment cuts.
    // CAF Files used:
    //   NearDet:
    //    prod_caf_S15-05-22a_nd_genie_fhc_nonswap_ndnewpos
    //   FarDet:
    //    prod_caf_S15-05-22a_fd_genie_fhc_nonswap_fdfluxv08 
    //   Unknown for NDOS

    switch(fDetId){
    case novadaq::cnv::kFARDET:   return TVector3(-6.833e-05,
                                                 +6.388e-02,
                                                 +9.980e-1).Unit();
    case novadaq::cnv::kNDOS:     return TVector3(-3.845e-4,
                                                 +6.517e-2,
                                                 +9.967e-1).Unit();
    case novadaq::cnv::kNEARDET:  return TVector3(-8.424e-04,
                                                 -6.174e-02,
                                                 +9.981e-1).Unit();
    case novadaq::cnv::kTESTBEAM: return TVector3(0,
                                                  0,
                                                  1).Unit();
    default:
      throw cet::exception("BadGeoConfig")
        << "Bad geometry configuration: " << fDetId << "\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }
  }

const PlaneGeo * geo::GeometryBase::Plane ( unsigned int  i ) const

inherited

Return the geometry description of the ith plane in the detector.

Parameters

i	: input plane number, starting from 0

Returns

plane geometry for ith plane

Exceptions

geo::Exception

if iplane is outside allowed range

Definition at line 749 of file GeometryBase.cxx.

References geo::GeometryBase::fPlanes, and MECModelEnuComparisons::i.

Referenced by dif::DiFShowerFinder::adjustPlane(), dif::DiFShowerFinder::adjustPlane_end(), geo::GeometryTest::analyze(), mcchk::DetAna::analyze(), dprf::ChannelPlots::analyze(), zcl::SMMTriggerAna::analyze(), dprf::TrackPlots::analyze(), slicer::S4DParamCalc::analyze(), mcchk::CosmicAna::analyze(), mcchk::RockAna::analyze(), numue::NumuEAna::analyze(), align::Alignment::analyze(), sn::SNSlicerAna::analyze(), calib::AssessCalib::analyze(), calib::HitEfficiency::analyze(), mcchk::LeptonAna::analyze(), comi::NumiFilteringAna::analyze(), zcl::FmmTrackerAna::analyze(), zcl::FmmTrackerValidation::analyze(), trk::KalmanTrackAna::analyze(), showere::ShowerEnergyAna::analyze(), remid::ReMIdDedxRock::analyze(), remid::ReMIdDedxFD::analyze(), remid::ReMIdDedxStudies::analyze(), LightLevels::analyze(), calib::CosmicTrends::analyze(), htk::HoughTrack::analyze(), trk::CosmicTrackAna::analyze(), comi::Leana::AnaNHit(), calib::CalibAna::ApproximateHitPos(), chaninfo::BadChanValidate::BadChansFromDB(), trk::CosmicTrackAna::beginJob(), chaninfo::CosmicEff::beginRun(), dprf::ChannelPlots::beginRun(), mcchk::DetAna::beginRun(), dprf::TrackPlots::beginRun(), calib::DriftResponseCalc::beginRun(), rsim::ReadoutSim::beginRun(), calib::CosmicTrends::beginRun(), chaninfo::BadChanValidate::beginSubRun(), calib::BestPathEstimates(), noe::build_cell_lookup_table(), geo::GeometryBase::BuildMaps(), chaninfo::BadChanList::CacheBadChannels(), slid::ParticleIDAlg::CalcAsymIneria(), slid::ParticleIDAlg::CalcCellPlaneTransverseDedx(), slid::ParticleIDAlg::CalcPlaneTransverseDedx(), slid::ParticleIDAlg::CalcTrkHitPath(), trk::KalmanTrackMerge::CanJoinTracks(), evd::RecoBaseDrawer::CellHit2D(), evd::RecoBaseDrawer::CellHit3D(), geo::GeometryBase::CellInfo(), geo::GeometryBase::CellTpos(), murem::TrackCleanUpAlg::CleanUpTrack(), geo::GeometryBase::ClosestApproach(), dt::Cand::ClosestToEndCell(), sn::SNMichelAnalyzer::ClusterAroundTrackEnds(), remid::ReMIdTrain::ContainedEvent(), remid::ReMIdDedx::ContainedEvent(), remid::ReMIdDedxRock::ContainedEvent(), remid::ReMIdDedxStudies::ContainedEvent(), skim::SkimmingUtils::CopyProngs(), skim::SkimmingUtils::CopyShowers(), skim::SkimmingUtils::CopySlice(), skim::SkimmingUtils::CopyTracks(), geo::GeometryBase::CountCellsOnLine(), trk::KalmanGeoHelper::CountMissedCellsOnLine(), util::CountXY(), LightLevels::CreateHistos(), rsim::ReadoutSim::CreateRawDigits(), murem::TrackCleanUpAlg::DeDxInPlane(), lem::DefaultVertex(), dif::DiFShowerFinder::distancefromtrack(), sn::TrackActivityRemover::DistanceHitToPoint(), sn::SNMichelAnalyzer::DistanceHitToPoint(), sn::TrackActivityRemover::DistanceHitToTrackBody(), slid::DeconvolveAlg::DistanceToCore(), evd::GeometryDrawer::DrawCells2D(), evd::RecoBaseDrawer::DrawCluster2D(), evd::GeometryDrawer::DrawHighlightCell(), evd::RecoBaseDrawer::DrawTrack2D(), skim::ParametersNumu::EarliestLatestHitPos(), mcchk::ShowerAnaCheck::EMMoliereRadius(), mcchk::ShowerAnaCheck::EMShowerWidth(), chaninfo::CosmicEff::endJob(), calib::BetheBlochFit::endJob(), calib::AttenFit::endJob(), calib::FindOutliers::endRun(), comi::ChanOcc::endSubRun(), rsim::RecordNoiseSpectrumFile::endSubRun(), LightLevels::endSubRun(), dif::DiFShowerFinder::eparm(), calib::MakeAttenuationProfiles::event(), util::EventBox(), calib::AttenFit::expFit(), calib::AttenuationFit::expFit(), rsim::IFPGAAlgorithm::FetchThresholds(), geo::LiveGeometry::FillBadBoxes(), caf::FillHadClustVars(), calib::StopperThreshold::FillHist(), calib::CosmicTrends::FillHistograms(), bpfit::BreakPoint::FillHitList(), trk::KalmanGeoHelper::fillprivategeo(), caf::FillSliceInfo(), caf::FillSliceVars(), trk::CosmicTrackAna::FillTrackHistograms(), calib::StopperThreshold::FillTree(), calib::CosmicCalib::filter(), comi::NumiFiltering::filter(), bsf::BremShowerFilter::filter(), showere::ShowerEnergyFilterMC::filter(), calib::FindBelowThresholdCalibCandidates(), trk::WindowTrackingAlg::FindEndPoint(), bsf::BremShowerFilter::findShowerByReco(), dif::DiFShowerFinder::findShowerByReco(), trk::KalmanTrack::FindTracks(), trident::DimuonSliceFilter::FindVertex(), bpfit::DimuonFitter::FindVertexZ(), trk::CosmicTrackUtilities::FindZBoundaries(), calib::FindZBoundaries(), bpfit::BreakPoint::FitTracks(), bpfit::DimuonFitter::FitView(), trk::WindowTrackingAlg::FitWindow(), evd::SimulationDrawer::FLSHit2D(), evd::SimulationDrawer::FLSHit3D(), align::SplitTracks::GeoEstimatePos(), mono::Monopole::get_position(), mono::SlowMonopoleAna::get_position(), numusand::NumuSandFxs::getActivity(), calib::Calibrator::GetAttenCurve(), cosrej::CosRejFxs::getBBC(), evd::SliceNavigator::GetBounds(), evd::GeometryDrawer::GetBox(), evd::SliceNavigator::GetBox(), jmshower::RecoJMShower::GetCellDistToPoint(), jmshower::RecoJMShower::GetCellDistToTrk(), align::SplitTracks::GetCellEndpoints(), slid::ParticleIDAlg::GetCellNodePos(), jmshower::RecoJMShower::GetCellNodePos(), jmshower::RecoJMShower::GetCellTransDedx(), jmshower::RecoJMShower::GetCentroid(), calib::CosmicCalib::GetChannelHists(), calib::MakeAttenuationProfiles::GetChannelHists(), evd::RecoBaseDrawer::GetClusterOutlinesConvexHull(), evd::GeometryDrawer::GetDCMBoxes(), nuesand::FillNueSandbox::GetdEdx(), geo::LiveGeometry::GetDetectorEdges(), skim::ParametersNumu::GetDetectorEdges(), cosrej::CosRejFxs::getFits(), airshower::AirSlicer::GetHitPos(), ndreco::NDRecoFxs::getProngActivity(), ndreco::NDRecoFxs::getProngDedx(), dt::Chunk::GetSegs(), ndreco::NDRecoFxs::getTrackActivity(), ndreco::NDRecoFxs::getTrackDedx(), jmshower::JMShower::GetTransHits(), jmshower::RecoJMShower::GetTrkCPlaneCell(), jmshower::RecoJMShower::GetTrkHitPath(), jmshower::RecoJMShower::GetTrkHitPos(), jmshower::RecoJMShower::GetTrkPlaneCell(), jmshower::RecoJMShower::GetTrkPlaneDistToEdge(), jmshower::RecoJMShower::GetTrkPlanePos(), calib::ThresholdAna::GetView(), calib::CosmicCalib::GetView(), calib::MakeAttenuationProfiles::GetView(), calib::AttenFit::GetView(), calib::AttenuationFit::GetView(), calib::AttenCache::GetView(), calib::GetXYZD(), slid::NuEEnergyAlg::HadronicDepEnergy(), calib::HasXYAdjacents(), dt::DiscreteTracker::HighActivity(), rb::HitList::HitList(), dt::Chunk::HitsOnLine(), nerd::NERDProng::HitToHitDistance(), me::MEFinder::HitToHitDistance(), geo::GeometryBase::IdToCell(), geo::GeometryBase::IdToPlane(), bsf::BremShowerFilter::inFiducial(), dif::DiFShowerFinder::inFiducial(), remid::RecoMuon::IsTrackContained(), trk::CosmicTrackAlg::IsTrackDownstreamFromTiming(), skim::EvaluatorNue::KeepSlice(), calib::AttenCache::LoadFromDatabase(), photrans::ImprovedTransport::LoadHistos(), chaninfo::BadChanList::LoadNewData(), trk::WindowTrackingAlg::LookForBremsstrahlungHits(), trk::WindowTrackingAlg::Make3DTrack(), calib::Calibrator::MakeCellHit(), evtsum::EventSummary::MakeOutput(), trk::CosmicTrackAlg::MakeTrack(), trk::WindowTrackingAlg::MakeTrack(), trk::WindowTrackingAlg::MakeViewTrack(), calib::MakeZBoundaryMap(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), lem::FindLEMMatches::MatchToVertex(), lem::LEM::MatchToVertex(), trk::KalmanGeoHelper::MatchTrajectoryToPlane(), trk::KalmanGeoHelper::MatchTrajectoryToPlaneInView(), rb::Cluster::MinMaxMeanXYZ(), zcl::SMMCluster::MissC(), dif::DiFShowerFinder::muonstub(), numue::NumuEAlg::NumuEAlg(), evd::RecoBaseDrawer::OfflineChans2D(), skim::ParametersNumu::ParametersNumu(), trk::CosmicTrackUtilities::PathLengthInCell(), calib::PathLengthInCell(), calib::CosmicCalib::PCHitPE(), calib::MakeAttenuationProfiles::PCHitPE(), rb::Track::PlaneDirMap(), slid::ParticleIDAlg::PlaneHitCell(), slid::ParticleIDAlg::PlaneHitXYZ(), slid::ParticleIDAlg::PlaneRadius(), ddthelpers::print(), calib::PCHitsList::ProcessTrackForBelowThresholdHits(), calib::PCHitsList::ProcessTrackForFLSHits(), calib::PCHitsList::ProcessTrackTrajectory(), photrans::PhotonSystRewriter::produce(), hv::HoughVertex::produce(), trk::KalmanTrack::produce(), trk::KalmanTrackMerge::produce(), slid::Recluster::produce(), jmshower::JMTrackMerge::produce(), numusand::FillSandbox::produce(), lem::AlignLibToVtx::produce(), lem::MakeLibrary::produce(), zcl::SMMCluster::produce(), photrans::PhotonTransport::produce(), vdt::VertexDT::produce(), jmshower::NueSel::produce(), cvn::RegCVNMapper::produce(), jmshower::NueSelLID::produce(), cheat::MCCheater::produce(), remid::RecoMuon::produce(), rsim::ReadoutSim::produce(), caf::CAFMaker::produce(), slid::ParticleIDAlg::Radius(), evd::RawDataDrawer::RawDigit2D(), lem::GenFromLib::readNext(), jmshower::RecoJMShower::RecoShowers(), calib::AttenCache::RepresentativePlane(), trk::RLFit::RLFit(), hough::MultiHoughT::Scrub(), earms::ElasticArmsHS::Scrub(), vdt::VertexDT::Scrub(), fuzz::FuzzyKVertex::Scrub(), geo::GeometryBase::SetDetectorSize(), chaninfo::BadChanList::SetRandomBadChannels(), trk::WindowTrackingAlg::SetTrackEndPoints(), trk::KalmanTrackMerge::ShiftInterpolationPoints(), trk::WindowTrackingAlg::ShortTrack(), trk::WindowTrackingAlg::ShortTrackExtraPlane(), trk::WindowTrackingAlg::ShortViewTrack(), skim::SkimmingUtils::SkimmingUtils(), lem::MakeLibrary::SliceMeanPosEWeighted(), murem::TrackCleanUpAlg::SortByDistFromTrack(), photrans::ImprovedTransport::StepAlongHit(), geo::GeometryTest::testCellId(), geo::GeometryTest::testCellIdFromPos(), geo::GeometryTest::testFindCell(), geo::GeometryTest::testFindPlanes(), calib::StopperThreshold::testPath(), geo::GeometryTest::testStepping(), geo::GeometryTest::testUniqueId(), rb::Prong::TotalLength(), dt::Chain::ToTrack(), zcl::FMMTracker::Tracking2D(), trk::CosmicTrackUtilities::TrackResiduals(), calib::AttenCache::TranslateChannelMC(), trk::CosmicTrackUtilities::TruePathLengthInCell(), calib::PCHitsList::TrueW(), slid::NuEEnergyAlg::VertexEnergy(), lem::VertexToPlaneAndCell(), dt::View::View(), trk::KalmanTrackMerge::ViewMergeTracks(), rb::Prong::W(), rb::Track::W(), calib::AttenuationFit::writeResults(), and calib::ZBounds().

                                                            {
    if (i>=fPlanes.size()) {
      throw cet::exception("BAD_PLANE_NUMBER")
        << "Plane index " << i << " out of range " << fPlanes.size() << "\n"
        << __FILE__ << ":" << __LINE__ << "\n";
      return 0;
    }
    return fPlanes[i];
  }//end of Plane

void geo::Geometry::preBeginRun

(

art::Run const &

run

)

Definition at line 59 of file Geometry_service.cc.

References om::cout, geo::GeometryBase::DetId(), allTimeWatchdog::endl, geo::GeometryBase::fGDMLFile, geo::GeometryBase::fGDMLFromFCL, geo::GeometryBase::FindGDMLFile(), fLoadedGDMLInfo, geo::GeometryBase::Params::ForceUseFCLOnly, geo::GeometryBase::fParams, geo::GeometryBase::fROOTFile, geo::GeometryBase::fRunNumber, art::DataViewImpl::getManyByType(), geo::GeometryBase::LoadGeometryFile(), geo::GeometryBase::MakeTmpFile(), geo::GeometryBase::RemoveTmpFile(), art::DataViewImpl::run(), and string.

Referenced by Geometry().

  {
    // check the current run number, do nothing if it is the same

    if(run.run() == fRunNumber) return;
    fRunNumber = run.run();

     // check here to see if we need to load a new geometry.
      // get the detector id from the run object
    std::vector< art::Handle<sumdata::RunData> > rdcol;
    run.getManyByType(rdcol);
    if(!rdcol.empty() && !fParams.ForceUseFCLOnly()){
      const std::string newGDMLInfo = rdcol[0]->DetGDMLInfo();

      if(newGDMLInfo == fLoadedGDMLInfo){
        std::cout << "Geometry for run " << fRunNumber
                  << " unchanged from previous run." << std::endl;
        // So there's nothing left to do
        return;
      }

      //if there is a gdml string in the file, use that
      if(!newGDMLInfo.empty()){

        //if geometry for last run came from a temp file, remove it
        if (fGDMLFile.find("TEMP_GEO_GDML") != std::string::npos){
          this->RemoveTmpFile(fGDMLFile);
        }

        //gdml info is nonzero, attempt to load geometry from here
        if (this->MakeTmpFile(newGDMLInfo) == 0){
          this->LoadGeometryFile(fGDMLFile, rdcol[0]->DetId());
          fLoadedGDMLInfo = newGDMLInfo;
          std::cout<<"Successfully loaded geometry from in file "<<fGDMLFile<<std::endl;
          std::cout<<"The geometry stored in the file was originally named: "<<rdcol[0]->DetGeometryBaseName()<<std::endl;
          fROOTFile = fGDMLFile; //hook left for root file
          return;
        }
        else {
          std::cout<<"WARNING: Could not load geometry information "
              <<"from within the file. The temporary file was not "
              <<"successfully created, something is wrong.  Will "
              <<"attempt to load from the gdml file directly."
              <<std::endl;
        }
      }

      // get strings with the relative paths for the detector geometry files
      std::string relpathgdml("Geometry/gdml/");
      relpathgdml += rdcol[0]->DetGeometryBaseName();
      relpathgdml += ".gdml";

      if(fGDMLFile.find(relpathgdml) != std::string::npos) return;

      // constructor decides if initialized value is a path or an environment variable
      fGDMLFile = FindGDMLFile(relpathgdml);

      mf::LogWarning("LoadNewGeometry") << "loading new geometry files\n"
                                        << fGDMLFile << "\n";
  
      this->LoadGeometryFile(fGDMLFile);
      fROOTFile = fGDMLFile; //hook left for root file  
      return;
    }
    else{
      //if we made it this far without a load, check to see if the original
      //gdml file specified in the fcl file is still stored in fGDMLFile
      //if this is not the case, make it so and load the geometry
      if (fGDMLFile.find(fGDMLFromFCL) == std::string::npos){
        fGDMLFile = fGDMLFromFCL;
        this->LoadGeometryFile(fGDMLFile);
      }
      mf::LogWarning("LoadNewGeometry") << "cannot find sumdata::RunData object to grab detid\n"
                                        << "or user has determined to only use the value "
                                        << "from the fcl file\n"
                                        << "using currently configured geometry:\n"
                                        << "\t" << fGDMLFile << "\n";
      fROOTFile = fGDMLFile; //hook left for root file
      return;
    }
    
  }// end of preBeginRun

void geo::GeometryBase::RemoveTmpFile ( std::string  fileName )

protectedinherited

Method to remove temporary file that holds detector geometry file.

Definition at line 175 of file GeometryBase.cxx.

Referenced by preBeginRun(), ~Geometry(), and geo::GeometryBase::~GeometryBase().

                                                    {
    remove(fileName.c_str());
  }

std::string geo::GeometryBase::ROOTFile ( ) const

inlineinherited

Definition at line 236 of file GeometryBase.h.

{ return fROOTFile; }

TGeoManager * geo::GeometryBase::ROOTGeoManager ( ) const

inherited

Definition at line 1508 of file GeometryBase.cxx.

Referenced by red::EWCosmics::analyze(), geo::GeometryBase::calculateMassesLong(), calib::BetheBlochFit::endJob(), red::OverburdenCalculator::onBeginRun(), supernova::SnovaGen::SnovaGen(), and g4n::FLSHitListAction::SteppingAction().

{ return gGeoManager;}

int geo::GeometryBase::Run ( ) const

inlineinherited

Definition at line 239 of file GeometryBase.h.

{ return fRunNumber; }

void geo::GeometryBase::setDetectorBigBox ( double  detector_bigbox_range )

protectedinherited

Get the Detector size

Dimensions of the DetectorBigBox in cm

Definition at line 1350 of file GeometryBase.cxx.

References geo::GeometryBase::DetectorEnclosureBox(), geo::GeometryBase::fDetectorBigBoxXHi, geo::GeometryBase::fDetectorBigBoxXLo, geo::GeometryBase::fDetectorBigBoxYHi, geo::GeometryBase::fDetectorBigBoxYLo, geo::GeometryBase::fDetectorBigBoxZHi, geo::GeometryBase::fDetectorBigBoxZLo, geo::GeometryBase::fIsDetectorBigBoxUsed, x1, submit_syst::x2, y1, and submit_syst::y2.

Referenced by geo::GeometryBase::LoadGeometryFile().

                                                                   {
    fIsDetectorBigBoxUsed = true;

    double x1, x2;
    double y1, y2;
    double z1, z2;

    /// Get the Detector size
    this->DetectorEnclosureBox(&x1, &x2, &y1, &y2, &z1, &z2);

    /// Dimensions of the DetectorBigBox in cm
    fDetectorBigBoxXHi = x2 + detector_bigbox_range;
    fDetectorBigBoxXLo = x1 - detector_bigbox_range;
    fDetectorBigBoxYHi = y2 + detector_bigbox_range;
    fDetectorBigBoxYLo = y1 - detector_bigbox_range;
    fDetectorBigBoxZHi = z2 + detector_bigbox_range;
    fDetectorBigBoxZLo = z1 - detector_bigbox_range;

  }// end of setDetectorBigBox

void geo::Geometry::setDetectorID ( novadaq::cnv::DetId  det_id )

overrideprotectedvirtual

Method to set DetectorID.

Method to set DetectorID

Reimplemented from geo::GeometryBase.

Definition at line 198 of file Geometry_service.cc.

References DEFINE_ART_SERVICE, geo::GeometryBase::fDetId, ds::DetectorService::SetDetector(), and geo::GeometryBase::setDetectorID().

Referenced by Geometry().

  {
    // Intercept all calls to setDetectorID. First set the parameter in
    // GeometryBase as usual.
    GeometryBase::setDetectorID(detid);

    // Make sure to also inform the detector service
    art::ServiceHandle<ds::DetectorService> ds;
    ds->SetDetector(fDetId);
  }

void geo::GeometryBase::SetDetectorSize ( )

protectedinherited

Calculate the detector size based on the geometry. Requires a little bit of work, so do this once when the geometry is updated

If there no Vertical or horizontal planes, then we go get the size from the Detector Enclosure It used to be that once this condition is true, then it we would just throw an exception. However, now Giuseppe Ferone and Denis Perevalov are studying interactions in simple detector, which don't have planes in it. Thus, they don't want the code breaking.

Getting Detector enclosure volume

Lo and Hi ranges of each of the axes. To be filled later

Fill them

Set detector size

Definition at line 1274 of file GeometryBase.cxx.

References APDHVSetting::dummy, geo::GeometryBase::fDetHalfHeight, geo::GeometryBase::fDetHalfWidth, geo::GeometryBase::fDetLength, geo::GeometryBase::fPlanes, and geo::GeometryBase::Plane().

Referenced by geo::GeometryBase::LoadGeometryFile().

  {
    double dummy;

    // Compute the height and width based on the sizes of the vertical
    // and horizontal planes. Remember, in the local plane frame z goes
    // along the length of the modules. Hence the use of axis=3 to get
    // the detector height (from the vertical planes) and width (from
    // the horizontal planes).
    const TGeoVolume* vvp = gGeoManager->GetVolume("vPlaneV");
    const TGeoVolume* hvp = gGeoManager->GetVolume("vPlaneH");

    // Handle the case of the IPND where the volumes are names slightly
    // different accounting for the different widths of the modules.
    if (vvp==0) vvp = gGeoManager->GetVolume("vPlaneVNDOS");
    if (hvp==0) hvp = gGeoManager->GetVolume("vPlaneHNDOS");
    if (vvp==0 && hvp==0) {
      vvp = gGeoManager->GetVolume("vPlaneVND");
      hvp = gGeoManager->GetVolume("vPlaneHND");
    }

    /// If there no Vertical or horizontal planes, then we go get the size from the Detector Enclosure
    /// It used to be that once this condition is true, then it we would just throw an exception.
    /// However, now Giuseppe Ferone and Denis Perevalov are studying interactions in simple detector,
    /// which don't have planes in it. Thus, they don't want the code breaking.
    if (vvp==0 || hvp==0) {
      mf::LogWarning("SetDetectorSize") 
        << "Warning! BAD_GEO_CONFIG. Unable to find shapes to set detector size\n"
        << "Now getting the detector size from vDetEnclosure\n";
      
      /// Getting Detector enclosure volume
      const TGeoVolume* vdet_enclosure = gGeoManager->GetVolume("vDetEnclosure");
      
      /// Lo and Hi ranges of each of the axes. To be filled later
      double x_lo, x_hi;
      double y_lo, y_hi;
      double z_lo, z_hi;
      
      /// Fill them
      vdet_enclosure->GetShape()->GetAxisRange(1, x_lo, x_hi);
      vdet_enclosure->GetShape()->GetAxisRange(2, y_lo, y_hi);
      vdet_enclosure->GetShape()->GetAxisRange(3, z_lo, z_hi);
      
      /// Set detector size
      fDetHalfWidth  = x_hi;
      fDetHalfHeight = y_hi;
      fDetLength     = z_hi - z_lo;
      
      return;
    }
    
    const TGeoShape* vp = vvp->GetShape();
    const TGeoShape* hp = hvp->GetShape();
    vp->GetAxisRange(3,dummy,fDetHalfHeight);
    hp->GetAxisRange(3,dummy,fDetHalfWidth);
    
    // Remember -- in the local plane coordinate frame, the "depth" is
    // along y, hence axis=2. In the world coordinate this is along z.
    double vplanez1, vplanez2;
    vp->GetAxisRange(2,vplanez1,vplanez2);
    
    // Compute the detector length based on the center positions of
    // cells in the first and last plane. Adjust to account for the
    // plane widths. This ajustment works if the vertical and horizontal
    // planes have the same depth, or if the detector begins and ends
    // with vertical planes. As of the TDR both of these conditions was
    // true for all the NOvA detectors.
    double xyz1[3] = {0,0,0};
    double xyz2[3] = {0,0,0};
    this->Plane(0)->               Cell(0)->GetCenter(xyz1);
    this->Plane(fPlanes.size()-1)->Cell(0)->GetCenter(xyz2);
    fDetLength = (xyz2[2]-xyz1[2])+(vplanez2-vplanez1);
  }// end of SetDetectorSize

void geo::GeometryBase::SetDrawOptions ( )

protectedinherited

Definition at line 691 of file GeometryBase.cxx.

Referenced by geo::GeometryBase::LoadGeometryFile().

{ }

double geo::GeometryBase::SurfaceY ( ) const

inherited

A typical y-position value at the surface (where earth meets air) for this detector site

Returns

typical y position at surface in units of cm

Exceptions

geo::Exception

if detector ID is not set properly

Todo:

: Only far det number is reasonable accurate. Near and IPND need some thought

Definition at line 500 of file GeometryBase.cxx.

References geo::GeometryBase::fDetId, makeTrainCVSamples::int, novadaq::cnv::kFARDET, novadaq::cnv::kNDOS, novadaq::cnv::kNEARDET, and novadaq::cnv::kTESTBEAM.

Referenced by evgen::CosmicsGen::produce(), evgen::CosmicPionGen::produce(), and evgen::CosmicsGen::ProjectCosmicsToSurface().

                                      {

    switch(fDetId){
    case novadaq::cnv::kFARDET:
      return  9.94E2;  // far : surface ~10m  above det. center

    case novadaq::cnv::kNDOS:
      return -1.966E2; // NDOS: surface ~2m   below det. center

    case novadaq::cnv::kNEARDET:
      return 130.0E2;  // near: surface ~130m above det. center

    case novadaq::cnv::kTESTBEAM:
      return -1.966E2; // NDOS: surface ~2m   below det. center

    default:
      throw cet::exception("BadGeoConfig")
        << "Bad geometry configuration: " << (int)fDetId << "\n"
        << __FILE__ << ":" << __LINE__ << "\n";
    }
  }// end of SurfaceY

double geo::GeometryBase::TotalMass ( const char *  volume = "vDetEnclosure" ) const

inherited

Return the total mass of the detector

Definition at line 1036 of file GeometryBase.cxx.

References d, Munits::nm, file_size_ana::node, and string.

Referenced by geo::GeometryTest::analyze().

  {
    double mass = 0.;

    // the detector resides in the vDetEnclosure volume so start there
    const TGeoVolume *enclosure = gGeoManager->FindVolumeFast(volume);

    int old_verbosity = gGeoManager->GetVerboseLevel();
    gGeoManager->SetVerboseLevel(0);
    // loop over the daughters in the volume and let ROOT calculate the
    // mass in kg for you for the daughter nodes that have vBlock in
    // their names this only goes one level down, so no worries of
    // double counting non-superblock blocks
    //
    for(int d = 0; d < enclosure->GetNdaughters(); ++d){
      const TGeoNode* node = enclosure->GetNode(d);
      const char* nm = node->GetName();
      if(strncmp(nm,"vSuperBlock",11) == 0){
        std::string name = nm;
        name.erase(name.find('_'));
        mass += gGeoManager->FindVolumeFast(name.c_str())->Weight();
      }
    }
    gGeoManager->SetVerboseLevel(old_verbosity);
    return mass;
  }// end of TotalMass

void geo::GeometryBase::WorldBox ( double *  xlo, double *  xhi, double *  ylo, double *  yhi, double *  zlo, double *  zhi  ) const

inherited

Return the ranges of x,y and z for the "world volume" that the entire geometry lives in. If any pointers are 0, then those coordinates are ignored.

Parameters

xlo	: On return, lower bound on x positions
xhi	: On return, upper bound on x positions
ylo	: On return, lower bound on y positions
yhi	: On return, upper bound on y positions
zlo	: On return, lower bound on z positions
zhi	: On return, upper bound on z positions

Definition at line 595 of file GeometryBase.cxx.

References ana::assert(), x1, submit_syst::x2, y1, and submit_syst::y2.

Referenced by evgen::CORSIKAGen::GetSample(), g4n::FLSHitListAction::ParticleProjection(), evgen::CosmicsGen::ProjectCosmicsToSurface(), and geo::GeometryTest::testProject().

                                                               {
    const TGeoShape* s = gGeoManager->GetVolume("vWorld")->GetShape();
    assert(s);

    if (xlo || xhi) {
      double x1, x2;
      s->GetAxisRange(1,x1,x2); if (xlo) *xlo = x1; if (xhi) *xhi = x2;
    }
    if (ylo || yhi) {
      double y1, y2;
      s->GetAxisRange(2,y1,y2); if (ylo) *ylo = y1; if (yhi) *yhi = y2;
    }
    if (zlo || zhi) {
      double z1, z2;
      s->GetAxisRange(3,z1,z2); if (zlo) *zlo = z1; if (zhi) *zhi = z2;
    }
  }// end of WorldBox

Member Data Documentation

std::set<unsigned int> geo::GeometryBase::fAllPlanes

protectedinherited

List of all planes.

Definition at line 295 of file GeometryBase.h.

Referenced by geo::GeometryBase::BuildMaps(), and geo::GeometryBase::GetPlanesByView().

double geo::GeometryBase::fBigBoxRange

protectedinherited

Range of big box.

Definition at line 311 of file GeometryBase.h.

Referenced by geo::GeometryBase::LoadGeometryFile(), and geo::GeometryBase::setDetectorID().

fhicl::ParameterSet geo::GeometryBase::fConfigFD

protectedinherited

Definition at line 321 of file GeometryBase.h.

fhicl::ParameterSet geo::GeometryBase::fConfigND

protectedinherited

Definition at line 321 of file GeometryBase.h.

fhicl::ParameterSet geo::GeometryBase::fConfigNDOS

protectedinherited

Definition at line 321 of file GeometryBase.h.

fhicl::ParameterSet geo::GeometryBase::fConfigTB

protectedinherited

Definition at line 321 of file GeometryBase.h.

CoordinateTransformation geo::GeometryBase::fCoordinateTransformation

protectedinherited

Coordinate Transformation class.

Definition at line 289 of file GeometryBase.h.

Referenced by geo::GeometryBase::setDetectorID().

double geo::GeometryBase::fDetectorBigBoxXHi

protectedinherited

Dimensions of the DetectorBigBox in cm.

Definition at line 314 of file GeometryBase.h.

Referenced by geo::GeometryBase::DetectorBigBox(), geo::GeometryBase::IntersectsBigBox(), geo::GeometryBase::LoadGeometryFile(), and geo::GeometryBase::setDetectorBigBox().

double geo::GeometryBase::fDetectorBigBoxXLo

protectedinherited

Definition at line 315 of file GeometryBase.h.

Referenced by geo::GeometryBase::DetectorBigBox(), geo::GeometryBase::IntersectsBigBox(), geo::GeometryBase::isInsideDetectorBigBox(), geo::GeometryBase::LoadGeometryFile(), and geo::GeometryBase::setDetectorBigBox().

double geo::GeometryBase::fDetectorBigBoxYHi

protectedinherited

Definition at line 316 of file GeometryBase.h.

Referenced by geo::GeometryBase::DetectorBigBox(), geo::GeometryBase::IntersectsBigBox(), geo::GeometryBase::LoadGeometryFile(), and geo::GeometryBase::setDetectorBigBox().

double geo::GeometryBase::fDetectorBigBoxYLo

protectedinherited

Definition at line 317 of file GeometryBase.h.

Referenced by geo::GeometryBase::DetectorBigBox(), geo::GeometryBase::IntersectsBigBox(), geo::GeometryBase::isInsideDetectorBigBox(), geo::GeometryBase::LoadGeometryFile(), and geo::GeometryBase::setDetectorBigBox().

double geo::GeometryBase::fDetectorBigBoxZHi

protectedinherited

Definition at line 318 of file GeometryBase.h.

Referenced by geo::GeometryBase::DetectorBigBox(), geo::GeometryBase::IntersectsBigBox(), geo::GeometryBase::LoadGeometryFile(), and geo::GeometryBase::setDetectorBigBox().

double geo::GeometryBase::fDetectorBigBoxZLo

protectedinherited

Definition at line 319 of file GeometryBase.h.

Referenced by geo::GeometryBase::DetectorBigBox(), geo::GeometryBase::IntersectsBigBox(), geo::GeometryBase::isInsideDetectorBigBox(), geo::GeometryBase::LoadGeometryFile(), and geo::GeometryBase::setDetectorBigBox().

double geo::GeometryBase::fDetHalfHeight

protectedinherited

Detector 1/2 height (cm)

Definition at line 300 of file GeometryBase.h.

Referenced by geo::GeometryBase::calculateMassesLong(), geo::GeometryBase::DEdge(), geo::GeometryBase::DetHalfHeight(), geo::GeometryBase::IntersectsDetector(), geo::GeometryBase::LoadGeometryFile(), geo::GeometryBase::setDetectorID(), and geo::GeometryBase::SetDetectorSize().

double geo::GeometryBase::fDetHalfWidth

protectedinherited

Detector 1/2 width (cm)

Definition at line 301 of file GeometryBase.h.

Referenced by geo::GeometryBase::calculateMassesLong(), geo::GeometryBase::DEdge(), geo::GeometryBase::DetHalfWidth(), geo::GeometryBase::IntersectsDetector(), geo::GeometryBase::LoadGeometryFile(), geo::GeometryBase::setDetectorID(), and geo::GeometryBase::SetDetectorSize().

novadaq::cnv::DetId geo::GeometryBase::fDetId

protectedinherited

id for the detector being used

Definition at line 285 of file GeometryBase.h.

Referenced by geo::GeometryBase::BoosterBeamDirection(), geo::GeometryBase::DetHalfHeight(), geo::GeometryBase::DetHalfWidth(), geo::GeometryBase::DetLength(), Geometry(), geo::GeometryBase::NPlanes(), geo::GeometryBase::NuMIBeamDirection(), setDetectorID(), geo::GeometryBase::setDetectorID(), and geo::GeometryBase::SurfaceY().

double geo::GeometryBase::fDetLength

protectedinherited

Detector length (cm)

Definition at line 299 of file GeometryBase.h.

Referenced by geo::GeometryBase::calculateMassesLong(), geo::GeometryBase::DEdge(), geo::GeometryBase::DetLength(), geo::GeometryBase::IntersectsDetector(), geo::GeometryBase::LoadGeometryFile(), geo::GeometryBase::setDetectorID(), and geo::GeometryBase::SetDetectorSize().

double geo::GeometryBase::fFiducialVolumeXHi

protectedinherited

Definition at line 323 of file GeometryBase.h.

Referenced by geo::GeometryBase::FiducialBox(), and geo::GeometryBase::setDetectorID().

double geo::GeometryBase::fFiducialVolumeXLo

protectedinherited

Definition at line 324 of file GeometryBase.h.

Referenced by geo::GeometryBase::FiducialBox(), geo::GeometryBase::isInsideFiducialVolume(), and geo::GeometryBase::setDetectorID().

double geo::GeometryBase::fFiducialVolumeYHi

protectedinherited

Definition at line 325 of file GeometryBase.h.

Referenced by geo::GeometryBase::FiducialBox(), and geo::GeometryBase::setDetectorID().

double geo::GeometryBase::fFiducialVolumeYLo

protectedinherited

Definition at line 326 of file GeometryBase.h.

Referenced by geo::GeometryBase::FiducialBox(), geo::GeometryBase::isInsideFiducialVolume(), and geo::GeometryBase::setDetectorID().

double geo::GeometryBase::fFiducialVolumeZHi

protectedinherited

Definition at line 327 of file GeometryBase.h.

Referenced by geo::GeometryBase::FiducialBox(), and geo::GeometryBase::setDetectorID().

double geo::GeometryBase::fFiducialVolumeZLo

protectedinherited

Definition at line 328 of file GeometryBase.h.

Referenced by geo::GeometryBase::FiducialBox(), geo::GeometryBase::isInsideFiducialVolume(), and geo::GeometryBase::setDetectorID().

std::string geo::GeometryBase::fGDMLFile

protectedinherited

gdml file holding the geometry

Definition at line 281 of file GeometryBase.h.

Referenced by geo::GeometryBase::ExtractGDML(), geo::GeometryBase::FileBaseName(), geo::GeometryBase::GeometryBase(), geo::GeometryBase::MakeTmpFile(), preBeginRun(), ~Geometry(), and geo::GeometryBase::~GeometryBase().

std::string geo::GeometryBase::fGDMLFromFCL

protectedinherited

keep track of original fcl parameter

Definition at line 283 of file GeometryBase.h.

Referenced by geo::GeometryBase::GeometryBase(), and preBeginRun().

UIDMap geo::GeometryBase::fIdMap

protectedinherited

Unique ID -> Plane,Cell.

Definition at line 293 of file GeometryBase.h.

Referenced by geo::GeometryBase::BuildMaps(), geo::GeometryBase::getPlaneCellMap(), and geo::GeometryBase::LoadGeometryFile().

bool geo::GeometryBase::fIsDetectorBigBoxUsed

protectedinherited

Do we need to use the BigBox cut?

Definition at line 310 of file GeometryBase.h.

Referenced by geo::GeometryBase::LoadGeometryFile(), geo::GeometryBase::setDetectorBigBox(), and geo::GeometryBase::setDetectorID().

std::string geo::Geometry::fLoadedGDMLInfo

protected

Definition at line 44 of file Geometry.h.

Referenced by preBeginRun().

Params geo::GeometryBase::fParams

protectedinherited

Definition at line 279 of file GeometryBase.h.

Referenced by geo::GeometryBase::GeometryBase(), geo::GeometryBase::MakeTmpFile(), preBeginRun(), and geo::GeometryBase::setDetectorID().

std::vector<PlaneGeo*> geo::GeometryBase::fPlanes

protectedinherited

The detector planes.

Definition at line 291 of file GeometryBase.h.

Referenced by geo::GeometryBase::BuildMaps(), geo::GeometryBase::CountCellsOnLineFast(), geo::GeometryBase::FirstPlaneInMuonCatcher(), geo::GeometryBase::LoadGeometryFile(), geo::GeometryBase::MakePlane(), geo::GeometryBase::NextPlaneInView(), geo::GeometryBase::NextPlaneOtherView(), geo::GeometryBase::NPlanes(), geo::GeometryBase::Plane(), geo::GeometryBase::SetDetectorSize(), and geo::GeometryBase::~GeometryBase().

std::vector<PlaneGeo*> geo::GeometryBase::fPlanesInMuonCatcher

protectedinherited

Same pointers as fPlanes.

Definition at line 292 of file GeometryBase.h.

Referenced by geo::GeometryBase::FirstPlaneInMuonCatcher(), geo::GeometryBase::LoadGeometryFile(), and geo::GeometryBase::MakePlane().

std::string geo::GeometryBase::fROOTFile

protectedinherited

root file holding the geometry

Definition at line 282 of file GeometryBase.h.

Referenced by geo::GeometryBase::GeometryBase(), and preBeginRun().

unsigned int geo::GeometryBase::fRunNumber

protectedinherited

Run number of configuration.

Definition at line 332 of file GeometryBase.h.

Referenced by preBeginRun().

std::string geo::GeometryBase::fStoreTempGeo

protectedinherited

folder to store temp geo file

Definition at line 330 of file GeometryBase.h.

std::set<unsigned int> geo::GeometryBase::fXplanes

protectedinherited

List of X measuring planes.

Definition at line 296 of file GeometryBase.h.

Referenced by geo::GeometryBase::BuildMaps(), and geo::GeometryBase::GetPlanesByView().

std::set<unsigned int> geo::GeometryBase::fYplanes

protectedinherited

List of Y measuring planes.

Definition at line 297 of file GeometryBase.h.

Referenced by geo::GeometryBase::BuildMaps(), and geo::GeometryBase::GetPlanesByView().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-14/Geometry/Geometry.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-04-14/Geometry/Geometry_service.cc