cms::DDCMSDetElementCreator Class Reference

Inheritance diagram for cms::DDCMSDetElementCreator:

Classes
struct	Count
struct	Data

Public Member Functions
	DDCMSDetElementCreator (dd4hep::Detector &)
int	operator() (dd4hep::PlacedVolume volume, int level) override
	Callback to output PlacedVolume information of an single Placement. More...
int	process (dd4hep::PlacedVolume volume, int level, bool recursive) override
	Callback to output PlacedVolume information of an entire Placement. More...
	~DDCMSDetElementCreator () override

Private Types
using	Counters = std::map< dd4hep::DetElement, Count >
using	Detectors = std::map< std::string, dd4hep::DetElement >
using	LeafCount = std::map< std::pair< dd4hep::DetElement, int >, std::pair< int, int > >
using	VolumeStack = std::vector< Data >

Private Member Functions
dd4hep::DetElement	addSubdetector (const std::string &nam, dd4hep::PlacedVolume volume, bool valid)
dd4hep::DetElement	createElement (const char *debugTag, dd4hep::PlacedVolume volume, int id)
void	createTopLevelDetectors (dd4hep::PlacedVolume volume)

Private Attributes
std::map< dd4hep::PlacedVolume, std::pair< int, int > >	m_allPlacements
Counters	m_counters
dd4hep::DetElement	m_currentDetector
dd4hep::SensitiveDetector	m_currentSensitive
dd4hep::Detector &	m_description
LeafCount	m_leafCount
dd4hep::Atom	m_silicon
VolumeStack	m_stack
Detectors	m_subdetectors
dd4hep::DetElement	m_tracker

Detailed Description

Definition at line 19 of file DDCMSDetElementCreator.cc.

Member Typedef Documentation

Counters

using cms::DDCMSDetElementCreator::Counters = std::map<dd4hep::DetElement, Count>

private

Definition at line 60 of file DDCMSDetElementCreator.cc.

Detectors

using cms::DDCMSDetElementCreator::Detectors = std::map<std::string, dd4hep::DetElement>

private

Definition at line 59 of file DDCMSDetElementCreator.cc.

LeafCount

using cms::DDCMSDetElementCreator::LeafCount = std::map<std::pair<dd4hep::DetElement, int>, std::pair<int, int> >

private

Definition at line 61 of file DDCMSDetElementCreator.cc.

VolumeStack

using cms::DDCMSDetElementCreator::VolumeStack = std::vector<Data>

private

Definition at line 62 of file DDCMSDetElementCreator.cc.

Constructor & Destructor Documentation

DDCMSDetElementCreator()

DDCMSDetElementCreator::DDCMSDetElementCreator

(

dd4hep::Detector &

desc

)

Definition at line 89 of file DDCMSDetElementCreator.cc.

                                                                   : m_description(desc) {
  dd4hep::DetectorHelper helper(m_description);
  m_silicon = helper.element("SI");
  if (!m_silicon.isValid()) {
    except("DDCMSDetElementCreator", "++ Failed to extract SILICON from the element table.");
  }
  m_stack.reserve(32);
}

References m_description, m_silicon, and m_stack.

~DDCMSDetElementCreator()

DDCMSDetElementCreator::~DDCMSDetElementCreator ( )

override

Create ID Descriptors and readout configurations

Definition at line 98 of file DDCMSDetElementCreator.cc.

                                                {
  Count total;
  stringstream str, id_str;
 
  printout(INFO, "DDCMSDetElementCreator", "+++++++++++++++ Summary of sensitve elements  ++++++++++++++++++++++++");
  for (const auto& c : m_counters) {
    printout(INFO,
             "DDCMSDetElementCreator",
             "++ Summary: SD: %-24s %7d DetElements %7d sensitives out of %7d volumes",
             (c.first.name() + string(":")).c_str(),
             c.second.elements,
             c.second.sensitives,
             c.second.volumes);
    total.elements += c.second.elements;
    total.sensitives += c.second.sensitives;
    total.volumes += c.second.volumes;
  }
  printout(INFO,
           "DDCMSDetElementCreator",
           "++ Summary:     %-24s %7d DetElements %7d sensitives out of %7d volumes",
           "Grand Total:",
           total.elements,
           total.sensitives,
           total.volumes);
  printout(INFO, "DDCMSDetElementCreator", "+++++++++++++++ Summary of geometry depth analysis  ++++++++++++++++++");
  int totalCount = 0, totalDepth = 0;
  map<dd4hep::DetElement, vector<pair<int, int> > > fields;
  for (const auto& l : m_leafCount) {
    dd4hep::DetElement de = l.first.first;
    printout(INFO,
             "DDCMSDetElementCreator",
             "++ Summary: SD: %-24s system:%04X Lvl:%3d Sensitives: %6d [Max: %6d].",
             (de.name() + string(":")).c_str(),
             de.id(),
             l.first.second,
             l.second.second,
             l.second.first);
    fields[de].push_back(make_pair(l.first.second, l.second.first));
    totalDepth += l.second.second;
    ++totalCount;
  }
  printout(INFO, "DDCMSDetElementCreator", "++ Summary:     %-24s  %d.", "Total DetElements:", totalCount);
  printout(INFO, "DDCMSDetElementCreator", "+++++++++++++++ Readout structure generation  ++++++++++++++++++++++++");
  str << endl;
  for (const auto& f : fields) {
    string roName = f.first.name() + string("Hits");
    int num_bits = 8;
    id_str.str("");
    id_str << "system:" << num_bits;
    for (const auto& q : f.second) {
      int bits = 0;
      if (q.second < 1 << 0)
        bits = 1;
      else if (q.second < 1 << 1)
        bits = 1;
      else if (q.second < 1 << 2)
        bits = 2;
      else if (q.second < 1 << 3)
        bits = 3;
      else if (q.second < 1 << 4)
        bits = 4;
      else if (q.second < 1 << 5)
        bits = 5;
      else if (q.second < 1 << 6)
        bits = 6;
      else if (q.second < 1 << 7)
        bits = 7;
      else if (q.second < 1 << 8)
        bits = 8;
      else if (q.second < 1 << 9)
        bits = 9;
      else if (q.second < 1 << 10)
        bits = 10;
      else if (q.second < 1 << 11)
        bits = 11;
      else if (q.second < 1 << 12)
        bits = 12;
      else if (q.second < 1 << 13)
        bits = 13;
      else if (q.second < 1 << 14)
        bits = 14;
      else if (q.second < 1 << 15)
        bits = 15;
      bits += 1;
      id_str << ",Lv" << q.first << ":" << bits;
      num_bits += bits;
    }
    string idspec = id_str.str();
    str << "<readout name=\"" << roName << "\">" << endl
        << "\t<id>" << idspec << "</id>  <!-- Number of bits: " << num_bits << " -->" << endl
        << "</readout>" << endl;
 
    IDDescriptor dsc(roName, idspec);
    m_description.addIDSpecification(dsc);
    Readout ro(roName);
    ro.setIDDescriptor(dsc);
    m_description.addReadout(ro);
    dd4hep::SensitiveDetector sd = m_description.sensitiveDetector(f.first.name());
    sd.setHitsCollection(ro.name());
    sd.setReadout(ro);
    printout(INFO,
             "DDCMSDetElementCreator",
             "++ Setting up readout for subdetector:%-24s id:%04X",
             f.first.name(),
             f.first.id());
  }
  printout(INFO, "DDCMSDetElementCreator", "+++++++++++++++ ID Descriptor generation  ++++++++++++++++++++++++++++");
  printout(INFO, "", str.str().c_str());
  char volId[32];
  for (auto& p : m_allPlacements) {
    dd4hep::PlacedVolume place = p.first;
    dd4hep::Volume volume = place.volume();
    ::snprintf(volId, sizeof(volId), "Lv%d", p.second.first);
    printout(DEBUG,
             "DDCMSDetElementCreator",
             "++ Set volid (%-24s): %-6s = %3d  -> %s  (%p)",
             volume.isSensitive() ? volume.sensitiveDetector().name() : "Not Sensitive",
             volId,
             p.second.second,
             place.name(),
             place.ptr());
    place.addPhysVolID(volId, p.second.second);
  }
  printout(ALWAYS,
           "DDCMSDetElementCreator",
           "++ Instrumented %ld subdetectors with %d DetElements %d sensitives out of %d volumes and %ld sensitive "
           "placements.",
           fields.size(),
           total.elements,
           total.sensitives,
           total.volumes,
           m_allPlacements.size());
}

References HltBtagPostValidation_cff::c, DEBUG, f, l1GtPatternGenerator_cfi::fields, unpackBuffers-CaloStage2::INFO, cmsLHEtoEOSManager::l, m_allPlacements, m_counters, m_description, m_leafCount, AlCaHLTBitMon_ParallelJobs::p, submitPVResolutionJobs::q, sd, str, AlCaHLTBitMon_QueryRunRegistry::string, dqmMemoryStats::total, and simpleEdmComparison::totalCount.

Member Function Documentation

addSubdetector()

dd4hep::DetElement DDCMSDetElementCreator::addSubdetector ( const std::string &  nam, dd4hep::PlacedVolume  volume, bool  valid  )

private

Definition at line 249 of file DDCMSDetElementCreator.cc.

                                                                      {
  auto idet = m_subdetectors.find(nam);
  if (idet == m_subdetectors.end()) {
    dd4hep::DetElement det(nam, m_subdetectors.size() + 1);
    det.setPlacement(volume);
    if (valid) {
      det.placement().addPhysVolID("system", det.id());
    }
    idet = m_subdetectors.insert(make_pair(nam, det)).first;
    m_description.add(det);
  }
  return idet->second;
}

References m_description, m_subdetectors, and RunInfoPI::valid.

Referenced by createTopLevelDetectors().

createElement()

dd4hep::DetElement DDCMSDetElementCreator::createElement ( const char *  debugTag, dd4hep::PlacedVolume  volume, int  id  )

private

Definition at line 233 of file DDCMSDetElementCreator.cc.

                                                                                               {
  string name = detElementName(volume);
  dd4hep::DetElement det(name, id);
  det.setPlacement(volume);
  return det;
}

References cms::detElementName(), and Skims_PA_cff::name.

Referenced by process().

createTopLevelDetectors()

void DDCMSDetElementCreator::createTopLevelDetectors ( dd4hep::PlacedVolume  volume )

private

Definition at line 239 of file DDCMSDetElementCreator.cc.

                                                                        {
  auto& data = m_stack.back();
  if (m_stack.size() == 2) {  // Main subssystem: tracker:Tracker
    data.element = m_tracker = addSubdetector(cms::detElementName(volume), volume, false);
    m_tracker->SetTitle("compound");
  } else if (m_stack.size() == 3) {  // Main subsystem detector: TIB, TEC, ....
    data.element = m_currentDetector = addSubdetector(cms::detElementName(volume), volume, true);
  }
}

References addSubdetector(), data, cms::detElementName(), m_currentDetector, m_stack, and m_tracker.

Referenced by process().

operator()()

int DDCMSDetElementCreator::operator() ( dd4hep::PlacedVolume  volume, int  level  )

override

Callback to output PlacedVolume information of an single Placement.

Definition at line 265 of file DDCMSDetElementCreator.cc.

                                                                                 {
  double fracSi = volume.volume().material().fraction(m_silicon);
  if (fracSi > 90e-2) {
    Data& data = m_stack.back();
    data.sensitive = true;
    data.hasSensitive = true;
    ++data.volumeCount;
    int idx = volume->GetMotherVolume()->GetIndex(volume.ptr()) + 1;
    auto& cnt = m_leafCount[make_pair(m_currentDetector, volumeLevel)];
    cnt.first = std::max(cnt.first, idx);
    ++cnt.second;
    m_allPlacements[volume] = make_pair(volumeLevel, idx);
    return 1;
  }
  return 0;
}

References data, MillePedeFileConverter_cfg::e, heavyIonCSV_trainingSettings::idx, m_allPlacements, m_currentDetector, m_leafCount, m_silicon, m_stack, and SiStripPI::max.

process()

int DDCMSDetElementCreator::process ( dd4hep::PlacedVolume  volume, int  level, bool  recursive  )

override

Callback to output PlacedVolume information of an entire Placement.

Complete structures if the m_stack size is > 3!

If this volume is sensitve, we must attach a sensitive detector handle

Now the cleanup kicks in....

Definition at line 282 of file DDCMSDetElementCreator.cc.

                                                                                        {
  m_stack.push_back(Data(volume));
  if (m_stack.size() <= 3) {
    createTopLevelDetectors(volume);
  }
  int ret = dd4hep::PlacedVolumeProcessor::process(volume, level, recursive);
 
  if (m_stack.size() > 3) {
    // Note: short-cuts to entries in the m_stack MUST be local and
    // initialized AFTER the call to "process"! The vector may be resized!
    auto& data = m_stack.back();
    auto& parent = m_stack[m_stack.size() - 2];
    auto& counts = m_counters[m_currentDetector];
    if (data.sensitive) {
      if (!m_currentSensitive.isValid()) {
        dd4hep::SensitiveDetector sd = m_description.sensitiveDetector(m_currentDetector.name());
        if (!sd.isValid()) {
          sd = dd4hep::SensitiveDetector(m_currentDetector.name(), "tracker");
          m_currentDetector->flag |= DetElement::Object::HAVE_SENSITIVE_DETECTOR;
          m_description.add(sd);
        }
        m_currentSensitive = sd;
      }
      volume.volume().setSensitiveDetector(m_currentSensitive);
      ++counts.sensitives;
    }
    ++counts.volumes;
    bool added = false;
    if (data.volumeCount > 0) {
      parent.daughterCount += data.volumeCount;
      parent.daughterCount += data.daughterCount;
      data.hasSensitive = true;
    } else {
      parent.daughterCount += data.daughterCount;
      data.hasSensitive = (data.daughterCount > 0);
    }
 
    if (data.hasSensitive) {
      // If we have sensitive elements at this level or below,
      // we must complete the DetElement hierarchy
      if (!data.element.isValid()) {
        data.element = createElement("Element", data.volume, m_currentDetector.id());
        ++counts.elements;
      }
      if (!parent.element.isValid()) {
        parent.element = createElement("Parent ", parent.volume, m_currentDetector.id());
        ++counts.elements;
      }
      printout(DEBUG,
               "DDCMSDetElementCreator",
               "++ Assign detector element: %s (%p, %ld children) to %s (%p) with %ld vols",
               data.element.name(),
               data.element.ptr(),
               data.element.children().size(),
               parent.element.name(),
               parent.element.ptr(),
               data.volumeCount);
 
      // Trickle up the tree only for sensitive pathes. Forget the passive rest
      // This should automatically omit non-sensitive pathes
      parent.hasSensitive = true;
      parent.element.add(data.element);
      added = true;
      // It is simpler to collect the volumes and later assign the volids
      // rather than checking if the volid already exists.
      int volumeLevel = level;
      int idx = data.volume->GetMotherVolume()->GetIndex(data.volume.ptr()) + 1;
      m_allPlacements[data.volume] = make_pair(volumeLevel, idx);  // 1...n
      // Update counters
      auto& cnt_det = m_leafCount[make_pair(m_currentDetector, volumeLevel)];
      cnt_det.first = std::max(cnt_det.first, idx);
      cnt_det.second += 1;
    }
    if (!added && data.element.isValid()) {
      printout(WARNING,
               "MEMORY-LEAK",
               "Level:%3d Orpahaned DetElement:%s Daugthers:%d Parent:%s",
               int(m_stack.size()),
               data.element.name(),
               data.volumeCount,
               parent.volume.name());
    }
  }
  if (m_stack.size() == 3) {
    m_currentSensitive = SensitiveDetector();
    m_currentDetector = DetElement();
    ret = 0;
  }
  m_stack.pop_back();
  return ret;
}

References dqmiodumpmetadata::counts, createElement(), createTopLevelDetectors(), data, DEBUG, heavyIonCSV_trainingSettings::idx, personalPlayback::level, m_allPlacements, m_counters, m_currentDetector, m_currentSensitive, m_description, m_leafCount, m_stack, SiStripPI::max, class-composition::parent, LaserDQM_cfg::process, postprocess-scan-build::recursive, runTheMatrix::ret, sd, and remoteMonitoring_LED_IterMethod_cfg::WARNING.

Member Data Documentation

m_allPlacements

std::map<dd4hep::PlacedVolume, std::pair<int, int> > cms::DDCMSDetElementCreator::m_allPlacements

private

Definition at line 64 of file DDCMSDetElementCreator.cc.

Referenced by operator()(), process(), and ~DDCMSDetElementCreator().

m_counters

Counters cms::DDCMSDetElementCreator::m_counters

private

Definition at line 66 of file DDCMSDetElementCreator.cc.

Referenced by process(), and ~DDCMSDetElementCreator().

m_currentDetector

dd4hep::DetElement cms::DDCMSDetElementCreator::m_currentDetector

private

Definition at line 70 of file DDCMSDetElementCreator.cc.

Referenced by createTopLevelDetectors(), operator()(), and process().

m_currentSensitive

dd4hep::SensitiveDetector cms::DDCMSDetElementCreator::m_currentSensitive

private

Definition at line 71 of file DDCMSDetElementCreator.cc.

Referenced by process().

m_description

dd4hep::Detector& cms::DDCMSDetElementCreator::m_description

private

Definition at line 72 of file DDCMSDetElementCreator.cc.

Referenced by addSubdetector(), DDCMSDetElementCreator(), process(), and ~DDCMSDetElementCreator().

m_leafCount

LeafCount cms::DDCMSDetElementCreator::m_leafCount

private

Definition at line 67 of file DDCMSDetElementCreator.cc.

Referenced by operator()(), process(), and ~DDCMSDetElementCreator().

m_silicon

dd4hep::Atom cms::DDCMSDetElementCreator::m_silicon

private

Definition at line 73 of file DDCMSDetElementCreator.cc.

Referenced by DDCMSDetElementCreator(), and operator()().

m_stack

VolumeStack cms::DDCMSDetElementCreator::m_stack

private

Definition at line 68 of file DDCMSDetElementCreator.cc.

Referenced by createTopLevelDetectors(), DDCMSDetElementCreator(), operator()(), and process().

m_subdetectors

Detectors cms::DDCMSDetElementCreator::m_subdetectors

private

Definition at line 69 of file DDCMSDetElementCreator.cc.

Referenced by addSubdetector().

m_tracker

dd4hep::DetElement cms::DDCMSDetElementCreator::m_tracker

private

Definition at line 70 of file DDCMSDetElementCreator.cc.

Referenced by createTopLevelDetectors().