#include <NavVolume6Faces.h>

Inheritance diagram for NavVolume6Faces:

Public Member Functions
VolumeCrossReturnType	crossToNextVolume (const TrajectoryStateOnSurface &currentState, const Propagator &prop) const override
	Cross this volume and point at the next. More...

const std::vector< VolumeSide > &	faces () const override
	Access to volume faces. More...

bool	inside (const GlobalPoint &gp, double tolerance) const override

bool	isIron () const override
	Access to Iron/Air information: More...

	NavVolume6Faces (const PositionType &pos, const RotationType &rot, DDSolidShape shape, const std::vector< NavVolumeSide > &faces, const MagneticFieldProvider< float > *mfp)

	NavVolume6Faces (const MagVolume &magvol, const bool isIron=false)
	A NavVolume6Faces that corresponds exactly to a MagVolume. More...

Container	nextSurface (const NavVolume::LocalPoint &pos, const NavVolume::LocalVector &mom, double charge, PropagationDirection propDir=alongMomentum) const override
	Give a sorted list of possible surfaces to propagate to. More...

Container	nextSurface (const NavVolume::LocalPoint &pos, const NavVolume::LocalVector &mom, double charge, PropagationDirection propDir, ConstReferenceCountingPointer< Surface > NotThisSurfaceP) const override
	Same, giving lowest priority to the surface we are on now (=NotThisSurface) More...

Public Member Functions inherited from NavVolume
	NavVolume (const PositionType &pos, const RotationType &rot, DDSolidShape shape, const MagneticFieldProvider< float > *mfp)

	~NavVolume () override

Public Member Functions inherited from MagVolume
LocalVector	fieldInTesla (const LocalPoint &lp) const

GlobalVector	fieldInTesla (const GlobalPoint &lp) const

virtual bool	inside (const LocalPoint &lp, double tolerance=0.) const

::GlobalVector	inTesla (const ::GlobalPoint &gp) const override

bool	isIron () const
	Temporary hack to pass information on material. Will eventually be replaced! More...

	MagVolume (const PositionType &pos, const RotationType &rot, DDSolidShape shape, const MagneticFieldProvider< float > *mfp, double sf=1.)

void	ownsFieldProvider (bool o)

const MagneticFieldProvider< float > *	provider () const

void	setIsIron (bool iron)

DDSolidShape	shapeType () const

	~MagVolume () override

Public Member Functions inherited from GloballyPositioned< float >
float	eta () const

	GloballyPositioned ()

	GloballyPositioned (const PositionType &pos, const RotationType &rot)

void	move (const GlobalVector &displacement)

float	phi () const

const PositionType &	position () const

void	rotate (const RotationType &rotation)

const RotationType &	rotation () const

GlobalPoint	toGlobal (const LocalPoint &lp) const

Point3DBase< U, GlobalTag >	toGlobal (const Point3DBase< U, LocalTag > &lp) const

GlobalVector	toGlobal (const LocalVector &lv) const

Vector3DBase< U, GlobalTag >	toGlobal (const Vector3DBase< U, LocalTag > &lv) const

LocalPoint	toLocal (const GlobalPoint &gp) const

Point3DBase< U, LocalTag >	toLocal (const Point3DBase< U, GlobalTag > &gp) const

LocalVector	toLocal (const GlobalVector &gv) const

Vector3DBase< U, LocalTag >	toLocal (const Vector3DBase< U, GlobalTag > &gv) const

virtual	~GloballyPositioned ()

Public Member Functions inherited from MagneticField
virtual MagneticField *	clone () const

GlobalVector	inInverseGeV (const GlobalPoint &gp) const
	Field value ad specified global point, in 1/Gev. More...

GlobalVector	inKGauss (const GlobalPoint &gp) const
	Field value ad specified global point, in KGauss. More...

virtual GlobalVector	inTesla (const GlobalPoint &gp) const =0
	Field value ad specified global point, in Tesla. More...

virtual GlobalVector	inTeslaUnchecked (const GlobalPoint &gp) const

virtual bool	isDefined (const GlobalPoint &) const
	True if the point is within the region where the concrete field. More...

	MagneticField ()

	MagneticField (const MagneticField &orig)

int	nominalValue () const
	The nominal field value for this map in kGauss. More...

virtual	~MagneticField ()

Private Member Functions
void	computeBounds (const std::vector< NavVolumeSide > &faces)

Bounds *	computeBounds (int index, const std::vector< const Plane * > &bpc)

Bounds *	computeBounds (int index, const std::vector< NavVolumeSide > &faces)

Private Attributes
bool	isThisIron

std::vector< VolumeSide >	theFaces

Container	theNavSurfaces

Additional Inherited Members
Public Types inherited from NavVolume
typedef std::vector< SurfaceAndBounds >	Container

Public Types inherited from MagVolume
typedef GloballyPositioned< float >	Base

typedef GloballyPositioned< float >::GlobalPoint	GlobalPoint

typedef GloballyPositioned< float >::GlobalVector	GlobalVector

typedef GloballyPositioned< float >::LocalPoint	LocalPoint

typedef GloballyPositioned< float >::LocalVector	LocalVector

Public Types inherited from GloballyPositioned< float >
typedef Point3DBase< float, GlobalTag >	GlobalPoint

typedef Vector3DBase< float, GlobalTag >	GlobalVector

typedef Point3DBase< float, LocalTag >	LocalPoint

typedef Vector3DBase< float, LocalTag >	LocalVector

typedef Point3DBase< float, GlobalTag >	PositionType

typedef TkRotation< float >	RotationType

typedef float	Scalar

Static Public Member Functions inherited from GloballyPositioned< float >
static float	iniEta ()

static float	iniPhi ()

Detailed Description

Definition at line 13 of file NavVolume6Faces.h.

Constructor & Destructor Documentation

NavVolume6Faces::NavVolume6Faces	(	const PositionType &	pos,
		const RotationType &	rot,
		DDSolidShape	shape,
		const std::vector< NavVolumeSide > &	faces,
		const MagneticFieldProvider< float > *	mfp
	)

Definition at line 25 of file NavVolume6Faces.cc.

References computeBounds(), mps_fire::i, and theFaces.

                                                             :
   NavVolume(pos,rot,shape,mfp) 
 {
   for (std::vector<NavVolumeSide>::const_iterator i=faces.begin();
        i != faces.end(); i++) {
     theFaces.push_back( VolumeSide( const_cast<Surface*>(&(i->surface().surface())), 
                     i->globalFace(), i->surfaceSide()));
     //  LogDebug("NavGeometry") << " or actually this is where we have side " << i->surfaceSide() << " and face " << i->globalFace() ;
   }
 
 
 
     computeBounds(faces);
 }

NavVolume6Faces::NavVolume6Faces	(	const MagVolume &	magvol,
		const bool	isIron = `false`
	)

explicit

A NavVolume6Faces that corresponds exactly to a MagVolume.

Definition at line 44 of file NavVolume6Faces.cc.

References computeBounds(), MagVolume::faces(), and mps_fire::i.

                                                                       :
   NavVolume( magvol.position(), magvol.rotation(), magvol.shapeType(), magvol.provider()),
   theFaces(magvol.faces()), isThisIron(isIron)
 {
   std::vector<NavVolumeSide> navSides;
   std::vector<VolumeSide> magSides( magvol.faces());
   NavSurfaceBuilder navBuilder;
 
   for (std::vector<VolumeSide>::const_iterator i=magSides.begin();
        i != magSides.end(); i++) {
     NavSurface* navSurface = navBuilder.build( i->surface());
     navSides.push_back( NavVolumeSide( navSurface, i->globalFace(), i->surfaceSide()));
   }
   computeBounds(navSides);
 }

Member Function Documentation

void NavVolume6Faces::computeBounds ( const std::vector< NavVolumeSide > & faces )

private

Definition at line 73 of file NavVolume6Faces.cc.

References NavSurface::addVolume(), NavSurface::bounds(), mps_fire::i, and theNavSurfaces.

Referenced by NavVolume6Faces().

 {
     bool allPlanes = true;
     // bool allPlanes = false; // for TESTS ONLY!!!
     std::vector<const Plane*> planes;
     for (std::vector<NavVolumeSide>::const_iterator iface=faces.begin(); iface!=faces.end(); iface++) {
     const Plane* plane = dynamic_cast<const Plane*>(&(iface->surface()));
     if (plane != nullptr) {
         planes.push_back(plane);
     }
     else allPlanes = false;
     }
     
     for (unsigned int i=0; i<faces.size(); i++) {
 
       // FIXME: who owns the new NavSurface? memory leak???
 
     NavSurface& navSurf = faces[i].mutableSurface();
     Bounds* myBounds = nullptr;
     if (allPlanes) {    
         myBounds = computeBounds( i, planes);
     }
     else {
         myBounds = computeBounds( i, faces);
     }
     navSurf.addVolume( this, myBounds, faces[i].surfaceSide());
     delete myBounds; // since navSurf now owns a copy
 
 // this is tricky: we want to avoid multiple copies of the Bounds; the NavSurface owns
 // a copy of Bounds for each touching Volume (instantiated in the call to addVolume).
 // We would like to keep a pointer to the same Bounds in the NavVolume, so we have to ASK
 // the NavSurface for the Bounds* of the Bounds we just gave it!
     //LogDebug("NavGeometry") << "Adding a Volume Side with center " << navSurf.surface().position() << " side "<< faces[i].surfaceSide() << " and face " << faces[i].globalFace();
     theNavSurfaces.push_back( SurfaceAndBounds(&navSurf, navSurf.bounds(this), faces[i].surfaceSide(), faces[i].globalFace()));
     }
 }

Bounds * NavVolume6Faces::computeBounds	(	int	index,
		const std::vector< const Plane * > &	bpc
	)

private

Definition at line 110 of file NavVolume6Faces.cc.

References TCMET_cfi::corner, gather_cfg::cout, ThreePlaneCrossing::crossing(), MillePedeFileConverter_cfg::e, mps_fire::i, Plane::localZ(), Plane::normalVector(), SurfaceOrientation::onSurface, GloballyPositioned< T >::position(), Plane::side(), and GloballyPositioned< T >::toLocal().

 {
   const Plane* plane( bpc[index]);
 
   // find the 4 intersecting planes
   int startIndex = 2*(1+index/2); // 2, 4, 6
   std::vector<const Plane*> crossed; crossed.reserve(4);
   for (int j = startIndex; j <  startIndex+4; j++) {
     crossed.push_back(bpc[j%6]);
   }
 
   // compute intersection corners of the plane triplets
   std::vector<GlobalPoint> corners; corners.reserve(4);
   ThreePlaneCrossing crossing;
   for ( int i=0; i<2; i++) {
     for ( int j=2; j<4; j++) {
       GlobalPoint corner( crossing.crossing( *plane, *crossed[i], *crossed[j]));
       corners.push_back(corner);
 
 #ifdef DEBUG
       std::cout << "Crossing of planes is " << corner << std::endl;
       std::cout << "NormalVectors of the planes are " << plane->normalVector()
        << " " << crossed[i]->normalVector() << " " << crossed[j]->normalVector() << std::endl;
       std::cout << "Positions of planes are " << plane->position()
        << " " << crossed[i]->position() << " " << crossed[j]->position() << std::endl;
       if (plane->side( corner, 1.e-5) == SurfaceOrientation::onSurface &&
       crossed[i]->side( corner, 1.e-5) == SurfaceOrientation::onSurface &&
       crossed[j]->side( corner, 1.e-5) == SurfaceOrientation::onSurface) {
       std::cout << "Crossing is really on all three surfaces" << std::endl;
       }
       else {
       std::cout << "CROSSING IS NOT ON SURFACES!!!" << std::endl;
       std::cout << plane->localZ(corner) << std::endl;
       std::cout << crossed[i]->localZ(corner) << std::endl;
       std::cout << crossed[j]->localZ(corner) << std::endl;
        }
 #endif
 
     }
   }
 
   // put corners in cyclic sequence (2 and 3 swapped)
   return new FourPointPlaneBounds( plane->toLocal( corners[0]), plane->toLocal( corners[1]),
                    plane->toLocal( corners[3]), plane->toLocal( corners[2]));
 }

Bounds * NavVolume6Faces::computeBounds	(	int	index,
		const std::vector< NavVolumeSide > &	faces
	)

private

Definition at line 157 of file NavVolume6Faces.cc.

References NavVolumeSide::surface(), NavSurface::surface(), and NavVolumeSide::surfaceSide().

 {
     typedef GeneralNSurfaceDelimitedBounds::SurfaceAndSide         SurfaceAndSide;
     typedef GeneralNSurfaceDelimitedBounds::SurfaceContainer       SurfaceContainer;
 
   // find the 4 intersecting surfaces
   int startIndex = 2*(1+index/2); // 2, 4, 6
   SurfaceContainer crossed; crossed.reserve(4);
   for (int j = startIndex; j <  startIndex+4; j++) {
     const NavVolumeSide& face(faces[j%6]);
     crossed.push_back( SurfaceAndSide(&(face.surface().surface()), face.surfaceSide()));
   }
   return new GeneralNSurfaceDelimitedBounds( &(faces[index].surface().surface()), crossed);
 }

VolumeCrossReturnType NavVolume6Faces::crossToNextVolume	(	const TrajectoryStateOnSurface &	currentState,
		const Propagator &	prop
	)		const

overridevirtual

Cross this volume and point at the next.

Implements NavVolume.

Definition at line 250 of file NavVolume6Faces.cc.

References alongMomentum, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), LogDebug, nextSurface(), oppositeSide(), TrajectoryStateOnSurface::surface(), and GloballyPositioned< float >::toLocal().

 {
   typedef TrajectoryStateOnSurface TSOS;
   typedef std::pair<TSOS,double> TSOSwithPath;
 
   NavVolume::Container nsc = nextSurface( toLocal( startingState.globalPosition()), 
                       toLocal( startingState.globalMomentum()), -1,
                       alongMomentum, &(startingState.surface()));
   int itry = 0;
   VolumeCrossReturnType VolumeCrossResult( nullptr, startingState, 0.0);
 
   for (NavVolume::Container::const_iterator isur = nsc.begin(); isur!=nsc.end(); isur++) {
 
     LogDebug("NavGeometry") <<  "crossToNextVolume: trying Surface no. " << itry ;
     TSOSwithPath state;
     
     try {
       state = isur->surface().propagateWithPath( prop, startingState);
     }
     catch (MagVolumeOutsideValidity& except) {
       LogDebug("NavGeometry") << "Ohoh... failed to stay inside magnetic field !! skip this surface " ;
     ++itry;
       continue;
     }
     
     if (!state.first.isValid()) {
       ++itry;
       continue;
     }
 
     LogDebug("NavGeometry") <<  "crossToNextVolume: reached Valid State at Surface no. " << itry ;
     LogDebug("NavGeometry") << " --> local position of Valid state is " <<  state.first.localPosition()  ;
     LogDebug("NavGeometry") << " --> global position of Valid state is " <<  state.first.globalPosition()  ;
 
     if (isur->bounds().inside(state.first.localPosition())) {
       //LogDebug("NavGeometry") << "crossToNextVolume: Surface containing destination point found at try " << itry ;
       // Found the exit surface !! Get pointer to next volume and save exit state:
       //VolumeCrossResult.first = isur->surface().nextVolume(state.localPosition(),oppositeSide(isur->side()));
       //VolumeCrossResult.second = state;
       //      exitSurface = &( isur->surface().surface() );
       //if(VolumeCrossResult.path() < 0.01) {
       //    LogDebug("NavGeometry") << " Stuck at  " << state.first.globalPosition() ;
       //}
       return VolumeCrossReturnType ( isur->surface().nextVolume(state.first.localPosition(), oppositeSide(isur->side())),
                 state.first, state.second );
       
       break;
     }
     else {
       LogDebug("NavGeometry") << "crossToNextVolume: BUT not inside the Bounds !! " ;
       ++itry;
     }
   }
 
   return VolumeCrossResult;
 }

const std::vector<VolumeSide>& NavVolume6Faces::faces ( ) const

inlineoverridevirtual

Access to volume faces.

Implements MagVolume.

Definition at line 39 of file NavVolume6Faces.h.

References theFaces.

39 { return theFaces; }

NavVolume6Faces::theFaces

std::vector< VolumeSide > theFaces

Definition: NavVolume6Faces.h:49

bool NavVolume6Faces::inside	(	const GlobalPoint &	gp,
		double	tolerance
	)		const

overridevirtual

Implements MagVolume.

Definition at line 61 of file NavVolume6Faces.cc.

References mps_fire::i, SurfaceOrientation::onSurface, and theFaces.

Referenced by isIron().

 {
   // check if the point is on the correct side of all delimiting surfaces
   for (std::vector<VolumeSide>::const_iterator i=theFaces.begin(); i!=theFaces.end(); i++) {
     Surface::Side side = i->surface().side( gp, tolerance);
     if ( side != i->surfaceSide() && side != SurfaceOrientation::onSurface) return false;
   }
   return true;
 }

bool NavVolume6Faces::isIron ( ) const

inlineoverridevirtual

Access to Iron/Air information:

Implements NavVolume.

Definition at line 42 of file NavVolume6Faces.h.

References runTauDisplay::gp, MagVolume::inside(), inside(), isThisIron, and RecoTauCleanerPlugins::tolerance.

42 { return isThisIron; }

NavVolume6Faces::isThisIron

bool isThisIron

Definition: NavVolume6Faces.h:51

NavVolume::Container NavVolume6Faces::nextSurface	(	const NavVolume::LocalPoint &	pos,
		const NavVolume::LocalVector &	mom,
		double	charge,
		PropagationDirection	propDir = `alongMomentum`
	)		const

overridevirtual

Give a sorted list of possible surfaces to propagate to.

Implements NavVolume.

Definition at line 174 of file NavVolume6Faces.cc.

References alongMomentum, geometryDiff::epsilon, mps_fire::i, mps_fire::result, theNavSurfaces, and GloballyPositioned< float >::toGlobal().

Referenced by crossToNextVolume().

 {
     typedef std::map<double,SurfaceAndBounds>   SortedContainer;
 
     GlobalPoint  gpos( toGlobal(pos));
     GlobalVector gmom( toGlobal(mom));
     GlobalVector gdir = (propDir == alongMomentum ? gmom : -gmom);
 
     SortedContainer sortedSurfaces;
     Container       verycloseSurfaces; // reachable surface with dist < epsilon !!
     Container       unreachableSurfaces;
 
     double epsilon = 0.01; // should epsilon be hard-coded or a variable in NavVolume?
 
     for (Container::const_iterator i=theNavSurfaces.begin(); i!=theNavSurfaces.end(); i++) {
     std::pair<bool,double> dist = i->surface().distanceAlongLine( gpos, gdir);
     if (dist.first) {
       if (dist.second > epsilon) sortedSurfaces[dist.second] = *i;
       else verycloseSurfaces.push_back(*i);
     } 
     else unreachableSurfaces.push_back(*i);
     }
     NavVolume::Container result;
     for (SortedContainer::const_iterator i=sortedSurfaces.begin(); i!=sortedSurfaces.end(); ++i) {
     result.push_back(i->second);
     }
     result.insert( result.end(), unreachableSurfaces.begin(), unreachableSurfaces.end());
     result.insert( result.end(), verycloseSurfaces.begin(), verycloseSurfaces.end());
     return result;
 }

NavVolume::Container NavVolume6Faces::nextSurface	(	const NavVolume::LocalPoint &	pos,
		const NavVolume::LocalVector &	mom,
		double	charge,
		PropagationDirection	propDir,
		ConstReferenceCountingPointer< Surface >	NotThisSurfaceP
	)		const

overridevirtual

Same, giving lowest priority to the surface we are on now (=NotThisSurface)

Implements NavVolume.

Definition at line 208 of file NavVolume6Faces.cc.

References alongMomentum, geometryDiff::epsilon, mps_fire::i, mps_fire::result, theNavSurfaces, and GloballyPositioned< float >::toGlobal().

 {
     typedef std::map<double,SurfaceAndBounds>   SortedContainer;
 
     GlobalPoint  gpos( toGlobal(pos));
     GlobalVector gmom( toGlobal(mom));
     GlobalVector gdir = (propDir == alongMomentum ? gmom : -gmom);
 
     SortedContainer sortedSurfaces;
     Container       verycloseSurfaces; // reachable surface with dist < epsilon (if 6-surface check fails)
     Container       unreachableSurfaces;
 
     double epsilon = 0.01; // should epsilon be hard-coded or a variable in NavVolume?
     bool surfaceMatched = false;
 
     for (Container::const_iterator i=theNavSurfaces.begin(); i!=theNavSurfaces.end(); i++) {
       if (&(i->surface().surface()) == NotThisSurfaceP) surfaceMatched = true;
     }
     
     for (Container::const_iterator i=theNavSurfaces.begin(); i!=theNavSurfaces.end(); i++) {
     std::pair<bool,double> dist = i->surface().distanceAlongLine( gpos, gdir);
     if (dist.first) { 
       if ( &(i->surface().surface()) == NotThisSurfaceP || (!surfaceMatched && dist.second < epsilon)) 
         verycloseSurfaces.push_back(*i);
       else sortedSurfaces[dist.second] = *i;
     }
     else unreachableSurfaces.push_back(*i);
     }
 
     NavVolume::Container result;
     for (SortedContainer::const_iterator i=sortedSurfaces.begin(); i!=sortedSurfaces.end(); ++i) {
     result.push_back(i->second);
     }
     result.insert( result.end(), unreachableSurfaces.begin(), unreachableSurfaces.end());
     result.insert( result.end(), verycloseSurfaces.begin(), verycloseSurfaces.end());
     return result;
 }

Member Data Documentation

bool NavVolume6Faces::isThisIron

private

Definition at line 51 of file NavVolume6Faces.h.

Referenced by isIron().

std::vector<VolumeSide> NavVolume6Faces::theFaces

private

Definition at line 49 of file NavVolume6Faces.h.

Referenced by faces(), inside(), and NavVolume6Faces().

Container NavVolume6Faces::theNavSurfaces

private

Definition at line 50 of file NavVolume6Faces.h.

Referenced by computeBounds(), and nextSurface().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation