HcalTestNumberingTest.cc

Go to the documentation of this file.

// -*- C++ -*-
//
// Package:    HcalTestNumberingTester
// Class:      HcalTestNumberingTester
//
//
// Original Author:  Sunanda Banerjee
//         Created:  Mon 2013/12/26
//

// system include files
#include <memory>
#include <iostream>
#include <fstream>

// user include files
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#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 "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "Geometry/Records/interface/IdealGeometryRecord.h"
#include "Geometry/Records/interface/HcalSimNumberingRecord.h"
#include "Geometry/Records/interface/HcalRecNumberingRecord.h"
#include "Geometry/HcalCommonData/interface/HcalDDDSimConstants.h"
#include "Geometry/HcalCommonData/interface/HcalDDDRecConstants.h"
#include "Geometry/HcalCommonData/interface/HcalNumberingFromDDD.h"
#include "Geometry/HcalCommonData/interface/HcalHitRelabeller.h"
#include "SimG4CMS/Calo/interface/HcalNumberingScheme.h"
#include "SimG4CMS/Calo/interface/HcalTestNumberingScheme.h"

class HcalTestNumberingTester : public edm::one::EDAnalyzer<> {
public:
  explicit HcalTestNumberingTester(const edm::ParameterSet&);
  ~HcalTestNumberingTester() override = default;
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

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

private:
  const edm::ESGetToken<HcalDDDSimConstants, HcalSimNumberingRecord> tokSim_;
  const edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> tokReco_;
};

HcalTestNumberingTester::HcalTestNumberingTester(const edm::ParameterSet&)
    : tokSim_(esConsumes<HcalDDDSimConstants, HcalSimNumberingRecord>()),
      tokReco_(esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord>()) {}

void HcalTestNumberingTester::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  descriptions.add("hcalTestNumberingTest", desc);
}

// ------------ method called to produce the data  ------------
void HcalTestNumberingTester::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup) {
  const HcalDDDSimConstants* hcs = &iSetup.getData(tokSim_);
  const HcalDDDRecConstants* hcr = &iSetup.getData(tokReco_);
  HcalNumberingScheme* schme1 = new HcalNumberingScheme();
  HcalNumberingScheme* schme2 = dynamic_cast<HcalNumberingScheme*>(new HcalTestNumberingScheme(false));

  for (int type = 0; type < 2; ++type) {
    HcalSubdetector sub = (type == 0) ? HcalBarrel : HcalEndcap;
    for (int zs = 0; zs < 2; ++zs) {
      int zside = 2 * zs - 1;
      std::pair<int, int> etas = hcr->getEtaRange(type);
      for (int eta = etas.first; eta <= etas.second; ++eta) {
        std::vector<std::pair<int, double> > phis = hcr->getPhis(sub, eta);
        for (unsigned int k = 0; k < phis.size(); ++k) {
          int phi = phis[k].first;
          int lmin = (type == 1 && eta == 16) ? 8 : 1;
          int lmax = (type == 1) ? 19 : ((eta == 16) ? 7 : 17);
          for (int lay = lmin; lay <= lmax; ++lay) {
            std::pair<int, int> etd = hcs->getEtaDepth(sub, eta, phi, zside, 0, lay);
            HcalNumberingFromDDD::HcalID tmp(sub, zs, etd.second, etd.first, phi, phi, lay);
            uint32_t id1 = schme1->getUnitID(tmp);
            uint32_t id2 = schme2->getUnitID(tmp);
            DetId id0 = HcalHitRelabeller::relabel(id2, hcr);
            std::string ok = (id1 != id0.rawId()) ? " *** ERROR ***" : " ";
            edm::LogVerbatim("HcalSim") << "I/P " << sub << ":" << zside * eta << ":" << phi << ":" << lay << " Normal "
                                        << std::hex << id1 << std::dec << " " << HcalDetId(id1) << " Test " << std::hex
                                        << id2 << std::dec << " " << HcalDetId(id0) << ok;
          }
        }
      }
    }
  }
}

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