#include <MagBLayer.h>

Public Member Functions
const MagVolume *	findVolume (const GlobalPoint &gp, double tolerance) const
	Find the volume containing a point, with a given tolerance. More...

	MagBLayer (std::vector< MagBSector * > &sectors, double rMin)
	Constructor. More...

	MagBLayer (MagVolume *aVolume, double rMin)
	Constructor for a trivial layer consisting of one single volume. More...

double	minR () const
	Lowest radius of the layer. More...

virtual	~MagBLayer ()
	Destructor. More...

Private Attributes
PeriodicBinFinderInPhi< float > *	theBinFinder

double	theRMin

std::vector< MagBSector * >	theSectors

MagVolume *	theSingleVolume

Detailed Description

A barrel layer (MagBLayer) groups volumes at the same distance to the origin. It consists of either 1 single volume (a cylinder) or 12 sectors in phi (MagBSector). Each sector consists of one or more rods (MagBRods) of equal width in phi. Rods consist of one or more slabs (MagBSlab); each one consisting of one or, in few cases, several volumes with the same lenght in Z.

Author: N. Amapane - INFN Torino

Definition at line 24 of file MagBLayer.h.

Constructor & Destructor Documentation

MagBLayer::MagBLayer	(	std::vector< MagBSector * > &	sectors,
		double	rMin
	)

Constructor.

Definition at line 21 of file MagBLayer.cc.

References Geom::pi(), theBinFinder, and theSectors.

     : theSectors(sectors), theSingleVolume(nullptr), theRMin(rMin), theBinFinder(nullptr) {
   //   LogTrace("MagGeometry") << "Building MagBLayer with " << theSectors.size()
   //          << " sectors, minR " << theRMin << endl;
   //FIXME: PeriodicBinFinderInPhi gets *center* of first bin
   theBinFinder = new PeriodicBinFinderInPhi<float>(theSectors.front()->minPhi() + Geom::pi() / 12., 12);
 }

MagBLayer::MagBLayer	(	MagVolume *	aVolume,
		double	rMin
	)

Constructor for a trivial layer consisting of one single volume.

Definition at line 30 of file MagBLayer.cc.

                                                     : theSingleVolume(nullptr), theRMin(rMin), theBinFinder(nullptr) {
   //   LogTrace("MagGeometry") << "Building MagBLayer with " << 0
   //          << " sectors, minR " << theRMin << endl;
 }

MagBLayer::~MagBLayer ( )

virtual

Destructor.

Definition at line 35 of file MagBLayer.cc.

References theBinFinder, theSectors, and theSingleVolume.

                       {
   delete theBinFinder;
 
   delete theSingleVolume;
 
   for (vector<MagBSector*>::const_iterator isec = theSectors.begin(); isec != theSectors.end(); ++isec) {
     delete (*isec);
   }
 }

Member Function Documentation

const MagVolume * MagBLayer::findVolume	(	const GlobalPoint &	gp,
		double	tolerance
	)		const

Find the volume containing a point, with a given tolerance.

Definition at line 45 of file MagBLayer.cc.

References newFWLiteAna::bin, PeriodicBinFinderInPhi< T >::binIndex(), MagVolume::inside(), LogTrace, phi, PV3DBase< T, PVType, FrameType >::phi(), mps_fire::result, theBinFinder, theSectors, and theSingleVolume.

                                                                                     {
   const MagVolume* result = nullptr;
 
   //In case the layer is composed of a single volume...
   if (theSingleVolume) {
     //    LogTrace("MagGeometry") << "   Trying the unique volume " << endl;
     if (theSingleVolume->inside(gp, tolerance)) {
       result = theSingleVolume;
       //       LogTrace("MagGeometry") << "***In unique bsector"
       //          << (result==0? " failed " : " OK ") <<endl;
     }
     return result;
   }
 
   // Normal cases - query the sectors.
 
   Geom::Phi<float> phi = gp.phi();
 
   // FIXME assume sectors are sorted in phi!
   int bin = theBinFinder->binIndex(phi);
   LogTrace("MagGeometry") << "   Trying sector at phi " << theSectors[bin]->minPhi() << " " << phi << endl;
   result = theSectors[bin]->findVolume(gp, tolerance);
   LogTrace("MagGeometry") << "***In guessed bsector" << (result == nullptr ? " failed " : " OK ") << endl;
 
   if (result == nullptr) {  // If fails, can be in previous bin.
     LogTrace("MagGeometry") << "   Trying sector at phi " << theSectors[theBinFinder->binIndex(bin - 1)]->minPhi()
                             << " " << phi << endl;
 
     result = theSectors[theBinFinder->binIndex(bin - 1)]->findVolume(gp, tolerance);
     LogTrace("MagGeometry") << "***In previous bsector" << (result == nullptr ? " failed " : " OK ") << endl;
   }
   return result;
 }