Merge pull request #34391 from jkiesele/rechittools_CMSSW_12_0_0_pre2

Added scintillator information to RecHitTools
cms-sw · Jul 9, 2021 · 76b9b71 · 76b9b71
2 parents 91b871c + 3046eb5
commit 76b9b71
Show file tree

Hide file tree

Showing 2 changed files with 15 additions and 0 deletions.
diff --git a/RecoLocalCalo/HGCalRecAlgos/interface/RecHitTools.h b/RecoLocalCalo/HGCalRecAlgos/interface/RecHitTools.h
@@ -37,6 +37,8 @@ namespace hgcal {
     std::float_t getRadiusToSide(const DetId&) const;
     int getSiThickIndex(const DetId&) const;
 
+    std::pair<float, float> getScintDEtaDPhi(const DetId&) const;
+
     unsigned int getLayer(DetId::Detector type, bool nose = false) const;
     unsigned int getLayer(ForwardSubdetector type) const;
     unsigned int getLayer(const DetId&) const;
@@ -47,6 +49,7 @@ namespace hgcal {
     bool isHalfCell(const DetId&) const;
 
     bool isSilicon(const DetId&) const;
+    bool isScintillator(const DetId&) const;
 
     bool isOnlySilicon(const unsigned int layer) const;
 

diff --git a/RecoLocalCalo/HGCalRecAlgos/src/RecHitTools.cc b/RecoLocalCalo/HGCalRecAlgos/src/RecHitTools.cc
@@ -219,6 +219,16 @@ int RecHitTools::getSiThickIndex(const DetId& id) const {
   return thickIndex;
 }
 
+std::pair<float, float> RecHitTools::getScintDEtaDPhi(const DetId& id) const {
+  if (!isScintillator(id)) {
+    LogDebug("getScintDEtaDPhi::InvalidScintDetid")
+        << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal scintillator!";
+    return {0.f, 0.f};
+  }
+  auto cellGeom = getSubdetectorGeometry(id)->getGeometry(id);
+  return {cellGeom->etaSpan(), cellGeom->phiSpan()};
+}
+
 std::float_t RecHitTools::getRadiusToSide(const DetId& id) const {
   auto geom = getSubdetectorGeometry(id);
   std::float_t size(std::numeric_limits<std::float_t>::max());
@@ -417,6 +427,8 @@ bool RecHitTools::isSilicon(const DetId& id) const {
           (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)));
 }
 
+bool RecHitTools::isScintillator(const DetId& id) const { return id.det() == DetId::HGCalHSc; }
+
 bool RecHitTools::isOnlySilicon(const unsigned int layer) const {
   // HFnose TODO
   bool isonlysilicon = (layer % bhLastLayer_) < bhOffset_;