33 #include "CLHEP/Geometry/Point3D.h"
34 #include "CLHEP/Geometry/Vector3D.h"
35 #include "CLHEP/Units/GlobalSystemOfUnits.h"
36 #include "CLHEP/Units/GlobalPhysicalConstants.h"
55 std::vector<edm::ESGetToken<HGCalDDDConstants, IdealGeometryRecord>>
geomToken_;
80 geometrySource_ =
cfg.getUntrackedParameter<std::vector<std::string>>(
"geometrySource");
84 esConsumes<HGCalDDDConstants, IdealGeometryRecord, edm::Transition::BeginRun>(
edm::ESInputTag{
"",
name}));
170 std::vector<float> hitVtxX = infoLayer->
hitvtxX();
171 std::vector<float> hitVtxY = infoLayer->
hitvtxY();
172 std::vector<float> hitVtxZ = infoLayer->
hitvtxZ();
173 std::vector<unsigned int> hitDet = infoLayer->
hitDets();
174 std::vector<unsigned int> hitIdx = infoLayer->
hitIndex();
177 std::vector<float> edepLayerEE = infoLayer->
eehgcEdep();
178 std::vector<float> edepLayerHE = infoLayer->
hefhgcEdep();
179 std::vector<float> edepLayerHB = infoLayer->
hebhgcEdep();
182 for (
i = 0;
i < edepLayerEE.size();
i++) {
186 for (
i = 0;
i < edepLayerHE.size();
i++) {
190 for (
i = 0;
i < edepLayerHB.size();
i++) {
200 for (
unsigned int i = 0;
i < hitVtxX.size();
i++) {
201 hitVtxX.at(
i) =
mmtocm * hitVtxX.at(
i);
202 hitVtxY.at(
i) =
mmtocm * hitVtxY.at(
i);
203 hitVtxZ.at(
i) =
mmtocm * hitVtxZ.at(
i);
209 std::pair<float, float>
xy;
211 int subdet, wafer, celltype, cell;
252 }
else if (
dtype == 1) {