#include <volumeHandle.h>

Inheritance diagram for MagGeoBuilderFromDDD::volumeHandle:

Public Member Functions
volumeHandle	operator= (const volumeHandle &v)=delete

DDSolidShape	shape () const override
	Shape of the solid. More...

std::vector< VolumeSide >	sides () const override
	The surfaces and they orientation, as required to build a MagVolume. More...

	volumeHandle (const DDExpandedView &fv, bool expand2Pi=false, bool debugVal=false)

	volumeHandle (const volumeHandle &v)=delete

Public Member Functions inherited from magneticfield::BaseVolumeHandle
	BaseVolumeHandle (bool expand2Pi=false, bool debugVal=false)

	BaseVolumeHandle (const BaseVolumeHandle &v)=delete

const GlobalPoint &	center () const
	Return the center of the volume. More...

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

bool	isPlaneMatched (int which_side) const
	if the specified surface has been matched. More...

Geom::Phi< float >	maxPhi () const
	Maximum value of phi covered by the volume. More...

double	maxZ () const

Geom::Phi< float >	minPhi () const
	Minimum value of phi covered by the volume. More...

double	minR () const
	Minimum R for any point within the volume. More...

double	minZ () const
	Z limits. More...

const GloballyPositioned< float > *	placement () const
	Position and rotation. More...

int	references (int which_side) const

const double	RN () const
	Distance of (x,y) plane from origin. More...

bool	sameSurface (const Surface &s1, Sides which_side, float tolerance=0.01)
	Find out if two surfaces are the same physical surface. More...

bool	setSurface (const Surface &s1, Sides which_side)
	Assign a shared surface perorming sanity checks. More...

const Surface &	surface (int which_side) const
	Get the current surface on specified side. More...

const Surface &	surface (Sides which_side) const

bool	toExpand () const

virtual	~BaseVolumeHandle ()

Private Member Functions
void	buildBox ()

void	buildCons ()

void	buildPseudoTrap ()

void	buildTrap ()

void	buildTruncTubs ()

void	buildTubs ()

void	referencePlane (const DDExpandedView &fv)

Private Attributes
DDSolid	solid

Additional Inherited Members
Public Types inherited from magneticfield::BaseVolumeHandle
typedef Surface::GlobalPoint	GlobalPoint

typedef Surface::LocalPoint	LocalPoint

typedef Surface::LocalVector	LocalVector

typedef SurfaceOrientation::GlobalFace	Sides

Public Attributes inherited from magneticfield::BaseVolumeHandle
unsigned short	copyno
	copy number More...

std::string	magFile
	Name of magnetic field table file. More...

MagVolume6Faces *	magVolume
	Pointer to the final MagVolume (must be set from outside) More...

int	masterSector
	The sector for which an interpolator for this class of volumes should be built. More...

std::string	name
	Name of the volume. More...

unsigned short	volumeno
	volume number More...

Protected Types inherited from magneticfield::BaseVolumeHandle
typedef ConstReferenceCountingPointer< Surface >	RCPS

Protected Member Functions inherited from magneticfield::BaseVolumeHandle
void	buildPhiZSurf (double startPhi, double deltaPhi, double zhalf, double rCentr)

Protected Attributes inherited from magneticfield::BaseVolumeHandle
GlobalPoint	center_

const bool	debug

bool	expand

bool	isAssigned [6]

bool	isIronFlag

GloballyPositioned< float > *	refPlane

RCPS	surfaces [6]

Geom::Phi< float >	thePhiMin

double	theRMax

double	theRMin

double	theRN

Detailed Description

A temporary container to cache info on a six-surface volume during the processing. Used to sort, organise, and build shared planes. One instance is created for each DDVolume. The parameters of the boundary surfaces are calculated during construction.

Author: N. Amapane - INFN Torino (original developer)

Definition at line 23 of file volumeHandle.h.

Constructor & Destructor Documentation

MagGeoBuilderFromDDD::volumeHandle::volumeHandle	(	const DDExpandedView &	fv,
		bool	expand2Pi = `false`,
		bool	debugVal = `false`
	)

Definition at line 32 of file volumeHandle.cc.

     : magneticfield::BaseVolumeHandle(expand2Pi, debugVal) {
   name = fv.logicalPart().name().name();
   copyno = fv.copyno();
   solid = fv.logicalPart().solid();
   center_ = GlobalPoint(fv.translation().x() / cm, fv.translation().y() / cm, fv.translation().z() / cm);
 
   // ASSUMPTION: volume names ends with "_NUM" where NUM is the volume number
   string volName = name;
   volName.erase(0, volName.rfind('_') + 1);
   volumeno = boost::lexical_cast<unsigned short>(volName);
 
   for (int i = 0; i < 6; ++i) {
     isAssigned[i] = false;
   }
 
   if (debug) {
     cout.precision(7);
   }
 
   referencePlane(fv);
 
   if (solid.shape() == DDSolidShape::ddbox) {
     buildBox();
   } else if (solid.shape() == DDSolidShape::ddtrap) {
     buildTrap();
   } else if (solid.shape() == DDSolidShape::ddcons) {
     buildCons();
   } else if (solid.shape() == DDSolidShape::ddtubs) {
     buildTubs();
   } else if (solid.shape() == DDSolidShape::ddpseudotrap) {
     buildPseudoTrap();
   } else if (solid.shape() == DDSolidShape::ddtrunctubs) {
     buildTruncTubs();
   } else {
     cout << "volumeHandle ctor: Unexpected solid: " << DDSolidShapesName::name(solid.shape()) << endl;
   }
 
   // NOTE: Table name and master sector are no longer taken from xml!
   //   DDsvalues_type sv(fv.mergedSpecifics());
 
   //   { // Extract the name of associated field file.
   //     std::vector<std::string> temp;
   //     std::string pname = "table";
   //     DDValue val(pname);
   //     DDsvalues_type sv(fv.mergedSpecifics());
   //     if (DDfetch(&sv,val)) {
   //       temp = val.strings();
   //       if (temp.size() != 1) {
   //    cout << "*** WARNING: volume has > 1 SpecPar " << pname << endl;
   //       }
   //       magFile = temp[0];
 
   //       string find="[copyNo]";
   //       std::size_t j;
   //       for ( ; (j = magFile.find(find)) != string::npos ; ) {
   //    stringstream conv;
   //    conv << setfill('0') << setw(2) << copyno;
   //    string repl;
   //    conv >> repl;
   //    magFile.replace(j, find.length(), repl);
   //       }
 
   //     } else {
   //       cout << "*** WARNING: volume does not have a SpecPar " << pname << endl;
   //       cout << " DDsvalues_type:  " << fv.mergedSpecifics() << endl;
   //     }
   //   }
 
   //   { // Extract the number of the master sector.
   //     std::vector<double> temp;
   //     const std::string pname = "masterSector";
   //     DDValue val(pname);
   //     if (DDfetch(&sv,val)) {
   //       temp = val.doubles();
   //       if (temp.size() != 1) {
   //    cout << "*** WARNING: volume has > 1 SpecPar " << pname << endl;
   //       }
   //       masterSector = int(temp[0]+.5);
   //     } else {
   //       if (MagGeoBuilderFromDDD::debug) {
   //    cout << "Volume does not have a SpecPar " << pname
   //         << " using: " << copyno << endl;
   //    cout << " DDsvalues_type:  " << fv.mergedSpecifics() << endl;
   //       }
   //       masterSector = copyno;
   //     }
   //   }
 
   // Get material for this volume
   if (fv.logicalPart().material().name().name() == "Iron")
     isIronFlag = true;
 
   if (debug) {
     cout << " RMin =  " << theRMin << endl;
     cout << " RMax =  " << theRMax << endl;
 
     if (theRMin < 0 || theRN < theRMin || theRMax < theRN)
       cout << "*** WARNING: wrong RMin/RN/RMax , shape: " << DDSolidShapesName::name(shape()) << endl;
 
     cout << "Summary: " << name << " " << copyno << " Shape= " << DDSolidShapesName::name(shape()) << " trasl "
          << center() << " R " << center().perp() << " phi " << center().phi() << " magFile " << magFile
          << " Material= " << fv.logicalPart().material().name() << " isIron= " << isIronFlag
          << " masterSector= " << masterSector << std::endl;
 
     cout << " Orientation of surfaces:";
     std::string sideName[3] = {"positiveSide", "negativeSide", "onSurface"};
     for (int i = 0; i < 6; ++i) {
       cout << "  " << i << ":" << sideName[surfaces[i]->side(center_, 0.3)];
     }
     cout << endl;
   }
 }

MagGeoBuilderFromDDD::volumeHandle::volumeHandle ( const volumeHandle & v )

delete

Member Function Documentation

void MagGeoBuilderFromDDD::volumeHandle::buildBox ( )

private

Referenced by shape(), and volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildCons ( )

private

Referenced by shape(), and volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildPseudoTrap ( )

private

Referenced by shape(), and volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildTrap ( )

private

Referenced by shape(), and volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildTruncTubs ( )

private

Referenced by shape(), and volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildTubs ( )

private

Referenced by shape(), and volumeHandle().

volumeHandle MagGeoBuilderFromDDD::volumeHandle::operator= ( const volumeHandle & v )

delete

void MagGeoBuilderFromDDD::volumeHandle::referencePlane ( const DDExpandedView & fv )

private

Definition at line 146 of file volumeHandle.cc.

References magneticfield::BaseVolumeHandle::center_, gather_cfg::cout, ddpseudotrap, magneticfield::BaseVolumeHandle::debug, GloballyPositioned< T >::position(), magneticfield::BaseVolumeHandle::refPlane, DDExpandedView::rotation(), DDSolid::shape(), solid, GloballyPositioned< T >::toGlobal(), x, y, z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by shape(), and volumeHandle().

                                                                               {
   // The refPlane is the "main plane" for the solid. It corresponds to the
   // x,y plane in the DDD local frame, and defines a frame where the local
   // coordinates are the same as in DDD.
   // In the geometry version 85l_030919, this plane is normal to the
   // beam line for all volumes but pseudotraps, so that global R is along Y,
   // global phi is along -X and global Z along Z:
   //
   //   Global(for vol at pi/2)    Local
   //   +R (+Y)                    +Y
   //   +phi(-X)                   -X
   //   +Z                         +Z
   //
   // For pseudotraps the refPlane is parallel to beam line and global R is
   // along Z, global phi is along +-X and and global Z along Y:
   //
   //   Global(for vol at pi/2)    Local
   //   +R (+Y)                    +Z
   //   +phi(-X)                   +X
   //   +Z                         +Y
   //
   // Note that the frame is centered in the DDD volume center, which is
   // inside the volume for DDD boxes and (pesudo)trapezoids, on the beam line
   // for tubs, cons and trunctubs.
 
   // In geometry version 1103l, trapezoids have X and Z in the opposite direction
   // than the above.  Boxes are either oriented as described above or in some case
   // have opposite direction for Y and X.
 
   // The global position
   Surface::PositionType &posResult = center_;
 
   // The reference plane rotation
   DD3Vector x, y, z;
   fv.rotation().GetComponents(x, y, z);
   if (debug) {
     if (x.Cross(y).Dot(z) < 0.5) {
       cout << "*** WARNING: Rotation is not RH " << endl;
     }
   }
 
   // The global rotation
   Surface::RotationType rotResult(float(x.X()),
                                   float(x.Y()),
                                   float(x.Z()),
                                   float(y.X()),
                                   float(y.Y()),
                                   float(y.Z()),
                                   float(z.X()),
                                   float(z.Y()),
                                   float(z.Z()));
 
   refPlane = new GloballyPositioned<float>(posResult, rotResult);
 
   // Check correct orientation
   if (debug) {
     cout << "Refplane pos  " << refPlane->position() << endl;
 
     // See comments above for the conventions for orientation.
     LocalVector globalZdir(0., 0., 1.);  // Local direction of the axis along global Z
     if (solid.shape() == DDSolidShape::ddpseudotrap) {
       globalZdir = LocalVector(0., 1., 0.);
     }
     if (refPlane->toGlobal(globalZdir).z() < 0.) {
       globalZdir = -globalZdir;
     }
 
     float chk = refPlane->toGlobal(globalZdir).dot(GlobalVector(0, 0, 1));
     if (chk < .999)
       cout << "*** WARNING RefPlane check failed!***" << chk << endl;
   }
 }

DDSolidShape MagGeoBuilderFromDDD::volumeHandle::shape ( ) const

inlineoverridevirtual

Shape of the solid.

Implements magneticfield::BaseVolumeHandle.

Definition at line 32 of file volumeHandle.h.

References buildBox(), buildCons(), buildPseudoTrap(), buildTrap(), buildTruncTubs(), buildTubs(), referencePlane(), DDSolid::shape(), sides(), and solid.

Referenced by volumeHandle().

32 { return solid.shape(); }

DDSolid::shape

DDSolidShape shape(void) const

The type of the solid.

Definition: DDSolid.cc:119

MagGeoBuilderFromDDD::volumeHandle::solid

DDSolid solid

Definition: volumeHandle.h:55

std::vector< VolumeSide > MagGeoBuilderFromDDD::volumeHandle::sides ( void ) const

overridevirtual

The surfaces and they orientation, as required to build a MagVolume.

Implements magneticfield::BaseVolumeHandle.

Definition at line 219 of file volumeHandle.cc.

References magneticfield::BaseVolumeHandle::center_, ddtubs, magneticfield::BaseVolumeHandle::expand, mps_fire::i, SurfaceOrientation::inner, SurfaceOrientation::phiminus, SurfaceOrientation::phiplus, mps_fire::result, alignCSCRings::s, DDSolid::shape(), solid, magneticfield::BaseVolumeHandle::surfaces, and magneticfield::BaseVolumeHandle::theRMin.

Referenced by MagGeoBuilderFromDDD::buildInterpolator(), and shape().

                                                                     {
   std::vector<VolumeSide> result;
   for (int i = 0; i < 6; ++i) {
     // If this is just a master volume out of wich a 2pi volume
     // should be built (e.g. central cylinder), skip the phi boundaries.
     if (expand && (i == phiplus || i == phiminus))
       continue;
 
     // FIXME: Skip null inner degenerate cylindrical surface
     if (solid.shape() == DDSolidShape::ddtubs && i == SurfaceOrientation::inner && theRMin < 0.001)
       continue;
 
     ReferenceCountingPointer<Surface> s = const_cast<Surface *>(surfaces[i].get());
     result.push_back(VolumeSide(s, GlobalFace(i), surfaces[i]->side(center_, 0.3)));
   }
   return result;
 }

Member Data Documentation

DDSolid MagGeoBuilderFromDDD::volumeHandle::solid

private

Definition at line 55 of file volumeHandle.h.

Referenced by referencePlane(), shape(), sides(), and volumeHandle().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation