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.

Date:

2007/03/07 16:28:39

Revision:

1.3

Definition at line 20 of file CylinderBuilderFromDet.h.

Member Typedef Documentation

typedef GeomDet CylinderBuilderFromDet::Det

Definition at line 22 of file CylinderBuilderFromDet.h.

typedef Surface::PositionType CylinderBuilderFromDet::PositionType

Definition at line 23 of file CylinderBuilderFromDet.h.

typedef Surface::RotationType CylinderBuilderFromDet::RotationType

Definition at line 24 of file CylinderBuilderFromDet.h.

typedef PositionType::BasicVectorType CylinderBuilderFromDet::Vector

Definition at line 25 of file CylinderBuilderFromDet.h.

Constructor & Destructor Documentation

CylinderBuilderFromDet::CylinderBuilderFromDet ( )

inline

Definition at line 27 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 83 of file CylinderBuilderFromDet.cc.

References pos, and PV3DBase< T, PVType, FrameType >::z().

                                                   {
  
  PositionType pos( 0, 0, 0.5*(zmin+zmax));
  RotationType rot;      // only "barrel" orientation supported
  
  return new BoundCylinder( pos, rot, 
                SimpleCylinderBounds( rmin, rmax, 
                          zmin-pos.z(), zmax-pos.z()));
}

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

Definition at line 9 of file CylinderBuilderFromDet.cc.

References BoundingBox::corners(), i, prof2calltree::last, max(), min, Basic3DVector< T >::perp(), pos, csvReporter::r, detailsBasic3DVector::z, 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 BoundPlane&>((**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
  
  return new BoundCylinder( pos, rot, 
                SimpleCylinderBounds( rmin, rmax, 
                          zmin-pos.z(), zmax-pos.z()));
}

void CylinderBuilderFromDet::operator()

(

const Det &

det

)

Definition at line 64 of file CylinderBuilderFromDet.cc.

References BoundSurface::bounds(), max(), min, PV3DBase< T, PVType, FrameType >::perp(), GloballyPositioned< T >::position(), csvReporter::r, GeomDet::surface(), Bounds::thickness(), and detailsBasic3DVector::z.

                                                      {
  BoundingBox bb( dynamic_cast<const BoundPlane&>(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 42 of file CylinderBuilderFromDet.h.

float CylinderBuilderFromDet::rmin

private

Definition at line 41 of file CylinderBuilderFromDet.h.

float CylinderBuilderFromDet::zmax

private

Definition at line 44 of file CylinderBuilderFromDet.h.

float CylinderBuilderFromDet::zmin

private

Definition at line 43 of file CylinderBuilderFromDet.h.