#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 23 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(0),
   theRMin(rMin),
   theBinFinder(0)
 {
   // TOFIX
 //   if (verbose.debugOut) cout << "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 36 of file MagBLayer.cc.

                                                     :
   theSingleVolume(0),
   theRMin(rMin),
   theBinFinder(0) 
 {
   // TOFIX
 //   if (verbose.debugOut) cout << "Building MagBLayer with " << 0
 //        << " sectors, minR " << theRMin << endl;    
 }

MagBLayer::~MagBLayer ( )

virtual

Destructor.

Definition at line 48 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 60 of file MagBLayer.cc.

References newFWLiteAna::bin, PeriodicBinFinderInPhi< T >::binIndex(), gather_cfg::cout, verbose::debugOut, MagVolume::inside(), phi, PV3DBase< T, PVType, FrameType >::phi(), query::result, theBinFinder, theSectors, and theSingleVolume.

                                                                                      {
   const MagVolume * result = 0;
 
   //In case the layer is composed of a single volume...
   if (theSingleVolume) {
     // TOFIX
 //     if (verbose.debugOut) cout << "   Trying the unique volume " << endl;
     if (theSingleVolume->inside(gp, tolerance)) {
       result = theSingleVolume;
     // TOFIX
 //       if (verbose.debugOut) cout << "***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);
     // TOFIX
   if (verbose::debugOut) cout << "   Trying sector at phi " << theSectors[bin]->minPhi()
                   << " " << phi << endl ;
   result = theSectors[bin]->findVolume(gp, tolerance);
     // TOFIX
   if (verbose::debugOut) cout << "***In guessed bsector"
                  << (result==0? " failed " : " OK ") <<endl;
 
   if (result==0) { // If fails, can be in previous bin.
     // TOFIX
     if (verbose::debugOut) cout << "   Trying sector at phi "
                    << theSectors[theBinFinder->binIndex(bin-1)]->minPhi()
                    << " " << phi << endl ;
     
     result = theSectors[theBinFinder->binIndex(bin-1)]->findVolume(gp, tolerance);
     // TOFIX
     if (verbose::debugOut) cout << "***In previous bsector"
                    << (result==0? " failed " : " OK ") <<endl;
 
   }
   return result;
 
 }