DDLTubs.cc

Go to the documentation of this file.

#include "DetectorDescription/Parser/src/DDLTubs.h"
#include "DetectorDescription/Core/interface/DDSolid.h"
#include "DetectorDescription/Core/interface/ClhepEvaluator.h"
#include "DataFormats/Math/interface/GeantUnits.h"
#include "DetectorDescription/Parser/interface/DDLElementRegistry.h"
#include "DetectorDescription/Parser/src/DDLSolid.h"
#include "DetectorDescription/Parser/src/DDXMLElement.h"

#include <map>
#include <utility>

class DDCompactView;

using namespace geant_units::operators;

DDLTubs::DDLTubs(DDLElementRegistry* myreg) : DDLSolid(myreg) {}

// Upon encountering the end of a Tubs element, call DDCore.
void DDLTubs::processElement(const std::string& name, const std::string& nmspace, DDCompactView& cpv) {
  ClhepEvaluator& ev = myRegistry_->evaluator();
  DDXMLAttribute atts = getAttributeSet();
  bool cutInside(false);

  if (atts.find("cutInside") != atts.end()) {
    cutInside = (atts.find("cutInside")->second == "true") ? true : false;
  }

  if (name == "Tubs") {
    DDSolid myTubs = DDSolidFactory::tubs(getDDName(nmspace),
                                          ev.eval(nmspace, atts.find("dz")->second),
                                          ev.eval(nmspace, atts.find("rMin")->second),
                                          ev.eval(nmspace, atts.find("rMax")->second),
                                          ev.eval(nmspace, atts.find("startPhi")->second),
                                          ev.eval(nmspace, atts.find("deltaPhi")->second));
  } else if (name == "Tube") {
    DDSolid myTubs = DDSolidFactory::tubs(getDDName(nmspace),
                                          ev.eval(nmspace, atts.find("dz")->second),
                                          ev.eval(nmspace, atts.find("rMin")->second),
                                          ev.eval(nmspace, atts.find("rMax")->second),
                                          0_deg,
                                          360_deg);
  } else if (name == "TruncTubs") {
    DDSolid myTT = DDSolidFactory::truncTubs(getDDName(nmspace),
                                             ev.eval(nmspace, atts.find("zHalf")->second),
                                             ev.eval(nmspace, atts.find("rMin")->second),
                                             ev.eval(nmspace, atts.find("rMax")->second),
                                             ev.eval(nmspace, atts.find("startPhi")->second),  //0. // startPhi
                                             ev.eval(nmspace, atts.find("deltaPhi")->second),
                                             ev.eval(nmspace, atts.find("cutAtStart")->second),
                                             ev.eval(nmspace, atts.find("cutAtDelta")->second),
                                             cutInside);  // cutInside
  } else if (name == "CutTubs") {
    DDSolid myTubs = DDSolidFactory::cuttubs(getDDName(nmspace),
                                             ev.eval(nmspace, atts.find("dz")->second),
                                             ev.eval(nmspace, atts.find("rMin")->second),
                                             ev.eval(nmspace, atts.find("rMax")->second),
                                             ev.eval(nmspace, atts.find("startPhi")->second),
                                             ev.eval(nmspace, atts.find("deltaPhi")->second),
                                             ev.eval(nmspace, atts.find("lx")->second),
                                             ev.eval(nmspace, atts.find("ly")->second),
                                             ev.eval(nmspace, atts.find("lz")->second),
                                             ev.eval(nmspace, atts.find("tx")->second),
                                             ev.eval(nmspace, atts.find("ty")->second),
                                             ev.eval(nmspace, atts.find("tz")->second));
  } else {
    std::string msg = "\nDDLTubs::processElement could not process element.";
    throwError(msg);
  }
  DDLSolid::setReference(nmspace, cpv);
}