Tubs.cc

Go to the documentation of this file.

#include "DetectorDescription/Core/src/Tubs.h"
 
#include <cmath>
#include <vector>
 
#include "DetectorDescription/Core/interface/DDSolidShapes.h"
#include "DataFormats/Math/interface/GeantUnits.h"
#include "DetectorDescription/Core/src/Solid.h"
 
using std::fabs;
using namespace geant_units::operators;
 
DDI::Tubs::Tubs(double zhalf, double rIn, double rOut, double startPhi, double deltaPhi) : Solid(DDSolidShape::ddtubs) {
  p_.emplace_back(zhalf);
  p_.emplace_back(rIn);
  p_.emplace_back(rOut);
  p_.emplace_back(startPhi);
  p_.emplace_back(deltaPhi);
}
 
void DDI::Tubs::stream(std::ostream& os) const {
  os << " zhalf=" << convertMmToCm(p_[0]) << " rIn=" << convertMmToCm(p_[1]) << " rOut=" << convertMmToCm(p_[2])
     << " startPhi=" << convertRadToDeg(p_[3]) << " deltaPhi=" << convertRadToDeg(p_[4]);
}
 
double DDI::Tubs::volume() const {
  double volume = 0;
  double z = 2. * p_[0];
  double rIn = p_[1];
  double rOut = p_[2];
  double phi = p_[4];
 
  double volume1 = 1_pi * rIn * rIn * z;
  double volume2 = 1_pi * rOut * rOut * z;
 
  double slice = fabs(phi / (2_pi));
 
  volume = slice * (volume2 - volume1);
 
  return volume;
}