BoundSpan.cc

Go to the documentation of this file.

#include "DataFormats/GeometrySurface/interface/BoundSpan.h"
#include "DataFormats/GeometrySurface/interface/Surface.h"
#include "DataFormats/GeometrySurface/interface/TrapezoidalPlaneBounds.h"
#include "DataFormats/GeometrySurface/interface/RectangularPlaneBounds.h"

#include "DataFormats/GeometryVector/interface/VectorUtil.h"

void BoundSpan::compute(Surface const & plane) {
  const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
  const RectangularPlaneBounds* rectangularBounds( dynamic_cast<const RectangularPlaneBounds*>(&(plane.bounds())));  
  
  Surface::GlobalPoint corners[8];
  
  if (trapezoidalBounds) {
    std::array<const float, 4> const & parameters = (*trapezoidalBounds).parameters();
    
    float hbotedge = parameters[0];
    float htopedge = parameters[1];
    float hapothem = parameters[3];   
    float thickness =  (*trapezoidalBounds).thickness();
    
    corners[0] = plane.toGlobal( LocalPoint( -htopedge, hapothem,  thickness/2));
    corners[1] = plane.toGlobal( LocalPoint(  htopedge, hapothem,  thickness/2));
    corners[2] = plane.toGlobal( LocalPoint(  hbotedge, -hapothem, thickness/2));
    corners[3] = plane.toGlobal( LocalPoint( -hbotedge, -hapothem, thickness/2));
    corners[4] = plane.toGlobal( LocalPoint( -htopedge, hapothem,  -thickness/2));
    corners[5] = plane.toGlobal( LocalPoint(  htopedge, hapothem,  -thickness/2));
    corners[6] = plane.toGlobal( LocalPoint(  hbotedge, -hapothem, -thickness/2));
    corners[7] = plane.toGlobal( LocalPoint( -hbotedge, -hapothem, -thickness/2));
    
  }else if(rectangularBounds) {
    float length = rectangularBounds->length();
    float width  = rectangularBounds->width();   
    float thickness =  (*rectangularBounds).thickness();
      
    corners[0] = plane.toGlobal( LocalPoint( -width/2, -length/2, thickness/2));
    corners[1] = plane.toGlobal( LocalPoint( -width/2, +length/2, thickness/2));
    corners[2] = plane.toGlobal( LocalPoint( +width/2, -length/2, thickness/2));
    corners[3] = plane.toGlobal( LocalPoint( +width/2, +length/2, thickness/2));
    corners[4] = plane.toGlobal( LocalPoint( -width/2, -length/2, -thickness/2));
    corners[5] = plane.toGlobal( LocalPoint( -width/2, +length/2, -thickness/2));
    corners[6] = plane.toGlobal( LocalPoint( +width/2, -length/2, -thickness/2));
    corners[7] = plane.toGlobal( LocalPoint( +width/2, +length/2, -thickness/2));
  }else{  
    
  }
  
  float phimin = corners[0].barePhi(); float phimax = phimin;
  float zmin   = corners[0].z();    float zmax   = zmin;
  for ( int i = 1; i < 8; i++ ) {
    float cPhi = corners[i].barePhi();
    if ( Geom::phiLess( cPhi, phimin)) { phimin = cPhi; }
    if ( Geom::phiLess( phimax, cPhi)) { phimax = cPhi; }
    float z = corners[i].z();
    if ( z < zmin) zmin = z;
    if ( z > zmax) zmax = z;  
  }
  m_zSpan.first    = zmin;
  m_zSpan.second   = zmax;
  m_phiSpan.first  = phimin;
  m_phiSpan.second = phimax;
}