DDG4Builder Class Reference

#include <DDG4Builder.h>

Public Member Functions
G4LogicalVolume *	BuildGeometry (SensitiveDetectorCatalog &)
	DDG4Builder (const DDCompactView *, G4LogicalVolumeToDDLogicalPartMap &, bool check)
	~DDG4Builder ()

Private Member Functions
G4LogicalVolume *	convertLV (const DDLogicalPart &dLogical)
G4Material *	convertMaterial (const DDMaterial &dMaterial)
G4VSolid *	convertSolid (const DDSolid &dSolid)
double	getDouble (const std::string &s, const DDLogicalPart &dLogical)
int	getInt (const std::string &s, const DDLogicalPart &dLogical)

Private Attributes
bool	check_
const DDCompactView *	compactView_
std::map< DDLogicalPart, G4LogicalVolume * >	logs_
G4LogicalVolumeToDDLogicalPartMap &	map_
std::map< DDMaterial, G4Material * >	mats_
DDG4SolidConverter *	solidConverter_
std::map< DDSolid, G4VSolid * >	sols_
DDG4DispContainer *	theVectorOfDDG4Dispatchables_

Detailed Description

Definition at line 23 of file DDG4Builder.h.

Constructor & Destructor Documentation

DDG4Builder()

DDG4Builder::DDG4Builder	(	const DDCompactView *	cpv,
		G4LogicalVolumeToDDLogicalPartMap &	lvmap,
		bool	check
	)

Definition at line 25 of file DDG4Builder.cc.

References theVectorOfDDG4Dispatchables_.

    : solidConverter_(new DDG4SolidConverter), compactView_(cpv), map_(lvmap), check_(check) {
  theVectorOfDDG4Dispatchables_ = new DDG4DispContainer();
}

~DDG4Builder()

DDG4Builder::~DDG4Builder

(

)

Definition at line 30 of file DDG4Builder.cc.

References solidConverter_.

{ delete solidConverter_; }

Member Function Documentation

BuildGeometry()

G4LogicalVolume * DDG4Builder::BuildGeometry

(

SensitiveDetectorCatalog &

catalog

)

Definition at line 97 of file DDG4Builder.cc.

References getEcalConditions_orcoffint2r_cff::catalog, check_, compactView_, conv, convertLV(), DDBase< N, C >::ddname(), submitPVResolutionJobs::err, Exception, getInt(), DDCompactView::graph(), mps_fire::i, DDMapper< KeyType, ValueType >::insert(), DDBase< N, C >::isDefined(), logs_, map_, hltrates_dqm_sourceclient-live_cfg::offset, eostools::rm(), DDCompactView::root(), AlCaHLTBitMon_QueryRunRegistry::string, makeGlobalPositionRcd_cfg::tag, groupFilesInBlocks::temp, theVectorOfDDG4Dispatchables_, DDBase< N, C >::toString(), x, y, and z.

Referenced by DDDWorld::DDDWorld().

                                                                             {
  G4ReflectionFactory *refFact = G4ReflectionFactory::Instance();
  refFact->SetScalePrecision(100. * refFact->GetScalePrecision());

  using Graph = DDCompactView::Graph;
  const auto &gra = compactView_->graph();
  using adjl_iterator = Graph::const_adj_iterator;
  adjl_iterator git = gra.begin();
  adjl_iterator gend = gra.end();

  Graph::index_type i = 0;
  for (; git != gend; ++git) {
    const DDLogicalPart &ddLP = gra.nodeData(git);
    if (!(ddLP.isDefined().second)) {
      throw cms::Exception("SimG4CoreGeometry",
                           " DDG4Builder::BuildGeometry() has encountered an "
                           "undefined DDLogicalPart named " +
                               ddLP.toString());
    }
    G4LogicalVolume *g4LV = convertLV(ddLP);
    ++i;
    if (!git->empty()) {
      // ask for children of ddLP
      Graph::edge_list::const_iterator cit = git->begin();
      Graph::edge_list::const_iterator cend = git->end();
      for (; cit != cend; ++cit) {
        // fetch specific data
        const DDLogicalPart &ddcurLP = gra.nodeData(cit->first);
        if (!ddcurLP.isDefined().second) {
          std::string err = " DDG4Builder::BuildGeometry() in processing \"children\" has ";
          err += "encountered an undefined DDLogicalPart named " + ddcurLP.toString() + " is a child of " +
                 ddLP.toString();
          throw cms::Exception("SimG4CoreGeometry", err);
        }
        int offset = getInt("CopyNoOffset", ddcurLP);
        int tag = getInt("CopyNoTag", ddcurLP);
        DDRotationMatrix rm(gra.edgeData(cit->second)->rot());
        DD3Vector x, y, z;
        rm.GetComponents(x, y, z);
        if ((x.Cross(y)).Dot(z) < 0)
          edm::LogVerbatim("SimG4CoreGeometry")
              << "DDG4Builder: Reflection: " << gra.edgeData(cit->second)->ddrot()
              << ">>Placement d=" << gra.nodeData(cit->first).ddname() << " m=" << ddLP.ddname()
              << " cp=" << gra.edgeData(cit->second)->copyno() << " r=" << gra.edgeData(cit->second)->ddrot().ddname();
        G4ThreeVector tempTran(gra.edgeData(cit->second)->trans().X(),
                               gra.edgeData(cit->second)->trans().Y(),
                               gra.edgeData(cit->second)->trans().Z());
        G4Translate3D transl = tempTran;
        CLHEP::HepRep3x3 temp(x.X(), x.Y(), x.Z(), y.X(), y.Y(), y.Z(), z.X(), z.Y(), z.Z());  // matrix
        CLHEP::HepRotation hr(temp);
        edm::LogVerbatim("SimG4CoreGeometry")
            << "Position " << gra.nodeData(cit->first).name().name() << ":"
            << gra.edgeData(cit->second)->copyno() + offset + tag << " in " << g4LV->GetName() << " at " << tempTran
            << " with rotation matrix (" << x.X() << ", " << x.Y() << ", " << x.Z() << ", " << y.X() << ", " << y.Y()
            << ", " << y.Z() << ", " << z.X() << ", " << z.Y() << ", " << z.Z() << ")";

        // G3 convention of defining rot-matrices ...
        G4Transform3D trfrm = transl * G4Rotate3D(hr.inverse());  //.inverse();

        refFact->Place(trfrm,  // transformation containing a possible reflection
                       gra.nodeData(cit->first).name().name(),
                       convertLV(gra.nodeData(cit->first)),                 // daugther
                       g4LV,                                                // mother
                       false,                                               // 'ONLY'
                       gra.edgeData(cit->second)->copyno() + offset + tag,  // copy number
                       check_);
      }  // iterate over children
    }    // if (children)
  }      // iterate over graph nodes

  // Looking for in the G4ReflectionFactory secretly created reflected
  // G4LogicalVolumes
  std::map<DDLogicalPart, G4LogicalVolume *>::const_iterator ddg4_it = logs_.begin();
  for (; ddg4_it != logs_.end(); ++ddg4_it) {
    G4LogicalVolume *reflLogicalVolume = refFact->GetReflectedLV(ddg4_it->second);
    if (reflLogicalVolume) {
      DDLogicalPart ddlv = ddg4_it->first;
      map_.insert(reflLogicalVolume, ddlv);
      DDG4Dispatchable *disp = new DDG4Dispatchable(&(ddg4_it->first), reflLogicalVolume);
      theVectorOfDDG4Dispatchables_->push_back(disp);
      edm::LogVerbatim("SimG4CoreGeometry")
          << "DDG4Builder: dd=" << ddlv.ddname() << " g4=" << reflLogicalVolume->GetName();
    }
  }

  G4LogicalVolume *world = logs_[compactView_->root()];

  //
  //  needed for building sensitive detectors
  //
  DDG4SensitiveConverter conv;
  conv.upDate(*theVectorOfDDG4Dispatchables_, catalog);

  return world;
}

convertLV()

G4LogicalVolume * DDG4Builder::convertLV ( const DDLogicalPart &  dLogical )

private

Definition at line 32 of file DDG4Builder.cc.

References convertMaterial(), convertSolid(), DDMapper< KeyType, ValueType >::insert(), logs_, map_, mps_fire::result, and theVectorOfDDG4Dispatchables_.

Referenced by BuildGeometry().

                                                                 {
  edm::LogVerbatim("SimG4CoreGeometry") << "DDG4Builder::convertLV(): DDLogicalPart = " << part;
  G4LogicalVolume *result = logs_[part];
  if (!result) {
    G4VSolid *g4s = convertSolid(part.solid());
    G4Material *g4m = convertMaterial(part.material());
    result = new G4LogicalVolume(g4s, g4m, part.name().name());
    map_.insert(result, part);
    DDG4Dispatchable *disp = new DDG4Dispatchable(&part, result);
    theVectorOfDDG4Dispatchables_->push_back(disp);
    edm::LogVerbatim("SimG4CoreGeometry")
        << "DDG4Builder::convertLV(): new G4LogicalVolume " << part.name().name()
        << "\nDDG4Builder: newEvent: dd=" << part.ddname() << " g4=" << result->GetName();
    logs_[part] = result;  // DDD -> GEANT4
  }
  return result;
}

convertMaterial()

G4Material * DDG4Builder::convertMaterial ( const DDMaterial &  dMaterial )

private

Definition at line 59 of file DDG4Builder.cc.

References DDMaterial::a(), HltBtagPostValidation_cff::c, DDMaterial::constituent(), DDMaterial::density(), Exception, g, mps_fire::i, mats_, DDName::name(), DDBase< N, C >::name(), DDMaterial::noOfConstituents(), mps_fire::result, DDBase< N, C >::toString(), and DDMaterial::z().

Referenced by convertLV().

                                                                   {
  edm::LogVerbatim("SimG4CoreGeometry") << "DDDetConstr::ConvertMaterial: material=" << material;
  G4Material *result = nullptr;
  if (material) {
    // only if it's a valid DDD-material
    if ((result = mats_[material])) {
      edm::LogVerbatim("SimG4CoreGeometry") << "  is already converted";
      return result;
    }
  } else {
    // only if it's NOT a valid DDD-material
    throw cms::Exception("SimG4CoreGeometry",
                         " material is not valid from the Detector Description: " + material.toString());
  }
  int c = 0;
  if ((c = material.noOfConstituents())) {
    // it's a composite material
    edm::LogVerbatim("SimG4CoreGeometry")
        << "  creating a G4-composite material. c=" << c << " d=" << material.density() / CLHEP::g * CLHEP::mole;
    result = new G4Material(material.name().name(), material.density(), c);
    for (int i = 0; i < c; ++i) {
      // recursive building of constituents
      edm::LogVerbatim("SimG4CoreGeometry")
          << "  adding the composite=" << material.name() << " fm=" << material.constituent(i).second;
      result->AddMaterial(convertMaterial(material.constituent(i).first),
                          material.constituent(i).second);  // fractionmass
    }
  } else {
    // it's an elementary material
    edm::LogVerbatim("SimG4CoreGeometry") << "  building an elementary material"
                                          << " z=" << material.z() << " a=" << material.a() / CLHEP::g * CLHEP::mole
                                          << " d=" << material.density() / CLHEP::g * CLHEP::cm3;
    result = new G4Material(material.name().name(), material.z(), material.a(), material.density());
  }
  mats_[material] = result;
  return result;
}

convertSolid()

G4VSolid * DDG4Builder::convertSolid ( const DDSolid &  dSolid )

private

Definition at line 50 of file DDG4Builder.cc.

References DDG4SolidConverter::convert(), mps_fire::result, solidConverter_, and sols_.

Referenced by convertLV().

                                                        {
  G4VSolid *result = sols_[solid];
  if (!result) {
    result = solidConverter_->convert(solid);
    sols_[solid] = result;
  }
  return result;
}

getDouble()

double DDG4Builder::getDouble ( const std::string &  s, const DDLogicalPart &  dLogical  )

private

Definition at line 215 of file DDG4Builder.cc.

References DDfetch(), Exception, mps_fire::result, contentValuesCheck::ss, AlCaHLTBitMon_QueryRunRegistry::string, groupFilesInBlocks::temp, unit(), findQualityFiles::v, and heppy_batch::val.

                                                                            {
  DDValue val(ss);
  std::vector<const DDsvalues_type *> result = part.specifics();
  bool foundIt = false;
  for (auto stype : result) {
    foundIt = DDfetch(stype, val);
    if (foundIt)
      break;
  }
  if (foundIt) {
    std::vector<std::string> temp = val.strings();
    if (temp.size() != 1) {
      throw cms::Exception("SimG4CoreGeometry",
                           " DDG4Builder::getDouble() Problem with Region tags "
                           "- one and only one allowed: " +
                               ss);
    }
    double v;
    std::string unit;
    std::istringstream is(temp[0]);
    is >> v >> unit;
    v = v * G4UnitDefinition::GetValueOf(unit.substr(1, unit.size()));
    return v;
  } else
    return 0;
}

getInt()

int DDG4Builder::getInt ( const std::string &  s, const DDLogicalPart &  dLogical  )

private

Definition at line 193 of file DDG4Builder.cc.

References DDfetch(), Exception, createfilelist::int, mps_fire::result, contentValuesCheck::ss, groupFilesInBlocks::temp, and heppy_batch::val.

Referenced by BuildGeometry().

                                                                      {
  DDValue val(ss);
  std::vector<const DDsvalues_type *> result = part.specifics();
  bool foundIt = false;
  for (auto stype : result) {
    foundIt = DDfetch(stype, val);
    if (foundIt)
      break;
  }
  if (foundIt) {
    std::vector<double> temp = val.doubles();
    if (temp.size() != 1) {
      throw cms::Exception("SimG4CoreGeometry",
                           " DDG4Builder::getInt() Problem with Region tags - "
                           "one and only one allowed: " +
                               ss);
    }
    return int(temp[0]);
  } else
    return 0;
}

Member Data Documentation

check_

bool DDG4Builder::check_

private

Definition at line 44 of file DDG4Builder.h.

Referenced by BuildGeometry().

compactView_

const DDCompactView* DDG4Builder::compactView_

private

Definition at line 41 of file DDG4Builder.h.

Referenced by BuildGeometry().

logs_

std::map<DDLogicalPart, G4LogicalVolume *> DDG4Builder::logs_

private

Definition at line 39 of file DDG4Builder.h.

Referenced by BuildGeometry(), and convertLV().

map_

G4LogicalVolumeToDDLogicalPartMap& DDG4Builder::map_

private

Definition at line 42 of file DDG4Builder.h.

Referenced by BuildGeometry(), and convertLV().

mats_

std::map<DDMaterial, G4Material *> DDG4Builder::mats_

private

Definition at line 37 of file DDG4Builder.h.

Referenced by convertMaterial().

solidConverter_

DDG4SolidConverter* DDG4Builder::solidConverter_

private

Definition at line 36 of file DDG4Builder.h.

Referenced by convertSolid(), and ~DDG4Builder().

sols_

std::map<DDSolid, G4VSolid *> DDG4Builder::sols_

private

Definition at line 38 of file DDG4Builder.h.

Referenced by convertSolid().

theVectorOfDDG4Dispatchables_

DDG4DispContainer* DDG4Builder::theVectorOfDDG4Dispatchables_

private

Definition at line 43 of file DDG4Builder.h.

Referenced by BuildGeometry(), convertLV(), and DDG4Builder().