EllipticalTube.cc

Go to the documentation of this file.

#include "DetectorDescription/Core/src/EllipticalTube.h"

#include <ostream>

#include "CLHEP/Units/GlobalSystemOfUnits.h"
#include "CLHEP/Units/SystemOfUnits.h"
#include "DataFormats/GeometryVector/interface/Pi.h"

void DDI::EllipticalTube::stream(std::ostream & os) const
{
  os << " xSemiAxis[cm]=" << p_[0]/cm
     << " ySemiAxis[cm]=" << p_[1]/cm
     << " zHeight[cm]=" << p_[2]/cm;
}

double DDI::EllipticalTube::volume() const { 
  double volume(0.);
  // who cares major or minor axis? pi * a * b == pi * xhalf * yhalf
  // area of a slice.
  // we KNOW they are all cm... CMS uses cm
  double area(Geom::pi()*p_[0]*p_[1]);
  //volume is z * area.  I'm pretty sure :)
  volume = area * p_[2] * 2;
  return volume; 

}