#include <MagGeometry.h>

Public Types
typedef Surface::GlobalPoint	GlobalPoint

typedef Surface::GlobalVector	GlobalVector

Public Member Functions
const std::vector < MagVolume6Faces const * > &	barrelVolumes () const

const std::vector < MagVolume6Faces const * > &	endcapVolumes () const

GlobalVector	fieldInTesla (const GlobalPoint &gp) const
	Return field vector at the specified global point. More...

MagVolume const *	findVolume (const GlobalPoint &gp, double tolerance=0.) const
	Find a volume. More...

	MagGeometry (int geomVersion, const std::vector< MagBLayer * > &, const std::vector< MagESector * > &, const std::vector< MagVolume6Faces * > &, const std::vector< MagVolume6Faces * > &)
	Constructor. More...

	MagGeometry (int geomVersion, const std::vector< MagBLayer const * > &, const std::vector< MagESector const * > &, const std::vector< MagVolume6Faces const * > &, const std::vector< MagVolume6Faces const * > &)

	~MagGeometry ()
	Destructor. More...

Private Member Functions
MagVolume const *	findVolume1 (const GlobalPoint &gp, double tolerance=0.) const

bool	inBarrel (const GlobalPoint &gp) const

Private Attributes
bool	cacheLastVolume

int	geometryVersion

const int	me_

MagBinFinders::GeneralBinFinderInR < double > const *	theBarrelBinFinder

double	theBarrelRsq1

double	theBarrelRsq2

double	theBarrelZ0

double	theBarrelZ1

double	theBarrelZ2

std::vector< MagBLayer const * >	theBLayers

std::vector< MagVolume6Faces const * >	theBVolumes

PeriodicBinFinderInPhi< float > const *	theEndcapBinFinder

std::vector< MagESector const * >	theESectors

std::vector< MagVolume6Faces const * >	theEVolumes

Friends
class	MagGeometryExerciser

Detailed Description

Entry point to the geometry of magnetic volumes.

Author: N. Amapane - INFN Torino

Definition at line 23 of file MagGeometry.h.

Member Typedef Documentation

typedef Surface::GlobalPoint MagGeometry::GlobalPoint

Definition at line 26 of file MagGeometry.h.

typedef Surface::GlobalVector MagGeometry::GlobalVector

Definition at line 25 of file MagGeometry.h.

Constructor & Destructor Documentation

MagGeometry::MagGeometry	(	int	geomVersion,
		const std::vector< MagBLayer * > &	tbl,
		const std::vector< MagESector * > &	tes,
		const std::vector< MagVolume6Faces * > &	tbv,
		const std::vector< MagVolume6Faces * > &	tev
	)

Constructor.

Definition at line 32 of file MagGeometry.cc.

     : MagGeometry(geomVersion,
                   reinterpret_cast<std::vector<MagBLayer const*> const&>(tbl),
                   reinterpret_cast<std::vector<MagESector const*> const&>(tes),
                   reinterpret_cast<std::vector<MagVolume6Faces const*> const&>(tbv),
                   reinterpret_cast<std::vector<MagVolume6Faces const*> const&>(tev)) {}

MagGeometry::MagGeometry	(	int	geomVersion,
		const std::vector< MagBLayer const * > &	tbl,
		const std::vector< MagESector const * > &	tes,
		const std::vector< MagVolume6Faces const * > &	tbv,
		const std::vector< MagVolume6Faces const * > &	tev
	)

Definition at line 43 of file MagGeometry.cc.

References LogTrace, me_, Geom::pi(), theBarrelBinFinder, theBarrelRsq1, theBarrelRsq2, theBarrelZ0, theBarrelZ1, theBarrelZ2, theBLayers, theEndcapBinFinder, and theESectors.

     : me_(++instanceCounter),
       theBLayers(tbl),
       theESectors(tes),
       theBVolumes(tbv),
       theEVolumes(tev),
       cacheLastVolume(true),
       geometryVersion(geomVersion) {
   vector<double> rBorders;
 
   for (vector<MagBLayer const*>::const_iterator ilay = theBLayers.begin(); ilay != theBLayers.end(); ++ilay) {
     LogTrace("MagGeoBuilder") << "  Barrel layer at " << (*ilay)->minR() << endl;
     //FIXME assume layers are already sorted in minR
     rBorders.push_back((*ilay)->minR() * (*ilay)->minR());
   }
 
   theBarrelBinFinder = new MagBinFinders::GeneralBinFinderInR<double>(rBorders);
 
 #ifdef EDM_ML_DEBUG
   for (vector<MagESector const*>::const_iterator isec = theESectors.begin(); isec != theESectors.end(); ++isec) {
     LogTrace("MagGeoBuilder") << "  Endcap sector at " << (*isec)->minPhi() << endl;
   }
 #endif
 
   //FIXME assume sectors are already sorted in phi
   //FIXME: PeriodicBinFinderInPhi gets *center* of first bin
   int nEBins = theESectors.size();
   if (nEBins > 0)
     theEndcapBinFinder = new PeriodicBinFinderInPhi<float>(theESectors.front()->minPhi() + Geom::pi() / nEBins, nEBins);
 
   // Compute barrel dimensions based on geometry version
   // FIXME: it would be nice to derive these from the actual geometry in the builder, possibly adding some specification to the geometry.
   switch (geomVersion >= 120812 ? 0 : (geomVersion >= 90812 ? 1 : 2)) {
     case 0:  // since 120812
       theBarrelRsq1 = 172.400 * 172.400;
       theBarrelRsq2 = 308.735 * 308.735;
       theBarrelZ0 = 350.000;
       theBarrelZ1 = 633.290;
       theBarrelZ2 = 662.010;
       break;
     case 1:  // version 90812 (no longer in use)
       theBarrelRsq1 = 172.400 * 172.400;
       theBarrelRsq2 = 308.755 * 308.755;
       theBarrelZ0 = 350.000;
       theBarrelZ1 = 633.890;
       theBarrelZ2 = 662.010;
       break;
     case 2:  // versions 71212, 90322
       theBarrelRsq1 = 172.400 * 172.400;
       theBarrelRsq2 = 308.755 * 308.755;
       theBarrelZ0 = 350.000;
       theBarrelZ1 = 633.290;
       theBarrelZ2 = 661.010;
       break;
   }
 
   LogTrace("MagGeometry_cache") << "*** In MagGeometry ctor: me_=" << me_ << " instanceCounter=" << instanceCounter
                                 << endl;
 }

MagGeometry::~MagGeometry ( )

Destructor.

Definition at line 107 of file MagGeometry.cc.

References theBarrelBinFinder, theBLayers, theEndcapBinFinder, and theESectors.

                           {
   if (theBarrelBinFinder != nullptr)
     delete theBarrelBinFinder;
   if (theEndcapBinFinder != nullptr)
     delete theEndcapBinFinder;
 
   for (vector<MagBLayer const*>::const_iterator ilay = theBLayers.begin(); ilay != theBLayers.end(); ++ilay) {
     delete (*ilay);
   }
 
   for (vector<MagESector const*>::const_iterator ilay = theESectors.begin(); ilay != theESectors.end(); ++ilay) {
     delete (*ilay);
   }
 }

Member Function Documentation

const std::vector<MagVolume6Faces const*>& MagGeometry::barrelVolumes ( ) const

inline

Definition at line 50 of file MagGeometry.h.

References theBVolumes.

50 { return theBVolumes; }

MagGeometry::theBVolumes

std::vector< MagVolume6Faces const * > theBVolumes

Definition: MagGeometry.h:67

const std::vector<MagVolume6Faces const*>& MagGeometry::endcapVolumes ( ) const

inline

Definition at line 51 of file MagGeometry.h.

References theEVolumes.

51 { return theEVolumes; }

MagGeometry::theEVolumes

std::vector< MagVolume6Faces const * > theEVolumes

Definition: MagGeometry.h:68

GlobalVector MagGeometry::fieldInTesla ( const GlobalPoint & gp ) const

Return field vector at the specified global point.

Definition at line 123 of file MagGeometry.cc.

References MagVolume::fieldInTesla(), findVolume(), edm::isNotFinite(), PV3DBase< T, PVType, FrameType >::mag(), and findQualityFiles::v.

Referenced by VolumeBasedMagneticField::inTesla(), and VolumeBasedMagneticField::inTeslaUnchecked().

                                                                   {
   MagVolume const* v = nullptr;
 
   v = findVolume(gp);
   if (v != nullptr) {
     return v->fieldInTesla(gp);
   }
 
   // Fall-back case: no volume found
 
   if (edm::isNotFinite(gp.mag())) {
     LogWarning("MagneticField") << "Input value invalid (not a number): " << gp << endl;
 
   } else {
     LogWarning("MagneticField") << "MagGeometry::fieldInTesla: failed to find volume for " << gp << endl;
   }
   return GlobalVector();
 }

MagVolume const * MagGeometry::findVolume	(	const GlobalPoint &	gp,
		double	tolerance = `0.`
	)		const

Find a volume.

Definition at line 184 of file MagGeometry.cc.

References newFWLiteAna::bin, PeriodicBinFinderInPhi< T >::binIndex(), MagBinFinders::GeneralBinFinderInR< T >::binIndex(), cacheLastVolume, inBarrel(), LogTrace, SiStripPI::max, me_, PV3DBase< T, PVType, FrameType >::perp2(), phi, PV3DBase< T, PVType, FrameType >::phi(), mps_fire::result, mathSSE::sqrt(), theBarrelBinFinder, theBLayers, theEndcapBinFinder, and theESectors.

Referenced by fieldInTesla(), and VolumeBasedMagneticField::findVolume().

                                                                                       {
   // Clear volume cache if this is a new instance
   if (me_ != localInstance) {
     LogTrace("MagGeometry_cache") << "*** In MagGeometry::findVolume resetting cache: me=" << me_
                                   << " localInstance=" << localInstance << endl;
     localInstance = me_;
     lastVolume = nullptr;
   }
 
   if (lastVolume != nullptr && lastVolume->inside(gp)) {
     return lastVolume;
   }
 
   MagVolume const* result = nullptr;
   if (inBarrel(gp)) {  // Barrel
     double aRsq = gp.perp2();
     int bin = theBarrelBinFinder->binIndex(aRsq);
 
     // Search up to 3 layers inwards. This may happen for very thin layers.
     for (int bin1 = bin; bin1 >= max(0, bin - 3); --bin1) {
       LogTrace("MagGeometry") << "Trying layer at R " << theBLayers[bin1]->minR() << " " << sqrt(aRsq) << endl;
       result = theBLayers[bin1]->findVolume(gp, tolerance);
       LogTrace("MagGeometry") << "***In blayer " << bin1 - bin << " " << (result == nullptr ? " failed " : " OK ")
                               << endl;
       if (result != nullptr)
         break;
     }
 
   } else {  // Endcaps
     Geom::Phi<float> phi = gp.phi();
     if (theEndcapBinFinder != nullptr && !theESectors.empty()) {
       int bin = theEndcapBinFinder->binIndex(phi);
       LogTrace("MagGeometry") << "Trying endcap sector at phi " << theESectors[bin]->minPhi() << " " << phi << endl;
       result = theESectors[bin]->findVolume(gp, tolerance);
       LogTrace("MagGeometry") << "***In guessed esector " << (result == nullptr ? " failed " : " OK ") << endl;
     } else
       edm::LogError("MagGeometry") << "Endcap empty";
   }
 
   if (result == nullptr && tolerance < 0.0001) {
     // If search fails, retry with a 300 micron tolerance.
     // This is a hack for thin gaps on air-iron boundaries,
     // which will not be present anymore once surfaces are matched.
     LogTrace("MagGeometry") << "Increasing the tolerance to 0.03" << endl;
     result = findVolume(gp, 0.03);
   }
 
   if (cacheLastVolume)
     lastVolume = result;
 
   return result;
 }

MagVolume const * MagGeometry::findVolume1	(	const GlobalPoint &	gp,
		double	tolerance = `0.`
	)		const

private

Definition at line 143 of file MagGeometry.cc.

References newFWLiteAna::found, inBarrel(), theBVolumes, theEVolumes, and findQualityFiles::v.

                                                                                        {
   MagVolume6Faces const* found = nullptr;
 
   int errCnt = 0;
   if (inBarrel(gp)) {  // Barrel
     for (vector<MagVolume6Faces const*>::const_iterator v = theBVolumes.begin(); v != theBVolumes.end(); ++v) {
       if ((*v) == nullptr) {  //FIXME: remove this check
         LogError("MagGeometry") << endl << "***ERROR: MagGeometry::findVolume: MagVolume for barrel not set" << endl;
         ++errCnt;
         if (errCnt < 3)
           continue;
         else
           break;
       }
       if ((*v)->inside(gp, tolerance)) {
         found = (*v);
         break;
       }
     }
 
   } else {  // Endcaps
     for (vector<MagVolume6Faces const*>::const_iterator v = theEVolumes.begin(); v != theEVolumes.end(); ++v) {
       if ((*v) == nullptr) {  //FIXME: remove this check
         LogError("MagGeometry") << endl << "***ERROR: MagGeometry::findVolume: MagVolume for endcap not set" << endl;
         ++errCnt;
         if (errCnt < 3)
           continue;
         else
           break;
       }
       if ((*v)->inside(gp, tolerance)) {
         found = (*v);
         break;
       }
     }
   }
 
   return found;
 }

bool MagGeometry::inBarrel ( const GlobalPoint & gp ) const

private

Definition at line 237 of file MagGeometry.cc.

References PV3DBase< T, PVType, FrameType >::perp2(), theBarrelRsq1, theBarrelRsq2, theBarrelZ0, and PV3DBase< T, PVType, FrameType >::z().

Referenced by findVolume(), and findVolume1().

                                                       {
   double aZ = fabs(gp.z());
   double aRsq = gp.perp2();
 
   return ((aZ < theBarrelZ0) || (aZ < theBarrelZ1 && aRsq > theBarrelRsq1) ||
           (aZ < theBarrelZ2 && aRsq > theBarrelRsq2));
 }