DDAngular Class Reference

Inheritance diagram for DDAngular:

Public Member Functions
	DDAngular ()
void	execute (DDCompactView &cpv) override
void	initialize (const DDNumericArguments &nArgs, const DDVectorArguments &vArgs, const DDMapArguments &mArgs, const DDStringArguments &sArgs, const DDStringVectorArguments &vsArgs) override

Private Member Functions
DD3Vector	fUnitVector (double theta, double phi)

Private Attributes
std::vector< double >	m_center
std::pair< std::string, std::string >	m_childNmNs
double	m_delta
std::string	m_idNameSpace
int	m_incrCopyNo
int	m_n
double	m_radius
double	m_rangeAngle
std::vector< double >	m_rotateSolid
DDRotationMatrix	m_solidRot
double	m_startAngle
int	m_startCopyNo

Detailed Description

Definition at line 21 of file DDAngular.cc.

Constructor & Destructor Documentation

DDAngular::DDAngular ( )

inline

Definition at line 24 of file DDAngular.cc.

References cmsIgProf_Analysis::execute(), initialize(), and theta().

    : m_n( 1 ),
      m_startCopyNo( 1 ),
      m_incrCopyNo( 1 ),
      m_startAngle( 0. ),
      m_rangeAngle( 360._deg ),
      m_radius( 0. ),
      m_delta( 0. )
  {}

Member Function Documentation

void DDAngular::execute ( DDCompactView &  cpv )

override

Definition at line 136 of file DDAngular.cc.

References geant_units::operators::convertRadToDeg(), popcon2dropbox::copy(), funct::cos(), DDcreateRotationMatrix(), DDrot(), mps_fire::i, LogDebug, class-composition::parent, DDCompactView::position(), idealTransformation::rotation, funct::sin(), AlCaHLTBitMon_QueryRunRegistry::string, and theta().

{
  DDName mother = parent().name();
  DDName ddname( m_childNmNs.first, m_childNmNs.second );
  double theta  = 90._deg;
  int    copy   = m_startCopyNo;
  double phi    = m_startAngle;

  for( int i = 0; i < m_n; ++i )
  {
    double phix = phi;
    double phiy = phix + 90._deg;
    double phideg = convertRadToDeg( phix );

    std::string rotstr = m_childNmNs.first + "_" + std::to_string( phideg * 10.);
    DDRotation rotation = DDRotation( DDName( rotstr, m_idNameSpace ));
    if( !rotation )
    {
      LogDebug( "DDAlgorithm" ) << "DDAngular: Creating a new "
                << "rotation: " << rotstr << "\t90., " 
                << convertRadToDeg( phix ) << ", 90.," 
                << convertRadToDeg( phiy ) << ", 0, 0";
    
      rotation = DDrot( DDName( rotstr, m_idNameSpace ),
            std::make_unique<DDRotationMatrix>( *DDcreateRotationMatrix( theta, phix, theta, phiy,
                                             0., 0. ) * m_solidRot )); 
    }
      
    double xpos = m_radius * cos( phi ) + m_center[0];
    double ypos = m_radius * sin( phi ) + m_center[1];
    double zpos = m_center[2];
    DDTranslation tran( xpos, ypos, zpos );
    
    cpv.position( ddname, mother, copy, tran, rotation );
    LogDebug( "DDAlgorithm" ) << "DDAngular: child " << m_childNmNs.second << ":" << m_childNmNs.first << " number " 
                  << copy << " positioned in " << mother << " at "
                  << tran << " with " << rotation << "\n";
    copy += m_incrCopyNo;
    phi  += m_delta;
  }
}

DD3Vector DDAngular::fUnitVector ( double  theta, double  phi  )

private

Definition at line 179 of file DDAngular.cc.

References funct::cos(), DEFINE_EDM_PLUGIN, and funct::sin().

{
  return DD3Vector( cos( phi ) * sin( theta ),
            sin( phi ) * sin( theta ),
            cos( theta ));
}

void DDAngular::initialize ( const DDNumericArguments &  nArgs, const DDVectorArguments &  vArgs, const DDMapArguments &  mArgs, const DDStringArguments &  sArgs, const DDStringVectorArguments &  vsArgs  )

override

Definition at line 62 of file DDAngular.cc.

References funct::abs(), geant_units::operators::convertRadToDeg(), DDSplit(), mps_fire::i, createfilelist::int, LogDebug, DDCurrentNamespace::ns(), class-composition::parent, and groupFilesInBlocks::temp.

{
  m_n           = int(nArgs["N"]);
  m_startCopyNo = int(nArgs["StartCopyNo"]);
  m_incrCopyNo  = int(nArgs["IncrCopyNo"]);
  m_rangeAngle  = nArgs["RangeAngle"];
  m_startAngle  = nArgs["StartAngle"];
  m_radius      = nArgs["Radius"];
  m_center      = vArgs["Center"];
  m_rotateSolid = vArgs["RotateSolid"];

  m_solidRot   = DDRotationMatrix();
  
  if( std::abs( m_rangeAngle - 360.0_deg ) < 0.001_deg )
  { 
    m_delta = m_rangeAngle / double( m_n );
  }
  else
  {
    if( m_n > 1 )
    {
      m_delta = m_rangeAngle / double( m_n - 1 );
    }
    else
    {
      m_delta = 0.;
    }
  }  
  
  LogDebug( "DDAlgorithm" ) << "DDAngular: Parameters for position"
                << "ing:: n " << m_n << " Start, Range, Delta " 
                << convertRadToDeg( m_startAngle ) << " " 
                << convertRadToDeg( m_rangeAngle ) << " " << convertRadToDeg( m_delta )
                << " Radius " << m_radius << " Centre " << m_center[0] 
                << ", " << m_center[1] << ", " << m_center[2];
  
  //======= collect data concerning the rotation of the solid 
  using sz_type = std::map<std::string, std::vector<double> >::mapped_type::size_type;
  sz_type sz = m_rotateSolid.size();
  if( sz%3 )
  {
    LogDebug( "DDAlgorithm" ) << "\trotateSolid must occur 3*n times (defining n subsequent rotations)\n"
                  << "\t  currently it appears " << sz << " times!\n";
  }
  for( sz_type i = 0; i < sz; i += 3 )
  {
    if(( m_rotateSolid[i] > 180._deg ) || ( m_rotateSolid[i] < 0._deg ))
    {
      LogDebug( "DDAlgorithm" ) << "\trotateSolid \'theta\' must be in range [0,180*deg]\n"
                << "\t  currently it is " << convertRadToDeg( m_rotateSolid[i] ) 
                << "*deg in rotateSolid[" << double(i) << "]!\n";
    }
    DDAxisAngle temp( fUnitVector( m_rotateSolid[i], m_rotateSolid[i + 1] ),
              m_rotateSolid[i + 2] );
    LogDebug( "DDAlgorithm" ) << "  rotsolid[" << i <<  "] axis=" << temp.Axis() << " rot.angle=" << convertRadToDeg( temp.Angle() );
    m_solidRot = temp * m_solidRot;           
  }

  m_idNameSpace = DDCurrentNamespace::ns();
  m_childNmNs   = DDSplit( sArgs["ChildName"] );
  if( m_childNmNs.second.empty())
    m_childNmNs.second = DDCurrentNamespace::ns();

  DDName parentName = parent().name();
  LogDebug( "DDAlgorithm" ) << "DDAngular: Parent " << parentName 
                << "\tChild " << m_childNmNs.first << "\tNameSpace "
                << m_childNmNs.second;
}

Member Data Documentation

std::vector<double> DDAngular::m_center

private

Definition at line 51 of file DDAngular.cc.

std::pair<std::string, std::string> DDAngular::m_childNmNs

private

Definition at line 55 of file DDAngular.cc.

double DDAngular::m_delta

private

Definition at line 50 of file DDAngular.cc.

std::string DDAngular::m_idNameSpace

private

Definition at line 54 of file DDAngular.cc.

int DDAngular::m_incrCopyNo

private

Definition at line 46 of file DDAngular.cc.

int DDAngular::m_n

private

Definition at line 44 of file DDAngular.cc.

double DDAngular::m_radius

private

Definition at line 49 of file DDAngular.cc.

double DDAngular::m_rangeAngle

private

Definition at line 48 of file DDAngular.cc.

std::vector<double> DDAngular::m_rotateSolid

private

Definition at line 52 of file DDAngular.cc.

DDRotationMatrix DDAngular::m_solidRot

private

Definition at line 58 of file DDAngular.cc.

double DDAngular::m_startAngle

private

Definition at line 47 of file DDAngular.cc.

int DDAngular::m_startCopyNo

private

Definition at line 45 of file DDAngular.cc.