MagGeoBuilderFromDDD::volumeHandle Class Reference

#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 (double halfX, double halfY, double halfZ)
void	buildCons (double zhalf, double rInMinusZ, double rOutMinusZ, double rInPlusZ, double rOutPlusZ, double startPhi, double deltaPhi)
void	buildPseudoTrap (double x1, double x2, double y1, double y2, double halfZ, double radius, bool atMinusZ)
void	buildTrap (double x1, double x2, double x3, double x4, double y1, double y2, double theta, double phi, double halfZ, double alpha1, double alpha2)
void	buildTruncTubs (double zhalf, double rIn, double rOut, double startPhi, double deltaPhi, double cutAtStart, double cutAtDelta, bool cutInside)
void	buildTubs (double zhalf, double rIn, double rOut, double startPhi, double deltaPhi)
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 40 of file volumeHandle.cc.

References DDTrap::alpha1(), DDTrap::alpha2(), DDPseudoTrap::atMinusZ(), buildBox(), buildCons(), buildPseudoTrap(), buildTrap(), buildTruncTubs(), buildTubs(), magneticfield::BaseVolumeHandle::center(), magneticfield::BaseVolumeHandle::center_, magneticfield::BaseVolumeHandle::copyno, DDExpandedView::copyno(), gather_cfg::cout, DDTruncTubs::cutAtDelta(), DDTruncTubs::cutAtStart(), DDTruncTubs::cutInside(), ddbox, ddcons, ddpseudotrap, ddtrap, ddtrunctubs, ddtubs, magneticfield::BaseVolumeHandle::debug, srCondWrite_cfg::deltaPhi, DDTruncTubs::deltaPhi(), DDTubs::deltaPhi(), DDCons::deltaPhi(), DDBox::halfX(), DDBox::halfY(), DDTrap::halfZ(), DDPseudoTrap::halfZ(), DDBox::halfZ(), mps_fire::i, magneticfield::BaseVolumeHandle::isAssigned, magneticfield::BaseVolumeHandle::isIronFlag, DDExpandedView::logicalPart(), magneticfield::BaseVolumeHandle::magFile, magneticfield::BaseVolumeHandle::masterSector, DDLogicalPart::material(), DDSolidShapesName::name(), DDName::name(), DDBase< N, C >::name(), magneticfield::BaseVolumeHandle::name, PV3DBase< T, PVType, FrameType >::perp(), phi, PV3DBase< T, PVType, FrameType >::phi(), DDTrap::phi(), DDCons::phiFrom(), DDPseudoTrap::radius(), CosmicsPD_Skims::radius, referencePlane(), DDTruncTubs::rIn(), DDTubs::rIn(), DDCons::rInMinusZ(), DDCons::rInPlusZ(), DDTruncTubs::rOut(), DDTubs::rOut(), DDCons::rOutMinusZ(), DDCons::rOutPlusZ(), shape(), DDSolid::shape(), solid, DDLogicalPart::solid(), DDTruncTubs::startPhi(), DDTubs::startPhi(), AlCaHLTBitMon_QueryRunRegistry::string, magneticfield::BaseVolumeHandle::surfaces, magneticfield::BaseVolumeHandle::theRMax, magneticfield::BaseVolumeHandle::theRMin, magneticfield::BaseVolumeHandle::theRN, DDTrap::theta(), theta(), DDExpandedView::translation(), magneticfield::BaseVolumeHandle::volumeno, DDTrap::x1(), DDPseudoTrap::x1(), DDTrap::x2(), DDPseudoTrap::x2(), DDTrap::x3(), DDTrap::x4(), DDTrap::y1(), DDPseudoTrap::y1(), DDTrap::y2(), DDPseudoTrap::y2(), DDTruncTubs::zHalf(), DDTubs::zhalf(), and DDCons::zhalf().

    : 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 = std::stoul(volName);

  for (int i = 0; i < 6; ++i) {
    isAssigned[i] = false;
  }

  if (debug) {
    cout.precision(7);
  }

  referencePlane(fv);

  if (solid.shape() == DDSolidShape::ddbox) {
    DDBox box(solid);
    double halfX = convertUnits(box.halfX());
    double halfY = convertUnits(box.halfY());
    double halfZ = convertUnits(box.halfZ());
    buildBox(halfX, halfY, halfZ);
  } else if (solid.shape() == DDSolidShape::ddtrap) {
    DDTrap trap(solid);
    double x1 = convertUnits(trap.x1());
    double x2 = convertUnits(trap.x2());
    double x3 = convertUnits(trap.x3());
    double x4 = convertUnits(trap.x4());
    double y1 = convertUnits(trap.y1());
    double y2 = convertUnits(trap.y2());
    double theta = trap.theta();
    double phi = trap.phi();
    double halfZ = convertUnits(trap.halfZ());
    double alpha1 = trap.alpha1();
    double alpha2 = trap.alpha2();
    buildTrap(x1, x2, x3, x4, y1, y2, theta, phi, halfZ, alpha1, alpha2);
  } else if (solid.shape() == DDSolidShape::ddcons) {
    DDCons cons(solid);
    double zhalf = convertUnits(cons.zhalf());
    double rInMinusZ = convertUnits(cons.rInMinusZ());
    double rOutMinusZ = convertUnits(cons.rOutMinusZ());
    double rInPlusZ = convertUnits(cons.rInPlusZ());
    double rOutPlusZ = convertUnits(cons.rOutPlusZ());
    double startPhi = cons.phiFrom();
    double deltaPhi = cons.deltaPhi();
    buildCons(zhalf, rInMinusZ, rOutMinusZ, rInPlusZ, rOutPlusZ, startPhi, deltaPhi);
  } else if (solid.shape() == DDSolidShape::ddtubs) {
    DDTubs tubs(solid);
    double zhalf = convertUnits(tubs.zhalf());
    double rIn = convertUnits(tubs.rIn());
    double rOut = convertUnits(tubs.rOut());
    double startPhi = tubs.startPhi();
    double deltaPhi = tubs.deltaPhi();
    buildTubs(zhalf, rIn, rOut, startPhi, deltaPhi);
  } else if (solid.shape() == DDSolidShape::ddpseudotrap) {
    DDPseudoTrap ptrap(solid);
    double x1 = convertUnits(ptrap.x1());
    double x2 = convertUnits(ptrap.x2());
    double y1 = convertUnits(ptrap.y1());
    double y2 = convertUnits(ptrap.y2());
    double halfZ = convertUnits(ptrap.halfZ());
    double radius = convertUnits(ptrap.radius());
    bool atMinusZ = ptrap.atMinusZ();
    buildPseudoTrap(x1, x2, y1, y2, halfZ, radius, atMinusZ);
  } else if (solid.shape() == DDSolidShape::ddtrunctubs) {
    DDTruncTubs tubs(solid);
    double zhalf = convertUnits(tubs.zHalf());            // half of the z-Axis
    double rIn = convertUnits(tubs.rIn());                // inner radius
    double rOut = convertUnits(tubs.rOut());              // outer radius
    double startPhi = tubs.startPhi();                    // angular start of the tube-section
    double deltaPhi = tubs.deltaPhi();                    // angular span of the tube-section
    double cutAtStart = convertUnits(tubs.cutAtStart());  // truncation at begin of the tube-section
    double cutAtDelta = convertUnits(tubs.cutAtDelta());  // truncation at end of the tube-section
    bool cutInside = tubs.cutInside();  // true, if truncation is on the inner side of the tube-section
    buildTruncTubs(zhalf, rIn, rOut, startPhi, deltaPhi, cutAtStart, cutAtDelta, cutInside);
  } 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 ( double  halfX, double  halfY, double  halfZ  )

private

Referenced by volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildCons ( double  zhalf, double  rInMinusZ, double  rOutMinusZ, double  rInPlusZ, double  rOutPlusZ, double  startPhi, double  deltaPhi  )

private

Referenced by volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildPseudoTrap ( double  x1, double  x2, double  y1, double  y2, double  halfZ, double  radius, bool  atMinusZ  )

private

Referenced by volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildTrap ( double  x1, double  x2, double  x3, double  x4, double  y1, double  y2, double  theta, double  phi, double  halfZ, double  alpha1, double  alpha2  )

private

Referenced by volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildTruncTubs ( double  zhalf, double  rIn, double  rOut, double  startPhi, double  deltaPhi, double  cutAtStart, double  cutAtDelta, bool  cutInside  )

private

Referenced by volumeHandle().

void MagGeoBuilderFromDDD::volumeHandle::buildTubs ( double  zhalf, double  rIn, double  rOut, double  startPhi, double  deltaPhi  )

private

Referenced by volumeHandle().

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

delete

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

private

Definition at line 201 of file volumeHandle.cc.

References gather_cfg::cout, ddpseudotrap, MagGeoBuilderFromDDD::debug, Vector3DBase< T, FrameTag >::dot(), DDExpandedView::rotation(), x, y, and z.

Referenced by 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 DDSolid::shape(), and solid.

Referenced by volumeHandle().

{ return solid.shape(); }

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 274 of file volumeHandle.cc.

References ddtubs, mps_fire::i, SurfaceOrientation::inner, SurfaceOrientation::phiminus, SurfaceOrientation::phiplus, mps_fire::result, and alignCSCRings::s.

Referenced by MagGeoBuilderFromDDD::buildInterpolator().

                                                                    {
  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 78 of file volumeHandle.h.

Referenced by shape(), and volumeHandle().