VPolyhedra.cxx

Go to the documentation of this file.

// $Id$
 
// -----------------------------------------------------------------------
// The BaseVGM package of the Virtual Geometry Model
// Copyright (C) 2007, Ivana Hrivnacova
// All rights reserved.
//
// For the licensing terms see vgm/LICENSE.
// Contact: ivana@ipno.in2p3.fr
// -----------------------------------------------------------------------
 
//
// Class VPolyhedra
// ---------------
// The ABC for polyhedra solids.
//
// Author: Ivana Hrivnacova; IPN Orsay
 
#include "VGM/common/Math.h"
 
#include "BaseVGM/solids/VPolyhedra.h"
 
//_____________________________________________________________________________
std::ostream& operator<<(std::ostream& out, const VGM::IPolyhedra& polyhedra)
{
  const VGM::ISolid& polyhedraSolid = polyhedra;
  out << polyhedraSolid;
  return out;
}
 
//_____________________________________________________________________________
BaseVGM::VPolyhedra::VPolyhedra() : VGM::IPolyhedra()
{
}
 
//_____________________________________________________________________________
BaseVGM::VPolyhedra::~VPolyhedra()
{
  //
}
 
//_____________________________________________________________________________
double BaseVGM::VPolyhedra::ConvertRadiusFactor() const
{
  //
  double phiTotal = DeltaPhi() / 180. * M_PI;
  if ((phiTotal <= 0) || (phiTotal >= 2 * M_PI * (1 - 1e-16)))
    phiTotal = 2 * M_PI;
 
  return cos(0.5 * phiTotal / NofSides());
}
 
//_____________________________________________________________________________
std::ostream& BaseVGM::VPolyhedra::Put(std::ostream& out) const
{
  out << "  sphi = " << StartPhi() << "deg"
      << "  dphi = " << DeltaPhi() << "deg"
      << "  nsides = " << NofSides() << "  nz = " << NofZPlanes() << std::endl;
 
  double* zvalues = ZValues();
  double* rinvalues = InnerRadiusValues();
  double* routvalues = OuterRadiusValues();
 
  for (int i = 0; i < NofZPlanes(); i++) {
    out << "   " << i << "th plane: "
        << " z = " << zvalues[i] << "mm"
        << " rin = " << rinvalues[i] << "mm"
        << " rout = " << routvalues[i] << "mm";
 
    if (i < NofZPlanes() - 1) out << std::endl;
  }
 
  return out;
}