CylinderBuilderFromDet Class Reference

#include <CylinderBuilderFromDet.h>

Public Types
typedef GeomDet	Det
typedef Surface::PositionType	PositionType
typedef Surface::RotationType	RotationType
typedef PositionType::BasicVectorType	Vector

Public Member Functions
BoundCylinder *	build () const
	CylinderBuilderFromDet ()
BoundCylinder *	operator() (std::vector< const Det * >::const_iterator first, std::vector< const Det * >::const_iterator last) const
void	operator() (const Det &det)

Private Attributes
float	rmax
float	rmin
float	zmax
float	zmin

Detailed Description

Given a container of GeomDets, constructs a cylinder of minimal dimensions that contains all of the Dets completely (all corners etc.) Useful for defining a BarrelDetLayer from a group of DetUnits.

Definition at line 18 of file CylinderBuilderFromDet.h.

Member Typedef Documentation

typedef GeomDet CylinderBuilderFromDet::Det

Definition at line 20 of file CylinderBuilderFromDet.h.

typedef Surface::PositionType CylinderBuilderFromDet::PositionType

Definition at line 21 of file CylinderBuilderFromDet.h.

typedef Surface::RotationType CylinderBuilderFromDet::RotationType

Definition at line 22 of file CylinderBuilderFromDet.h.

typedef PositionType::BasicVectorType CylinderBuilderFromDet::Vector

Definition at line 23 of file CylinderBuilderFromDet.h.

Constructor & Destructor Documentation

CylinderBuilderFromDet::CylinderBuilderFromDet ( )

inline

Definition at line 25 of file CylinderBuilderFromDet.h.

      : rmin(std::numeric_limits<float>::max()),
        rmax(0.0),
        zmin(std::numeric_limits<float>::max()),
        zmax(std::numeric_limits<float>::min()) {}

Member Function Documentation

BoundCylinder * CylinderBuilderFromDet::build

(

)

const

Definition at line 78 of file CylinderBuilderFromDet.cc.

References Cylinder::computeRadius(), makeMuonMisalignmentScenario::rot, and PV3DBase< T, PVType, FrameType >::z().

                                                   {
  PositionType pos(0, 0, 0.5 * (zmin + zmax));
  RotationType rot;  // only "barrel" orientation supported

  auto scp = new SimpleCylinderBounds(rmin, rmax, zmin - pos.z(), zmax - pos.z());
  return new Cylinder(Cylinder::computeRadius(*scp), pos, rot, scp);
}

BoundCylinder * CylinderBuilderFromDet::operator()	(	std::vector< const Det * >::const_iterator	first,
		std::vector< const Det * >::const_iterator	last
	)		const

Definition at line 8 of file CylinderBuilderFromDet.cc.

References Cylinder::computeRadius(), BoundingBox::corners(), mps_fire::i, dqmdumpme::last, SiStripPI::max, min(), Basic3DVector< T >::perp(), alignCSCRings::r, makeMuonMisalignmentScenario::rot, and PV3DBase< T, PVType, FrameType >::z().

                                                                                               {
  // find mean position and radius
  typedef PositionType::BasicVectorType Vector;
  Vector posSum(0, 0, 0);
  float rSum = 0;
  for (vector<const Det*>::const_iterator i = first; i != last; i++) {
    posSum += (**i).surface().position().basicVector();
    rSum += (**i).surface().position().perp();
  }
  float div(1 / float(last - first));
  PositionType meanPos(div * posSum);
  float meanR(div * rSum);

  // find max deviations from mean pos in Z and from mean R
  float rmin = meanR;
  float rmax = meanR;
  float zmin = meanPos.z();
  float zmax = meanPos.z();
  for (vector<const Det*>::const_iterator i = first; i != last; i++) {
    vector<GlobalPoint> corners = BoundingBox::corners(dynamic_cast<const Plane&>((**i).surface()));
    for (vector<GlobalPoint>::const_iterator ic = corners.begin(); ic != corners.end(); ic++) {
      float r = ic->perp();
      float z = ic->z();
      rmin = min(rmin, r);
      rmax = max(rmax, r);
      zmin = min(zmin, z);
      zmax = max(zmax, z);
    }
    // in addition to the corners we have to check the middle of the
    // det +/- thickness/2
    // , since the min  radius for some barrel dets is reached there
    float rdet = (**i).surface().position().perp();
    float halfThick = (**i).surface().bounds().thickness() / 2.F;
    rmin = min(rmin, rdet - halfThick);
    rmax = max(rmax, rdet + halfThick);
  }

  // the transverse position is zero by construction.
  // the Z position is the average between zmin and zmax, since the bounds
  // are symmetric
  // for the same reason the R is the average between rmin and rmax,
  // but this is done by the Bounds anyway.

  PositionType pos(0, 0, 0.5 * (zmin + zmax));
  RotationType rot;  // only "barrel" orientation supported

  auto scp = new SimpleCylinderBounds(rmin, rmax, zmin - pos.z(), zmax - pos.z());
  return new Cylinder(Cylinder::computeRadius(*scp), pos, rot, scp);
}

void CylinderBuilderFromDet::operator()

(

const Det &

det

)

Definition at line 59 of file CylinderBuilderFromDet.cc.

References Surface::bounds(), SiStripPI::max, min(), PV3DBase< T, PVType, FrameType >::perp(), GloballyPositioned< T >::position(), alignCSCRings::r, GeomDet::surface(), and Bounds::thickness().

                                                      {
  BoundingBox bb(dynamic_cast<const Plane&>(det.surface()));
  for (int nc = 0; nc < 8; ++nc) {
    float r = bb[nc].perp();
    float z = bb[nc].z();
    rmin = std::min(rmin, r);
    rmax = std::max(rmax, r);
    zmin = std::min(zmin, z);
    zmax = std::max(zmax, z);
  }
  // in addition to the corners we have to check the middle of the
  // det +/- thickness/2
  // , since the min  radius for some barrel dets is reached there
  float rdet = det.surface().position().perp();
  float halfThick = det.surface().bounds().thickness() / 2.F;
  rmin = std::min(rmin, rdet - halfThick);
  rmax = std::max(rmax, rdet + halfThick);
}

Member Data Documentation

float CylinderBuilderFromDet::rmax

private

Definition at line 40 of file CylinderBuilderFromDet.h.

float CylinderBuilderFromDet::rmin

private

Definition at line 39 of file CylinderBuilderFromDet.h.

float CylinderBuilderFromDet::zmax

private

Definition at line 42 of file CylinderBuilderFromDet.h.

float CylinderBuilderFromDet::zmin

private

Definition at line 41 of file CylinderBuilderFromDet.h.