Geometry/HGCalCommonData/test/HGCalNumberingTester.cc - master - githubmirror/cmssw

Geometry/HGCalCommonData/test/HGCalNumberingTester.cc
233 строки12 KB
Sunanda
Add afew new modes to enable cog computation for wafers as default for the future V19 version of HGCal
03 июл 2024, 10:03
03 июл 2024, 10:0369c5983
// -*- C++ -*-
//
// Package:    HGCalNumberingTester
// Class:      HGCalNumberingTester
//
/**\class HGCalNumberingTester HGCalNumberingTester.cc
 test/HGCalNumberingTester.cc

 Description: <one line class summary>

 Implementation:
     <Notes on implementation>
*/
//
// Original Author:  Sunanda Banerjee
//         Created:  Mon 2014/03/21
//
//

// system include files
#include <array>
#include <fstream>
#include <iostream>
#include <memory>
#include <string>
#include <vector>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "DetectorDescription/Core/interface/DDCompactView.h"
#include "DetectorDescription/Core/interface/DDExpandedView.h"
#include "DetectorDescription/Core/interface/DDSpecifics.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "Geometry/HGCalCommonData/interface/HGCalDDDConstants.h"
#include "Geometry/Records/interface/IdealGeometryRecord.h"

class HGCalNumberingTester : public edm::one::EDAnalyzer<> {
public:
  explicit HGCalNumberingTester(const edm::ParameterSet&);
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

  void beginJob() override {}
  void analyze(edm::Event const& iEvent, edm::EventSetup const&) override;
  void endJob() override {}

private:
  edm::ESGetToken<HGCalDDDConstants, IdealGeometryRecord> dddToken_;
  std::string nameSense_, nameDetector_;
  std::vector<double> positionX_, positionY_;
  int increment_, detType_;
  bool reco_;
};

HGCalNumberingTester::HGCalNumberingTester(const edm::ParameterSet& iC) {
  nameSense_ = iC.getParameter<std::string>("NameSense");
  nameDetector_ = iC.getParameter<std::string>("NameDevice");
  positionX_ = iC.getParameter<std::vector<double> >("LocalPositionX");
  positionY_ = iC.getParameter<std::vector<double> >("LocalPositionY");
  increment_ = iC.getParameter<int>("Increment");
  detType_ = iC.getParameter<int>("DetType");
  reco_ = iC.getParameter<bool>("Reco");

  dddToken_ = esConsumes<HGCalDDDConstants, IdealGeometryRecord>(edm::ESInputTag{"", nameSense_});

  std::string unit("mm");
  if (reco_) {
    for (unsigned int k = 0; k < positionX_.size(); ++k) {
      positionX_[k] /= CLHEP::cm;
      positionY_[k] /= CLHEP::cm;
    }
    unit = "cm";
  } else {
    for (unsigned int k = 0; k < positionX_.size(); ++k) {
      positionX_[k] /= CLHEP::mm;
      positionY_[k] /= CLHEP::mm;
    }
  }
  edm::LogVerbatim("HGCalGeom") << "Test numbering for " << nameDetector_ << " using constants of " << nameSense_
                                << " at " << positionX_.size() << " local positions "
                                << "for every " << increment_ << " layers for DetType " << detType_ << " and  RecoFlag "
                                << reco_;
  for (unsigned int k = 0; k < positionX_.size(); ++k)
    edm::LogVerbatim("HGCalGeom") << "Position[" << k << "] " << positionX_[k] << " " << unit << ", " << positionY_[k]
                                  << " " << unit;
}

void HGCalNumberingTester::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  std::vector<double> vecxy;
  edm::ParameterSetDescription desc;
  desc.add<std::string>("NameSense", "HGCalEESensitive");
  desc.add<std::string>("NameDevice", "HGCal EE");
  desc.add<std::vector<double> >("LocalPositionX", vecxy);
  desc.add<std::vector<double> >("LocalPositionY", vecxy);
  desc.add<int>("Increment", 19);
  desc.add<int>("DetType", 2);
  desc.add<bool>("Reco", false);
  descriptions.add("hgcalNumberingTesterEE", desc);
}

// ------------ method called to produce the data  ------------
void HGCalNumberingTester::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  const HGCalDDDConstants& hgdc = iSetup.getData(dddToken_);
  edm::LogVerbatim("HGCalGeom") << nameDetector_ << " Layers = " << hgdc.layers(reco_)
                                << " Sectors = " << hgdc.sectors() << " Minimum Slope = " << hgdc.minSlope();
  if (detType_ != 0) {
    edm::LogVerbatim("HGCalGeom") << "Minimum Wafer # " << hgdc.waferMin() << " Mamximum Wafer # " << hgdc.waferMax()
                                  << " Wafer counts " << hgdc.waferCount(0) << ":" << hgdc.waferCount(1);
    for (unsigned int i = 0; i < hgdc.layers(true); ++i) {
      int lay = i + 1;
      double z = hgdc.waferZ(lay, reco_);
      edm::LogVerbatim("HGCalGeom") << "Layer " << lay << " Wafers " << hgdc.wafers(lay, 0) << ":"
                                    << hgdc.wafers(lay, 1) << ":" << hgdc.wafers(lay, 2) << " Z " << z << " R "
                                    << hgdc.rangeR(z, reco_).first << ":" << hgdc.rangeR(z, reco_).second;
    }
  }
  edm::LogVerbatim("HGCalGeom") << std::endl;
  std::pair<float, float> xy;
  std::string flg;
  int subsec(0);
  int loff = hgdc.firstLayer();
  double scl = (reco_ ? 10.0 : 1.0);
  for (unsigned int k = 0; k < positionX_.size(); ++k) {
    float localx(positionX_[k]), localy(positionY_[k]);
    for (unsigned int i = 0; i < hgdc.layers(reco_); ++i) {
      double zpos = hgdc.waferZ(i + loff, reco_);
      int zside = (zpos > 0) ? 1 : -1;
      if (detType_ == 1) {
        std::pair<int, int> kxy, lxy;
        kxy = hgdc.assignCell(localx, localy, i + loff, subsec, reco_);
        xy = hgdc.locateCell(kxy.second, i + loff, kxy.first, reco_);
        lxy = hgdc.assignCell(xy.first, xy.second, i + loff, 0, reco_);
        flg = (kxy == lxy) ? " " : " ***** Error *****";
        edm::LogVerbatim("HGCalGeom") << "Input: (" << localx << "," << localy << "," << i + loff << ", " << subsec
                                      << "), assignCell o/p (" << kxy.first << ", " << kxy.second
                                      << ") locateCell o/p (" << xy.first << ", " << xy.second << "),"
                                      << " final (" << lxy.first << ", " << lxy.second << ")" << flg;
        kxy = hgdc.assignCell(-localx, -localy, i + loff, subsec, reco_);
        xy = hgdc.locateCell(kxy.second, i + loff, kxy.first, reco_);
        lxy = hgdc.assignCell(xy.first, xy.second, i + loff, 0, reco_);
        flg = (kxy == lxy) ? " " : " ***** Error *****";
        edm::LogVerbatim("HGCalGeom") << "Input: (" << -localx << "," << -localy << "," << i + loff << ", " << subsec
                                      << "), assignCell o/p (" << kxy.first << ", " << kxy.second
                                      << ") locateCell o/p (" << xy.first << ", " << xy.second << "), final ("
                                      << lxy.first << ", " << lxy.second << ")" << flg;
      } else if (detType_ == 0) {
        std::array<int, 3> kxy, lxy;
        kxy = hgdc.assignCellTrap(localx, localy, zpos, i + loff, reco_);
        xy = hgdc.locateCellTrap(zside, i + loff, kxy[0], kxy[1], reco_, false);
        lxy = hgdc.assignCellTrap(xy.first, xy.second, zpos, i + loff, reco_);
        flg = (kxy == lxy) ? " " : " ***** Error *****";
        edm::LogVerbatim("HGCalGeom") << "Input: (" << localx << "," << localy << "," << zpos << ", " << i + loff
                                      << "), assignCell o/p (" << kxy[0] << ":" << kxy[1] << ":" << kxy[2]
                                      << ") locateCell o/p (" << xy.first << ", " << xy.second << "), final (" << lxy[0]
                                      << ":" << lxy[1] << ":" << lxy[2] << ") Dist "
                                      << hgdc.distFromEdgeTrap(scl * localx, scl * localy, scl * zpos) << " " << flg;
        kxy = hgdc.assignCellTrap(-localx, -localy, zpos, i + loff, reco_);
        xy = hgdc.locateCellTrap(zside, i + loff, kxy[0], kxy[1], reco_, false);
        lxy = hgdc.assignCellTrap(xy.first, xy.second, zpos, i + loff, reco_);
        flg = (kxy == lxy) ? " " : " ***** Error *****";
        edm::LogVerbatim("HGCalGeom") << "Input: (" << -localx << "," << -localy << "," << zpos << ", " << i + loff
                                      << "), assignCell o/p (" << kxy[0] << ":" << kxy[1] << ":" << kxy[2]
                                      << ") locateCell o/p (" << xy.first << ", " << xy.second << "), final (" << lxy[0]
                                      << ":" << lxy[1] << ":" << lxy[2] << ") Dist "
                                      << hgdc.distFromEdgeTrap(scl * localx, scl * localy, scl * zpos) << " " << flg;
      } else {
        std::array<int, 5> kxy, lxy;
        kxy = hgdc.assignCellHex(localx, localy, zside, i + loff, reco_, false, false);
        xy = hgdc.locateCell(zside, i + loff, kxy[0], kxy[1], kxy[3], kxy[4], reco_, true, false, false, false);
        lxy = hgdc.assignCellHex(xy.first, xy.second, zside, i + loff, reco_, false, false);
        flg = (kxy == lxy) ? " " : " ***** Error *****";
        edm::LogVerbatim("HGCalGeom") << "Input: (" << localx << "," << localy << ", " << i + loff
                                      << "), assignCell o/p (" << kxy[0] << ":" << kxy[1] << ":" << kxy[2] << ":"
                                      << kxy[3] << ":" << kxy[4] << ") locateCell o/p (" << xy.first << ", "
                                      << xy.second << "), final (" << lxy[0] << ":" << lxy[1] << ":" << lxy[2] << ":"
                                      << lxy[3] << ":" << lxy[4] << ") Dist "
                                      << hgdc.distFromEdgeHex(scl * localx, scl * localy, scl * zpos) << " " << flg;
        kxy = hgdc.assignCellHex(-localx, -localy, zside, i + loff, reco_, false, false);
        xy = hgdc.locateCell(zside, i + loff, kxy[0], kxy[1], kxy[3], kxy[4], reco_, true, false, false, false);
        lxy = hgdc.assignCellHex(xy.first, xy.second, zside, i + loff, reco_, false, false);
        flg = (kxy == lxy) ? " " : " ***** Error *****";
        edm::LogVerbatim("HGCalGeom") << "Input: (" << -localx << "," << -localy << ", " << i + loff
                                      << "), assignCell o/p (" << kxy[0] << ":" << kxy[1] << ":" << kxy[2] << ":"
                                      << kxy[3] << ":" << kxy[4] << ") locateCell o/p (" << xy.first << ", "
                                      << xy.second << "), final (" << lxy[0] << ":" << lxy[1] << ":" << lxy[2] << ":"
                                      << lxy[3] << ":" << lxy[4] << ") Dist "
                                      << hgdc.distFromEdgeHex(scl * localx, scl * localy, scl * zpos) << " " << flg;
      }
      if (k == 0 && i == 0 && detType_ == 1) {
        std::vector<int> ncells = hgdc.numberCells(i + 1, reco_);
        edm::LogVerbatim("HGCalGeom") << "Layer " << i + 1 << " with " << ncells.size() << " rows";
        int ntot(0);
        for (unsigned int k = 0; k < ncells.size(); ++k) {
          ntot += ncells[k];
          edm::LogVerbatim("HGCalGeom") << "Row " << k << " with " << ncells[k] << " cells";
        }
        edm::LogVerbatim("HGCalGeom") << "Total Cells " << ntot << ":" << hgdc.maxCells(i + 1, reco_);
      }
      i += increment_;
    }
  }

  // Test the range variables
  if (hgdc.getParameter()->detectorType_ > 0) {
    unsigned int kk(0);
    for (auto const& zz : hgdc.getParameter()->zLayerHex_) {
      std::pair<double, double> rr = hgdc.rangeR(zz, true);
      edm::LogVerbatim("HGCalGeom") << "[" << kk << "]\t z = " << zz << "\t rMin = " << rr.first
                                    << "\t rMax = " << rr.second;
      ++kk;
    }
  }

  // For scintillators
  if (hgdc.tileTrapezoid()) {
    unsigned int kk(0);
    for (auto const& lay : hgdc.getParameter()->layer_) {
      auto rRange = hgdc.getREtaRange(lay);
      edm::LogVerbatim("HGCalGeom") << "[" << kk << "] Layer " << lay << " R/Eta " << rRange.first << ":"
                                    << rRange.second << " nPhi " << hgdc.getPhiBins(lay);
      ++kk;
    }
  }
}

// define this as a plug-in
DEFINE_FWK_MODULE(HGCalNumberingTester);