DDG4SolidConverter Class Reference

#include <DDG4SolidConverter.h>

Public Types
typedef G4VSolid *(*	FNPTR) (const DDSolid &)

Public Member Functions
G4VSolid *	convert (const DDSolid &)
	DDG4SolidConverter ()
	~DDG4SolidConverter ()

Static Private Member Functions
static G4VSolid *	box (const DDSolid &)
static G4VSolid *	cons (const DDSolid &)
static G4VSolid *	cuttubs (const DDSolid &)
static G4VSolid *	ellipsoid (const DDSolid &)
static G4VSolid *	ellipticaltube (const DDSolid &)
static G4VSolid *	extrudedpolygon (const DDSolid &)
static G4VSolid *	intersection (const DDSolid &)
static G4VSolid *	multiunion (const DDSolid &)
static G4VSolid *	orb (const DDSolid &)
static G4VSolid *	para (const DDSolid &)
static G4VSolid *	polycone_rrz (const DDSolid &)
static G4VSolid *	polycone_rz (const DDSolid &)
static G4VSolid *	polyhedra_rrz (const DDSolid &)
static G4VSolid *	polyhedra_rz (const DDSolid &)
static G4VSolid *	pseudotrap (const DDSolid &s)
static G4VSolid *	reflected (const DDSolid &)
static G4VSolid *	shapeless (const DDSolid &)
static G4VSolid *	sphere (const DDSolid &)
static G4VSolid *	subtraction (const DDSolid &)
static G4VSolid *	torus (const DDSolid &)
static G4VSolid *	trap (const DDSolid &)
static G4VSolid *	trunctubs (const DDSolid &)
static G4VSolid *	tubs (const DDSolid &)
static G4VSolid *	unionsolid (const DDSolid &)

Private Attributes
std::map< DDSolidShape, FNPTR >	convDispatch_

Static Private Attributes
static const std::vector< double > *	par_ = 0
static G4RotationMatrix *	rot = 0

Friends
class	testPseudoTrap
class	testTruncTubs

Detailed Description

Definition at line 14 of file DDG4SolidConverter.h.

Member Typedef Documentation

typedef G4VSolid*(* DDG4SolidConverter::FNPTR) (const DDSolid &)

Definition at line 19 of file DDG4SolidConverter.h.

Constructor & Destructor Documentation

DDG4SolidConverter::DDG4SolidConverter

(

)

Definition at line 13 of file DDG4SolidConverter.cc.

References box(), cons(), convDispatch_, cuttubs(), ddbox, ddcons, ddcuttubs, ddellipsoid, ddellipticaltube, ddextrudedpolygon, ddintersection, ddmultiunion, ddorb, ddparallelepiped, ddpolycone_rrz, ddpolycone_rz, ddpolyhedra_rrz, ddpolyhedra_rz, ddpseudotrap, ddreflected, ddsphere, ddsubtraction, ddtorus, ddtrap, ddtrunctubs, ddtubs, ddunion, ellipsoid(), ellipticaltube(), extrudedpolygon(), intersection(), multiunion(), orb(), para(), polycone_rrz(), polycone_rz(), polyhedra_rrz(), polyhedra_rz(), pseudotrap(), reflected(), sphere(), subtraction(), torus(), trap(), trunctubs(), tubs(), and unionsolid().

Referenced by intersection(), reflected(), subtraction(), and unionsolid().

                                       {
  // could also be done 'dynamically' from outside 
  // would then need to have a 'register' method ...
  convDispatch_[ddbox]            = DDG4SolidConverter::box; 
  convDispatch_[ddtubs]           = DDG4SolidConverter::tubs;
  convDispatch_[ddcuttubs]        = DDG4SolidConverter::cuttubs;
  convDispatch_[ddtrap]           = DDG4SolidConverter::trap;
  convDispatch_[ddcons]           = DDG4SolidConverter::cons;
  convDispatch_[ddpolycone_rrz]   = DDG4SolidConverter::polycone_rrz;
  convDispatch_[ddpolycone_rz]    = DDG4SolidConverter::polycone_rz;
  convDispatch_[ddpolyhedra_rrz]  = DDG4SolidConverter::polyhedra_rrz;
  convDispatch_[ddpolyhedra_rz]   = DDG4SolidConverter::polyhedra_rz;
  convDispatch_[ddextrudedpolygon]= DDG4SolidConverter::extrudedpolygon;
  convDispatch_[ddtorus]          = DDG4SolidConverter::torus;   
  convDispatch_[ddreflected]      = DDG4SolidConverter::reflected;
  convDispatch_[ddunion]          = DDG4SolidConverter::unionsolid;
  convDispatch_[ddintersection]   = DDG4SolidConverter::intersection;
  convDispatch_[ddsubtraction]    = DDG4SolidConverter::subtraction;
  convDispatch_[ddmultiunion]     = DDG4SolidConverter::multiunion;
  convDispatch_[ddpseudotrap]     = DDG4SolidConverter::pseudotrap;
  convDispatch_[ddtrunctubs]      = DDG4SolidConverter::trunctubs;
  convDispatch_[ddsphere]         = DDG4SolidConverter::sphere;   
  convDispatch_[ddorb]            = DDG4SolidConverter::orb;   
  convDispatch_[ddellipticaltube] = DDG4SolidConverter::ellipticaltube;   
  convDispatch_[ddellipsoid]      = DDG4SolidConverter::ellipsoid;   
  convDispatch_[ddparallelepiped] = DDG4SolidConverter::para;   
}

DDG4SolidConverter::~DDG4SolidConverter

(

)

Definition at line 41 of file DDG4SolidConverter.cc.

{}

Member Function Documentation

G4VSolid * DDG4SolidConverter::box ( const DDSolid &  solid )

staticprivate

Definition at line 64 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), and DDBase< N, C >::name().

Referenced by DDG4SolidConverter(), and trunctubs().

                                                        {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: box = " << solid;
  return new G4Box(solid.name().name(), (*par_)[0],(*par_)[1],(*par_)[2]);
}

G4VSolid * DDG4SolidConverter::cons ( const DDSolid &  solid )

staticprivate

Definition at line 111 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), and DDBase< N, C >::name().

Referenced by DDG4SolidConverter().

                                                         {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: cons = " << solid;  
  return new G4Cons(solid.name().name(), (*par_)[1],  // rmin -z
            (*par_)[2],  // rmax -z
            (*par_)[3],  // rmin +z
            (*par_)[4],  // rmax +z
            (*par_)[0],  // zHalf
            (*par_)[5],  // phistart
            (*par_)[6]); // deltaphi
}

G4VSolid * DDG4SolidConverter::convert

(

const DDSolid &

solid

)

Definition at line 43 of file DDG4SolidConverter.cc.

References convDispatch_, Exception, par_, DDSolid::parameters(), mps_fire::result, DDSolid::shape(), and DDBase< N, C >::toString().

Referenced by DDG4Builder::convertSolid().

                                                            {
  if ( !solid ) {
    edm::LogError("SimG4CoreGeometry") <<" DDG4SolidConverter::convert(..) found an undefined DDSolid " << solid.toString();
    throw cms::Exception("SimG4CoreGeometry", "DDG4SolidConverter::convert(..) found an undefined DDSolid " + solid.toString());
  }
  G4VSolid * result = nullptr;
  par_ = &(solid.parameters());
  std::map<DDSolidShape,FNPTR>::iterator it = convDispatch_.find(solid.shape());
  if (it != convDispatch_.end()) {
    result = it->second(solid);
  } else {
    throw cms::Exception("DetectorDescriptionFault") 
      <<  "DDG4SolidConverter::convert: conversion failed for s=" << solid
      << "\n solid.shape()=" << solid.shape()
      << std::endl;
  }
  return result;
}

G4VSolid * DDG4SolidConverter::cuttubs ( const DDSolid &  solid )

staticprivate

Definition at line 81 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), and DDBase< N, C >::name().

Referenced by DDG4SolidConverter().

                                                            {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: tubs = " << solid;  
  return new G4CutTubs(solid.name().name(), (*par_)[1], // rmin
               (*par_)[2], // rmax
               (*par_)[0], // dzHalf
               (*par_)[3], // phiStart
               (*par_)[4], // deltaPhi
               G4ThreeVector((*par_)[5],(*par_)[6],(*par_)[7]),
               G4ThreeVector((*par_)[8],(*par_)[9],(*par_)[10]));
}

G4VSolid * DDG4SolidConverter::ellipsoid ( const DDSolid &  solid )

staticprivate

Definition at line 598 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), DDBase< N, C >::name(), DDEllipsoid::xSemiAxis(), DDEllipsoid::ySemiAxis(), DDEllipsoid::zBottomCut(), DDEllipsoid::zSemiAxis(), and DDEllipsoid::zTopCut().

Referenced by DDG4SolidConverter().

                                                              {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: ellipsoid = " << solid;
  DDEllipsoid sp(solid);
  return new G4Ellipsoid(solid.name().name(),
             sp.xSemiAxis(),
             sp.ySemiAxis(),
             sp.zSemiAxis(),
             sp.zBottomCut(),
             sp.zTopCut());
}

G4VSolid * DDG4SolidConverter::ellipticaltube ( const DDSolid &  solid )

staticprivate

Definition at line 588 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), DDBase< N, C >::name(), DDEllipticalTube::xSemiAxis(), DDEllipticalTube::ySemiAxis(), and DDEllipticalTube::zHeight().

Referenced by DDG4SolidConverter().

                                                                   {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: ellipticaltube = " << solid;
  DDEllipticalTube sp(solid);
  return new G4EllipticalTube(solid.name().name(),
                  sp.xSemiAxis(),
                  sp.ySemiAxis(),
                  sp.zHeight());
}

G4VSolid * DDG4SolidConverter::extrudedpolygon ( const DDSolid &  solid )

staticprivate

Definition at line 236 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), DDBase< N, C >::name(), x, y, and z.

Referenced by DDG4SolidConverter().

                                                                    {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: extr_pgon = " << solid;
  std::vector<double> x = static_cast<DDExtrudedPolygon>(solid).xVec();
  std::vector<double> y = static_cast<DDExtrudedPolygon>(solid).yVec();
  std::vector<double> z = static_cast<DDExtrudedPolygon>(solid).zVec();
  std::vector<double> zx = static_cast<DDExtrudedPolygon>(solid).zxVec();
  std::vector<double> zy = static_cast<DDExtrudedPolygon>(solid).zyVec();
  std::vector<double> zs = static_cast<DDExtrudedPolygon>(solid).zscaleVec();

  std::vector<G4TwoVector> polygon;
  std::vector<G4ExtrudedSolid::ZSection> zsections;
  for( unsigned int it = 0; it < x.size(); ++it )
    polygon.emplace_back( x[it], y[it] );
  for( unsigned int it = 0; it < z.size(); ++it )
    zsections.emplace_back( z[it], G4TwoVector(zx[it], zy[it]), zs[it] );
  return new G4ExtrudedSolid( solid.name().name(), polygon, zsections );
}

G4VSolid * DDG4SolidConverter::intersection ( const DDSolid &  solid )

staticprivate

Definition at line 341 of file DDG4SolidConverter.cc.

References DDG4SolidConverter(), mps_fire::i, LogDebug, AnalysisDataFormats_SUSYBSMObjects::ms, multiunion(), DDName::name(), DDBase< N, C >::name(), DDRotation::rotation(), DDBooleanSolid::rotation(), DDBooleanSolid::solidA(), DDBooleanSolid::solidB(), DDMultiUnionSolid::solids(), and DDBooleanSolid::translation().

Referenced by DDG4SolidConverter().

                                                                 {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: intersection = " << solid;
  G4IntersectionSolid * us = nullptr;
  DDBooleanSolid bs(solid);
  if (bs) {
    G4VSolid * sa = DDG4SolidConverter().convert(bs.solidA());
    G4VSolid * sb = DDG4SolidConverter().convert(bs.solidB());
    LogDebug("SimG4CoreGeometry") << " name:" << solid.name() << " t=" << bs.translation() << std::flush;
    LogDebug("SimG4CoreGeometry") << " " << bs.rotation().rotation()->Inverse() << std::flush;
    std::vector<double> tdbl(9);
    bs.rotation().rotation()->Inverse().GetComponents(tdbl.begin(), tdbl.end());
    CLHEP::HepRep3x3 temprep(tdbl[0], tdbl[1], tdbl[2], tdbl[3], tdbl[4], tdbl[5], tdbl[6], tdbl[7], tdbl[8]);
    CLHEP::Hep3Vector temphvec(bs.translation().X(), bs.translation().Y(), bs.translation().Z()); 
    us = new G4IntersectionSolid(solid.name().name(),
                 sa,
                 sb,
                 new CLHEP::HepRotation(temprep),
                 temphvec);
  }
  return us;       
}

G4VSolid * DDG4SolidConverter::multiunion ( const DDSolid &  solid )

staticprivate

Definition at line 384 of file DDG4SolidConverter.cc.

Referenced by DDG4SolidConverter(), and intersection().

                                                                 {
  return nullptr;
}

G4VSolid * DDG4SolidConverter::orb ( const DDSolid &  solid )

staticprivate

Definition at line 581 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), DDBase< N, C >::name(), and DDOrb::radius().

Referenced by DDG4SolidConverter().

                                                        {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: orb = " << solid;  
  DDOrb sp(solid);
  return new G4Orb(solid.name().name(), sp.radius());
}

G4VSolid * DDG4SolidConverter::para ( const DDSolid &  solid )

staticprivate

Definition at line 610 of file DDG4SolidConverter.cc.

References DDParallelepiped::alpha(), LogDebug, DDName::name(), DDBase< N, C >::name(), DDParallelepiped::phi(), DDParallelepiped::theta(), DDParallelepiped::xHalf(), DDParallelepiped::yHalf(), and DDParallelepiped::zHalf().

Referenced by DDG4SolidConverter().

                                                         {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: parallelepiped = " << solid;  
  DDParallelepiped sp(solid);
  return new G4Para(solid.name().name(),
            sp.xHalf(),
            sp.yHalf(),
            sp.zHalf(),
            sp.alpha(),
            sp.theta(),
            sp.phi());
}

G4VSolid * DDG4SolidConverter::polycone_rrz ( const DDSolid &  solid )

staticprivate

Definition at line 154 of file DDG4SolidConverter.cc.

References KineDebug3::count(), mps_fire::i, LogDebug, DDName::name(), DDBase< N, C >::name(), and par_.

Referenced by DDG4SolidConverter().

                                                                 {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: pcon_rrz = " << solid;
  std::vector<double> z_p;
  std::vector<double> rmin_p;
  std::vector<double> rmax_p;
  std::vector<double>::const_iterator i = par_->begin()+2;
  int count = 0;
  for (; i!=par_->end(); ++i) {
    LogDebug("SimG4CoreGeometry") << "z=" << *i/CLHEP::cm;
    z_p.push_back(*i); ++i;
    LogDebug("SimG4CoreGeometry") << "rmin=" << *i/CLHEP::cm;
    rmin_p.push_back(*i); ++i;
    LogDebug("SimG4CoreGeometry") << "rmax=" << *i/CLHEP::cm;
    rmax_p.push_back(*i); 
    count++;
  }
  LogDebug("SimG4CoreGeometry") << "sp=" << (*par_)[0]/CLHEP::deg << " ep=" 
                << (*par_)[1]/CLHEP::deg;
  /*
    std::cout << "### Polycone_RRZ: " << "sp=" << (*par_)[0]/CLHEP::deg 
          << " ep=" << (*par_)[1]/CLHEP::deg 
          << " N= " << count << std::endl;
    for(int i=0; i<count; ++i) { 
      std::cout << " R1= " << rmin_p[i] << " R1= " << rmax_p[i] << " Z= " << z_p[i] << std::endl;
    }
  */
  return new G4Polycone(solid.name().name(), (*par_)[0], (*par_)[1], // start,delta-phi
                        count, // sections
            &(z_p[0]),
            &(rmin_p[0]),
            &(rmax_p[0]));
  
}

G4VSolid * DDG4SolidConverter::polycone_rz ( const DDSolid &  solid )

staticprivate

Definition at line 124 of file DDG4SolidConverter.cc.

References KineDebug3::count(), mps_fire::i, LogDebug, DDName::name(), DDBase< N, C >::name(), alignCSCRings::r, and z.

Referenced by DDG4SolidConverter().

                                                                {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: pcon_rz = " << solid;  
  std::vector<double> r;
  std::vector<double> z;
  std::vector<double>::const_iterator i = (*par_).begin()+2;
  int count=0;
  for(; i!=(*par_).end(); ++i) {
    LogDebug("SimG4CoreGeometry") << " z=" << *i/CLHEP::cm;
    z.push_back(*i); ++i;
    LogDebug("SimG4CoreGeometry") << " r=" << *i/CLHEP::cm;
    r.push_back(*i);
    count++;
  }
  LogDebug("SimG4CoreGeometry") << "sp=" << (*par_)[0]/CLHEP::deg << " ep=" 
                << (*par_)[1]/CLHEP::deg;
  /*
   std::cout << "### Polycone_RZ: " << "sp=" << (*par_)[0]/CLHEP::deg 
   << " ep=" << (*par_)[1]/CLHEP::deg 
   << " N= " << count << std::endl;
   for(int i=0; i<count; ++i) { 
   std::cout << " R= " << r[i] << " Z= " << z[i] << std::endl;
   }
  */
  return new G4Polycone(solid.name().name(), (*par_)[0], (*par_)[1], // start,delta-phi
            count, // numRZ
            &(r[0]),
            &(z[0]));
}

G4VSolid * DDG4SolidConverter::polyhedra_rrz ( const DDSolid &  solid )

staticprivate

Definition at line 210 of file DDG4SolidConverter.cc.

References KineDebug3::count(), mps_fire::i, createfilelist::int, LogDebug, DDName::name(), DDBase< N, C >::name(), and par_.

Referenced by DDG4SolidConverter().

                                                                  {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: phed_rrz = " << solid;
  std::vector<double> z_p;
  std::vector<double> rmin_p;
  std::vector<double> rmax_p;
  std::vector<double>::const_iterator i = par_->begin()+3;
  int count = 0;
  for (; i!=par_->end(); ++i) {
    LogDebug("SimG4CoreGeometry") << "z=" << *i/CLHEP::cm;
    z_p.push_back(*i); ++i;
    LogDebug("SimG4CoreGeometry") << "rmin=" << *i/CLHEP::cm;
    rmin_p.push_back(*i); ++i;
    LogDebug("SimG4CoreGeometry") << "rmax=" << *i/CLHEP::cm;
    rmax_p.push_back(*i); 
    count++;
  }
  LogDebug("SimG4CoreGeometry") << "sp=" << (*par_)[0]/CLHEP::deg << " ep=" 
                << (*par_)[1]/CLHEP::deg;
  return new G4Polyhedra(solid.name().name(), (*par_)[1], (*par_)[2], int((*par_)[0]), // start,delta-phi,sides
             count, // sections
             &(z_p[0]),
             &(rmin_p[0]),
             &(rmax_p[0]));  
}

G4VSolid * DDG4SolidConverter::polyhedra_rz ( const DDSolid &  solid )

staticprivate

Definition at line 190 of file DDG4SolidConverter.cc.

References KineDebug3::count(), mps_fire::i, createfilelist::int, LogDebug, DDName::name(), DDBase< N, C >::name(), par_, alignCSCRings::r, and z.

Referenced by DDG4SolidConverter().

                                                                 {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: phed_rz = " << solid;  
  std::vector<double> r;
  std::vector<double> z;
  std::vector<double>::const_iterator i = par_->begin()+3;
  int count=0;
    
  for(; i!=par_->end(); ++i) {
    z.push_back(*i); ++i;
    r.push_back(*i);
    count++;
  }
      
  return new G4Polyhedra(solid.name().name(), (*par_)[1], (*par_)[2], int((*par_)[0]),// start,delta-phi;sides
             count, // numRZ
             &(r[0]),
             &(z[0]));
}

G4VSolid * DDG4SolidConverter::pseudotrap ( const DDSolid &  s )

staticprivate

correct implementation, but fails to visualize due to G4/Iguana limitations

Definition at line 391 of file DDG4SolidConverter.cc.

References funct::abs(), DDPseudoTrap::atMinusZ(), delta, Exception, h, DDPseudoTrap::halfZ(), LogDebug, DDName::name(), dataset::name, DDBase< N, C >::name(), EnergyCorrector::pt, alignCSCRings::r, DDPseudoTrap::radius(), mps_fire::result, rot, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, trap(), tubs(), x, DDPseudoTrap::x1(), DDPseudoTrap::x2(), DDPseudoTrap::y1(), and DDPseudoTrap::y2().

Referenced by DDG4SolidConverter().

                                                               {
  if(nullptr == rot) {
    rot = new G4RotationMatrix;
    rot->rotateX(90.*deg);
  }    

  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: pseudoTrap = " << solid;
  G4Trd * trap = nullptr;
  G4Tubs * tubs = nullptr;
  G4VSolid * result = nullptr;
  DDPseudoTrap pt(solid); // pt...PseudoTrap
  double r = pt.radius();
  bool atMinusZ = pt.atMinusZ();
  double x = 0;
  double h = 0;
  bool intersec = false; // union or intersection solid
  if (pt.atMinusZ()) {
    x = pt.x1(); // tubs radius
  } 
  else {
    x = pt.x2(); // tubs radius
  }
  double openingAngle = 2.*asin(x/std::abs(r));
  //trap = new G4Trd(solid.name().name(), 
  double displacement=0;
  double startPhi=0;
  /* calculate the displacement of the tubs w.r.t. to the trap,
     determine the opening angle of the tubs */
  double delta = sqrt(r*r-x*x);
  if (r < 0 && std::abs(r) >= x) {
    intersec = true; // intersection solid
    h = pt.y1() < pt.y2() ? pt.y2() : pt.y1(); // tubs half height
    h += h/20.; // enlarge a bit - for subtraction solid
    if (atMinusZ) {
      displacement = - pt.halfZ() - delta; 
      startPhi = 270.*deg - openingAngle/2.;
    }
    else {
      displacement =   pt.halfZ() + delta;
      startPhi = 90.*deg - openingAngle/2.;
    }
  }
  else if ( r > 0 && std::abs(r) >= x )
    {
      if (atMinusZ) {
        displacement = - pt.halfZ() + delta;
    startPhi = 90.*deg - openingAngle/2.;
    h = pt.y1();
      }
      else {
        displacement =   pt.halfZ() - delta; 
    startPhi = 270.*deg - openingAngle/2.;
    h = pt.y2();
      }    
    }
  else {
    throw cms::Exception("DetectorDescriptionFault", "Check parameters of the PseudoTrap! name=" + pt.name().name());
   }
  G4ThreeVector displ(0.,0.,displacement); // displacement of the tubs w.r.t. trap
  LogDebug("SimG4CoreGeometry") << "DDSolidConverter::pseudotrap(): displacement=" << displacement 
                << " openingAngle=" << openingAngle/deg << " x=" << x << " h=" << h;
    
  // Now create two solids (trd & tubs), and a boolean solid out of them 
  std::string name=pt.name().name();
  trap = new G4Trd(name, pt.x1(), pt.x2(), pt.y1(), pt.y2(), pt.halfZ());
  tubs = new G4Tubs(name, 
            0., // rMin
            std::abs(r), // rMax
            h, // half height
            startPhi, // start angle
            openingAngle);
  if (intersec) {
    result = new G4SubtractionSolid(name, trap, tubs, rot, displ);
  }
  else {
    G4VSolid * tubicCap = new G4SubtractionSolid(name, 
                         tubs, 
                         new G4Box(name, 1.1*x, sqrt(r*r-x*x), 1.1*h),  
                         nullptr, 
                         G4ThreeVector());
    result = new G4UnionSolid(name, trap, tubicCap, rot, displ);
            
    // approximative implementation - also fails to visualize due to G4/Iguana limitations
    /*
      delete tubs;
      tubs = new G4Tubs(name, 
      sqrt(r*r-x*x), // rMin-approximation!
      std::abs(r), // rMax
      h, // half height
      startPhi, // start angle
      openingAngle);
      result = new G4UnionSolid(name, trap, tubs, rot, displ);
    */
  }                
  return result;
}

G4VSolid * DDG4SolidConverter::reflected ( const DDSolid &  solid )

staticprivate

Definition at line 271 of file DDG4SolidConverter.cc.

References DDG4SolidConverter(), LogDebug, DDName::name(), DDBase< N, C >::name(), and DDReflectionSolid::unreflected().

Referenced by DDG4SolidConverter().

                                                              {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: reflected = " << solid;
  G4ReflectedSolid * rs = nullptr;
  DDReflectionSolid rfs(solid); 
  if (rfs) {    
    rs = new G4ReflectedSolid(solid.name().name(), 
                              DDG4SolidConverter().convert(rfs.unreflected()), 
                  z_reflection);
    
  } // else ?
  return rs;
}

static G4VSolid* DDG4SolidConverter::shapeless ( const DDSolid &  )

staticprivate

G4VSolid * DDG4SolidConverter::sphere ( const DDSolid &  solid )

staticprivate

Definition at line 569 of file DDG4SolidConverter.cc.

References DDSphere::deltaPhi(), DDSphere::deltaTheta(), DDSphere::innerRadius(), LogDebug, DDName::name(), DDBase< N, C >::name(), DDSphere::outerRadius(), DDSphere::startPhi(), and DDSphere::startTheta().

Referenced by DDG4SolidConverter().

                                                           {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: sphere = " << solid;
  DDSphere sp(solid);
  return new G4Sphere(solid.name().name(), sp.innerRadius(),
              sp.outerRadius(),
              sp.startPhi(),
              sp.deltaPhi(),
              sp.startTheta(),
              sp.deltaTheta());
}

G4VSolid * DDG4SolidConverter::subtraction ( const DDSolid &  solid )

staticprivate

Definition at line 314 of file DDG4SolidConverter.cc.

References DDG4SolidConverter(), LogDebug, DDName::name(), DDBase< N, C >::name(), DDRotation::rotation(), DDBooleanSolid::rotation(), DDBooleanSolid::solidA(), DDBooleanSolid::solidB(), and DDBooleanSolid::translation().

Referenced by DDG4SolidConverter().

                                                                {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: subtraction = " << solid;
  G4SubtractionSolid * us = nullptr;
  DDBooleanSolid bs(solid);
  if (bs) {
    G4VSolid * sa = DDG4SolidConverter().convert(bs.solidA());
    G4VSolid * sb = DDG4SolidConverter().convert(bs.solidB());
    LogDebug("SimG4CoreGeometry") << " name:" << solid.name() << " t=" << bs.translation() << std::flush;
    //       stringstream sst;
    //       bs.rotation().rotation()->inverse().print(sst);
    //       LogDebug("SimG4CoreGeometry") << " " << sst.str() << std::flush;
    LogDebug("SimG4CoreGeometry") << " " << bs.rotation().rotation()->Inverse() << std::flush;
    std::vector<double> tdbl(9);
    bs.rotation().rotation()->Inverse().GetComponents(tdbl.begin(), tdbl.end());
    CLHEP::HepRep3x3 temprep(tdbl[0], tdbl[1], tdbl[2], tdbl[3], tdbl[4], tdbl[5], tdbl[6], tdbl[7], tdbl[8]);
    CLHEP::Hep3Vector temphvec(bs.translation().X(), bs.translation().Y(), bs.translation().Z()); 
    us = new G4SubtractionSolid(solid.name().name(),
                sa,
                sb,
                new CLHEP::HepRotation(temprep),
                temphvec);
  }
  return us;       
}

G4VSolid * DDG4SolidConverter::torus ( const DDSolid &  solid )

staticprivate

Definition at line 255 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), and DDBase< N, C >::name().

Referenced by DDG4SolidConverter().

                                                          {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: torus = " << solid; 
  return new G4Torus(solid.name().name(), (*par_)[0], // rmin
             (*par_)[1], // rmax
             (*par_)[2], // Rtor
             (*par_)[3], // phiStart
             (*par_)[4]);// deltaPhi
}

G4VSolid * DDG4SolidConverter::trap ( const DDSolid &  solid )

staticprivate

Definition at line 94 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), and DDBase< N, C >::name().

Referenced by DDG4SolidConverter(), and pseudotrap().

                                                         {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: trap = " << solid;
  return new G4Trap(solid.name().name(), (*par_)[0],  // pDz
            (*par_)[1],  // theta
            (*par_)[2],  // phi
            (*par_)[3],  // y1
            (*par_)[4],  // x1
            (*par_)[5],  // x2
            (*par_)[6],  // alpha1
            (*par_)[7],  // y2
            (*par_)[8],  // x3
            (*par_)[9],  // x4
            (*par_)[10]);// alpha2
}

G4VSolid * DDG4SolidConverter::trunctubs ( const DDSolid &  solid )

staticprivate

Definition at line 490 of file DDG4SolidConverter.cc.

References funct::abs(), alpha, box(), funct::cos(), DDTruncTubs::cutAtDelta(), DDTruncTubs::cutAtStart(), DDTruncTubs::cutInside(), hiPixelPairStep_cff::deltaPhi, DDTruncTubs::deltaPhi(), Exception, DDName::fullname(), LogDebug, DDSolidShapesName::name(), DDName::name(), dataset::name, DDBase< N, C >::name(), dttmaxenums::R, alignCSCRings::r, mps_fire::result, DDTruncTubs::rIn(), rot, DDTruncTubs::rOut(), DDSolid::shape(), funct::sin(), mathSSE::sqrt(), DDTruncTubs::startPhi(), AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::tt, tubs(), and DDTruncTubs::zHalf().

Referenced by DDG4SolidConverter().

                                                              {
  // truncated tube-section: a boolean subtraction solid:
  //                         from a tube-section a box is subtracted according to the  
  //                         given parameters
  LogDebug("SimG4CoreGeometry") << "MantisConverter: solidshape=" << DDSolidShapesName::name(solid.shape()) << " " << solid;
  LogDebug("SimG4CoreGeometry") << "before";
  DDTruncTubs tt(solid);
  LogDebug("SimG4CoreGeometry") << "after";
  double rIn(tt.rIn()), rOut(tt.rOut()), zHalf(tt.zHalf()),
    startPhi(tt.startPhi()), deltaPhi(tt.deltaPhi()), 
    cutAtStart(tt.cutAtStart()), cutAtDelta(tt.cutAtDelta());
  bool cutInside(bool(tt.cutInside()));
  std::string name=tt.name().name();

  // check the parameters
  if (rIn <= 0 || rOut <=0 || cutAtStart <=0 || cutAtDelta <= 0) {
    throw cms::Exception("DetectorDescriptionFault", "TruncTubs " + std::string(tt.name().fullname()) + ": 0 <= rIn,cutAtStart,rOut,cutAtDelta,rOut violated!");
  }
  if (rIn >= rOut) {
    throw cms::Exception("DetectorDescriptionFault", "TruncTubs " + std::string(tt.name().fullname()) + ": rIn<rOut violated!");
  }
  if (startPhi != 0.) {
    throw cms::Exception("DetectorDescriptionFault", "TruncTubs " + std::string(tt.name().fullname()) + ": startPhi != 0 not supported!");
  }
  //     if (cutInside != false) {
  //       throw cms::Exception("DetectorDescriptionFault", "TruncTubs " + std::string(tt.name()) + " cutInside == true not supported!");
  //     }
  
  startPhi=0.;
  double r(cutAtStart), R(cutAtDelta);
  G4VSolid * result(nullptr);
  G4VSolid * tubs = new G4Tubs(name,rIn,rOut,zHalf,startPhi,deltaPhi);
  LogDebug("SimG4CoreGeometry") << "G4Tubs: " << rIn/CLHEP::cm << ' ' << rOut/CLHEP::cm << ' ' << zHalf/CLHEP::cm << ' ' << startPhi/CLHEP::deg << ' ' << deltaPhi/CLHEP::deg;
  LogDebug("SimG4CoreGeometry") << solid;
  // length & hight of the box 
  double boxX(30.*rOut), boxY(20.*rOut); // exaggerate dimensions - does not matter, it's subtracted!
   
  // width of the box > width of the tubs
  double boxZ(1.1*zHalf);
   
  // angle of the box w.r.t. tubs
  double cath = r-R*cos(deltaPhi);
  double hypo = sqrt(r*r+R*R-2.*r*R*cos(deltaPhi));
  double cos_alpha = cath/hypo;

  double alpha = -acos(cos_alpha);
  LogDebug("SimG4CoreGeometry") << "cath=" << cath/CLHEP::cm;
  LogDebug("SimG4CoreGeometry") << "hypo=" << hypo/CLHEP::cm;
  LogDebug("SimG4CoreGeometry") << "al=" << acos(cath/hypo)/CLHEP::deg;
  LogDebug("SimG4CoreGeometry") << "deltaPhi=" << deltaPhi/CLHEP::deg << "\n"
                << "r=" << r/CLHEP::cm << "\n"
                <<  "R=" << R/CLHEP::cm;

  LogDebug("SimG4CoreGeometry") << "alpha=" << alpha/CLHEP::deg;
    
  // rotationmatrix of box w.r.t. tubs
  G4RotationMatrix * rot = new G4RotationMatrix;
  rot->rotateZ(-alpha);
  LogDebug("SimG4CoreGeometry") << (*rot);

  // center point of the box
  double xBox;
  if (!cutInside) {
    xBox = r+boxY/sin(std::abs(alpha));
  } else {
    xBox = -(boxY/sin(std::abs(alpha))-r);
  }

  G4ThreeVector trans(xBox,0.,0.);
  LogDebug("SimG4CoreGeometry") << "trans=" << trans;

  G4VSolid * box = new G4Box(name,boxX,boxY,boxZ);
  result = new G4SubtractionSolid(name,tubs,box,rot,trans);
      
  return result;

}

G4VSolid * DDG4SolidConverter::tubs ( const DDSolid &  solid )

staticprivate

Definition at line 71 of file DDG4SolidConverter.cc.

References LogDebug, DDName::name(), and DDBase< N, C >::name().

Referenced by DDG4SolidConverter(), pseudotrap(), and trunctubs().

                                                         {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: tubs = " << solid;  
  return new G4Tubs(solid.name().name(), (*par_)[1], // rmin
            (*par_)[2], // rmax
            (*par_)[0], // dzHalf
            (*par_)[3], // phiStart
            (*par_)[4]);// deltaPhi
}

G4VSolid * DDG4SolidConverter::unionsolid ( const DDSolid &  solid )

staticprivate

Definition at line 286 of file DDG4SolidConverter.cc.

References DDG4SolidConverter(), LogDebug, DDName::name(), DDBase< N, C >::name(), DDRotation::rotation(), DDBooleanSolid::rotation(), DDBooleanSolid::solidA(), DDBooleanSolid::solidB(), and DDBooleanSolid::translation().

Referenced by DDG4SolidConverter().

                                                               {
  LogDebug("SimG4CoreGeometry") << "DDG4SolidConverter: unionsolid = " << solid.name();
  G4UnionSolid * us = nullptr;
  DDBooleanSolid bs(solid);
  if (bs) {
    LogDebug("SimG4CoreGeometry") << "SolidA=" << bs.solidA();
    G4VSolid * sa = DDG4SolidConverter().convert(bs.solidA());
    LogDebug("SimG4CoreGeometry") << "SolidB=" << bs.solidB();
    G4VSolid * sb = DDG4SolidConverter().convert(bs.solidB());
    LogDebug("SimG4CoreGeometry") << " name:" << solid.name() << " t=" << bs.translation() << std::flush;
    LogDebug("SimG4CoreGeometry") << " " << bs.rotation().rotation()->Inverse() << std::flush;
    std::vector<double> tdbl(9);
    bs.rotation().rotation()->Inverse().GetComponents(tdbl.begin(), tdbl.end());
    CLHEP::HepRep3x3 temprep(tdbl[0], tdbl[1], tdbl[2], tdbl[3], tdbl[4], tdbl[5], tdbl[6], tdbl[7], tdbl[8]);
    CLHEP::Hep3Vector temphvec(bs.translation().X(), bs.translation().Y(), bs.translation().Z()); 
    us = new G4UnionSolid(solid.name().name(),
              sa,
              sb,
              new CLHEP::HepRotation(temprep),
              temphvec);
    
  } // else?
  return us;       
}

Friends And Related Function Documentation

friend class testPseudoTrap

friend

Definition at line 55 of file DDG4SolidConverter.h.

friend class testTruncTubs

friend

Definition at line 54 of file DDG4SolidConverter.h.

Member Data Documentation

std::map<DDSolidShape,FNPTR> DDG4SolidConverter::convDispatch_

private

Definition at line 50 of file DDG4SolidConverter.h.

Referenced by convert(), and DDG4SolidConverter().

const std::vector< double > * DDG4SolidConverter::par_ = 0

staticprivate

Definition at line 49 of file DDG4SolidConverter.h.

Referenced by convert(), polycone_rrz(), polyhedra_rrz(), and polyhedra_rz().

G4RotationMatrix * DDG4SolidConverter::rot = 0

staticprivate

Definition at line 52 of file DDG4SolidConverter.h.

Referenced by pseudotrap(), and trunctubs().