DDI::Polycone Class Reference

#include <Polycone.h>

Inheritance diagram for DDI::Polycone:

Public Member Functions
	Polycone (double startPhi, double deltaPhi, const std::vector< double > &z, const std::vector< double > &rmin, const std::vector< double > &rmax)
	Polycone (double startPhi, double deltaPhi, const std::vector< double > &z, const std::vector< double > &r)
void	stream (std::ostream &) const override
double	volume () const override
Public Member Functions inherited from DDI::Solid
const std::vector< double > &	parameters () const
void	setParameters (std::vector< double > const &p)
DDSolidShape	shape () const
	Solid ()
	Solid (DDSolidShape shape)
virtual	~Solid ()

Additional Inherited Members
Protected Attributes inherited from DDI::Solid
std::vector< double >	p_
DDSolidShape	shape_

Detailed Description

Definition at line 11 of file Polycone.h.

Constructor & Destructor Documentation

Polycone::Polycone	(	double	startPhi,
		double	deltaPhi,
		const std::vector< double > &	z,
		const std::vector< double > &	rmin,
		const std::vector< double > &	rmax
	)

Definition at line 15 of file Polycone.cc.

References ddpolycone_rrz, Exception, mps_fire::i, and DDI::Solid::p_.

    : Solid(DDSolidShape::ddpolycone_rrz) {
  p_.emplace_back(startPhi);
  p_.emplace_back(deltaPhi);
  if ((z.size() != rmin.size()) || (z.size() != rmax.size())) {
    throw cms::Exception("DDException") << "Polycone(..): std::vectors z,rmin,rmax not of same length";
  } else {
    for (unsigned int i = 0; i < z.size(); ++i) {
      p_.emplace_back(z[i]);
      p_.emplace_back(rmin[i]);
      p_.emplace_back(rmax[i]);
    }
  }
}

Polycone::Polycone	(	double	startPhi,
		double	deltaPhi,
		const std::vector< double > &	z,
		const std::vector< double > &	r
	)

Definition at line 34 of file Polycone.cc.

References ddpolycone_rz, Exception, mps_fire::i, and DDI::Solid::p_.

    : Solid(DDSolidShape::ddpolycone_rz) {
  p_.emplace_back(startPhi);
  p_.emplace_back(deltaPhi);
  if (z.size() != r.size()) {
    throw cms::Exception("DDException") << "Polycone(..): std::vectors z,rmin,rmax not of same length";
  } else {
    for (unsigned int i = 0; i < z.size(); ++i) {
      p_.emplace_back(z[i]);
      p_.emplace_back(r[i]);
    }
  }
}

Member Function Documentation

void DDI::Polycone::stream ( std::ostream &  os ) const

overridevirtual

Reimplemented from DDI::Solid.

Definition at line 108 of file Polycone.cc.

References angle_units::operators::convertMmToCm(), angle_units::operators::convertRadToDeg(), and isotrackApplyRegressor::k.

                                             {
  os << " startPhi[deg]=" << convertRadToDeg(p_[0]) << " dPhi[deg]=" << convertRadToDeg(p_[1]) << " Sizes[cm]=";
  for (unsigned k = 2; k < p_.size(); ++k)
    os << convertMmToCm(p_[k]) << " ";
}

double Polycone::volume ( void  ) const

overridevirtual

Reimplemented from DDI::Solid.

Definition at line 48 of file Polycone.cc.

References cms::cuda::assert(), ddpolycone_rrz, ddpolycone_rz, PVValHelper::dz, mps_fire::i, dqmiolumiharvest::j, heppy_loop::loop, DDI::Solid::p_, mps_fire::result, alignCSCRings::s, and DDI::Solid::shape_.

                              {
  double result = -1.;
  if (shape_ == DDSolidShape::ddpolycone_rrz) {
    unsigned int loop = (p_.size() - 2) / 3 - 1;
    assert(loop > 0);
    double sec = 0;
    int i = 2;
    for (unsigned int j = 2; j < (loop + 2); ++j) {
      double dz = std::fabs(p_[i] - p_[i + 3]);
      double v_min = dz * 1_pi / 3. * (p_[i + 1] * p_[i + 1] + p_[i + 4] * p_[i + 4] + p_[i + 1] * p_[i + 4]);
      double v_max = dz * 1_pi / 3. * (p_[i + 2] * p_[i + 2] + p_[i + 5] * p_[i + 5] + p_[i + 2] * p_[i + 5]);
      double s = v_max - v_min;
      //assert(s>=0);
      sec += s;
      i += 3;
    }
    result = sec * std::fabs(p_[1]) / 2._pi;
  }

  if (shape_ == DDSolidShape::ddpolycone_rz) {
    double volume = 0;
    double phiFrom = p_[0];
    double phiTo = p_[0] + p_[1];
    double slice = (std::fabs(phiFrom - phiTo)) / 2_pi;
    double zBegin = 0;
    double zEnd = 0;
    double rBegin = 0;
    double rEnd = 0;
    double z = 0;
    unsigned int i = 2;

    while (i < (p_.size() - 2)) {
      zBegin = p_[i];
      zEnd = p_[i + 2];
      rBegin = p_[i + 1];
      rEnd = p_[i + 3];
      z = zBegin - zEnd;

      /* for calculation of volume1 look at calculation of DDConsImpl of a volume1. Furthermore z can be smaller than zero. This makes sense since volumes we have to substract */
      double volume1 = (rEnd * rEnd + rBegin * rBegin + rBegin * rEnd) * z / 3;
      volume = volume + volume1;
      i = i + 2;
    }

    /* last line (goes from last z/r value to first */
    i = p_.size() - 2;
    zBegin = p_[i];
    zEnd = p_[2];
    rBegin = p_[i + 1];
    rEnd = p_[3];
    z = zBegin - zEnd;

    double volume2 = (rEnd * rEnd + rBegin * rBegin + rBegin * rEnd) * z / 3;
    volume = volume + volume2;
    volume = std::fabs(slice * 1_pi * volume);
    result = volume;
  }
  return result;
}