Torus.cc

Go to the documentation of this file.

#include "DetectorDescription/Core/src/Torus.h"
#include "DataFormats/Math/interface/GeantUnits.h"

#include <cmath>
#include <vector>

#include "DetectorDescription/Core/interface/DDSolidShapes.h"
#include "DetectorDescription/Core/src/Solid.h"

using std::sqrt;
using namespace geant_units;
using namespace geant_units::operators;

DDI::Torus::Torus( double pRMin,
           double pRMax,
           double pRTor,
           double pSPhi,
           double pDPhi
           )
 : Solid(DDSolidShape::ddtorus) 
{        
  p_.emplace_back(pRMin); // ......... 0
  p_.emplace_back(pRMax); // ......... 1
  p_.emplace_back(pRTor); // ......... 2
  p_.emplace_back(pSPhi); // ......... 3
  p_.emplace_back(pDPhi); // ......... 4
}


void DDI::Torus::stream(std::ostream & os) const
{
  os << " rMin=" << convertMmToCm( p_[0] )
     << " rMax=" << convertRadToDeg( p_[1] )
     << " rTor=" << convertRadToDeg( p_[2] )
     << " sPhi=" << convertMmToCm( p_[3] )
     << " dPhi=" << convertMmToCm( p_[4] );
}

double DDI::Torus::volume() const
{
  double volume=0;
  
  /* use notation as described in documentation about geant 4 shapes */
  
  // From Geant4: { fCubicVolume = fDPhi*pi*fRtor*(fRmax*fRmax-fRmin*fRmin);
  
  volume = p_[4] * piRadians * p_[2] * (p_[1] * p_[1] - p_[0] * p_[0]);
  
  return volume;
}