DDAngular.cc File Reference

#include "DD4hep/DetFactoryHelper.h"
#include "DataFormats/Math/interface/GeantUnits.h"
#include "DetectorDescription/DDCMS/interface/DDPlugins.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <Math/Rotation3D.h>
#include <Math/AxisAngle.h>
#include <Math/DisplacementVector3D.h>

Go to the source code of this file.

Typedefs
using	DD3Vector = ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > >
using	DDAxisAngle = ROOT::Math::AxisAngle

Functions
static long	algorithm (Detector &, cms::DDParsingContext &ctxt, xml_h e, SensitiveDetector &)
	DD4HEP_OPEN_PLUGIN (dd4hep, ddcms_det_element_DDCMS_global_DDAngular)

Typedef Documentation

using DD3Vector = ROOT::Math::DisplacementVector3D<ROOT::Math::Cartesian3D<double> >

Definition at line 15 of file DDAngular.cc.

using DDAxisAngle = ROOT::Math::AxisAngle

Definition at line 16 of file DDAngular.cc.

Function Documentation

static long algorithm ( Detector &  , cms::DDParsingContext &  ctxt, xml_h  e, SensitiveDetector &    )

static

Definition at line 28 of file DDAngular.cc.

References writedatasetfile::args, geant_units::operators::convertRadToDeg(), popcon2dropbox::copy(), funct::cos(), delta, cms::DDAlgoArguments::find(), mps_fire::i, LogDebug, cms::makeRotation3D(), gen::n, cms::DDNamespace::name(), NAMESPACE_SEP, cms::DDNamespace::nsName(), cms::DDAlgoArguments::parentName(), cms::DDNamespace::prepend(), TCMET_cfi::radius, idealTransformation::rotation, cms::DDNamespace::rotation(), cms::DDParsingContext::rotations, funct::sin(), groupFilesInBlocks::temp, theta(), cms::DDAlgoArguments::value(), and cms::DDNamespace::volume().

{
  cms::DDNamespace ns( ctxt, e, true );
  DDAlgoArguments  args( ctxt, e );

  int            n           = args.value<int>("N");
  int            startCopyNo = args.find("StartCopyNo") ? args.value<int>("StartCopyNo") : 1;
  int            incrCopyNo  = args.find("IncrCopyNo") ? args.value<int>("IncrCopyNo") : 1;
  double         rangeAngle  = args.value<double>("RangeAngle");
  double         startAngle  = args.value<double>("StartAngle");
  double         radius      = args.value<double>("Radius");
  vector<double> center      = args.value<vector<double> >("Center");
  vector<double> rotateSolid = args.value<vector<double> >("RotateSolid");
  Volume         mother      = ns.volume(args.parentName());
  string         childName   = args.value<string>("ChildName");
  if( strchr( childName.c_str(), NAMESPACE_SEP ) == nullptr )
    childName = ns.name() + childName;
  Volume         child       = ns.volume( childName );

  double         delta = 0e0;
  
  if( fabs( rangeAngle - 360.0_deg ) < 0.001_deg ) { 
    delta = rangeAngle/double( n );
  } else if( n > 1 ) {
    delta = rangeAngle/double( n - 1 );
  }

  LogDebug("DDAlgorithm") << "debug: Parameters for positioning:: n "
              << n << " Start, Range, Delta " 
              << convertRadToDeg( startAngle ) << " " 
              << convertRadToDeg( rangeAngle ) << " " << convertRadToDeg( delta )
              << " Radius " << radius << " Centre " << center[0]
              << ", " << center[1] << ", " << center[2]
              << ", Rotate solid " << rotateSolid[0] << ", " << rotateSolid[1]
              << ", " << rotateSolid[2];
  LogDebug("DDAlgorithm") << "debug: Parent " << mother.name() 
              << "\tChild " << child.name() << " NameSpace "
              << ns.name();

  Rotation3D solidRot = Rotation3D();  // Identity rotation
  auto sz = 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( unsigned int i = 0; i < sz; i += 3 )
  {
    if(( i > 180._deg ) || ( i < 0._deg ))
    {
      LogDebug( "DDAlgorithm" ) << "\trotateSolid \'theta\' must be in range [0,180*deg]\n"
                << "\t  currently it is " << convertRadToDeg( i )
                << "*deg in rotateSolid[" << double(i) << "]!\n";
    }
    DDAxisAngle temp( fUnitVector( rotateSolid[i], rotateSolid[i + 1] ),
              rotateSolid[i + 2] );
    LogDebug( "DDAlgorithm" ) << "  rotsolid[" << i <<  "] axis=" << temp.Axis() << " rot.angle=" << convertRadToDeg( temp.Angle() );
    solidRot = temp * solidRot;
  }

  double theta  = 90._deg;
  int    copy   = startCopyNo;
  double phi    = startAngle;
  for( int i = 0; i < n; ++i )
  {
    double phix = phi;
    double phiy = phix + 90._deg;
    double phideg = convertRadToDeg( phix );

    Rotation3D rotation = makeRotation3D( theta, phix, theta, phiy, 0., 0. ) * solidRot;
    string rotstr = ns.nsName( child.name()) + std::to_string( phideg * 10. );
    auto irot = ctxt.rotations.find( ns.prepend( rotstr ));
    if( irot != ctxt.rotations.end()) {
      rotation = ns.rotation( ns.prepend( rotstr ));
    }
    
    double xpos = radius * cos( phi ) + center[0];
    double ypos = radius * sin( phi ) + center[1];
    double zpos = center[2];
    Position tran( xpos, ypos, zpos );
    mother.placeVolume( child, copy, Transform3D( rotation, tran ));
    LogDebug("DDAlgorithm") << "test " << child.name() << " number " 
                << copy << " positioned in " << mother.name() << " at "
                << tran  << " with " << rotstr << ": " << rotation;
    copy += incrCopyNo;
    phi  += delta;
  }
  return 1;
}

DD4HEP_OPEN_PLUGIN	(	dd4hep	,
		ddcms_det_element_DDCMS_global_DDAngular
	)

Definition at line 122 of file DDAngular.cc.