SimMuon/CSCDigitizer/test/CSCNoiseMatrixTest.cc - master - githubmirror/cmssw

SimMuon/CSCDigitizer/test/CSCNoiseMatrixTest.cc
67 строк3 KB
Matti Kortelainen
Move RandomNumberGenerator service base class to FWCore/AbstractServices
28 мар 2025, 00:56
28 мар 2025, 00:56c4a96e9
#include "CalibMuon/CSCCalibration/interface/CSCConditions.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "Geometry/CSCGeometry/interface/CSCGeometry.h"
#include "Geometry/CSCGeometry/interface/CSCLayerGeometry.h"
#include "Geometry/Records/interface/MuonGeometryRecord.h"
#include "SimMuon/CSCDigitizer/src/CSCDbStripConditions.h"

#include "FWCore/AbstractServices/interface/RandomNumberGenerator.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/EDMException.h"

class CSCNoiseMatrixTest : public edm::one::EDAnalyzer<> {
public:
  CSCNoiseMatrixTest(const edm::ParameterSet &pset)
      : theDbConditions(pset, consumesCollector()), tokGeo_{esConsumes(edm::ESInputTag("", "idealForDigi"))} {
    edm::Service<edm::RandomNumberGenerator> rng;
    if (!rng.isAvailable()) {
      throw cms::Exception("Configuration") << "CSCNoiseMatrixTest requires the RandomNumberGeneratorService\n"
                                               "which is not present in the configuration file.  You must add "
                                               "the service\n"
                                               "in the configuration file or remove the modules that require it.";
    }
  }

  ~CSCNoiseMatrixTest() override {}
  // virtual void beginRun(edm::Run const& run, edm::EventSetup const&
  // eventSetup)
  void analyze(const edm::Event &e, const edm::EventSetup &eventSetup) override {
    edm::Service<edm::RandomNumberGenerator> rng;
    CLHEP::HepRandomEngine *engine = &rng->getEngine(e.streamID());

    theDbConditions.initializeEvent(eventSetup);
    // fake signal is going to get noisified 200000 times!
    const int nScaBins = 8;
    const float scaBinSize = 50.;
    std::vector<float> binValues(nScaBins, 0.);

    // find the geometry & conditions for this event
    const CSCGeometry *pGeom = &eventSetup.getData(tokGeo_);

    // try making a noisifier and using it
    const CSCGeometry::LayerContainer &layers = pGeom->layers();
    for (auto layerItr = layers.begin(); layerItr != layers.end(); ++layerItr) {
      unsigned nstrips = (**layerItr).geometry()->numberOfStrips();
      for (unsigned istrip = 1; istrip != nstrips; ++istrip) {
        CSCAnalogSignal signal(istrip, scaBinSize, binValues, 0., 0.);
        // theDbConditions.fetchNoisifier((**layerItr).id(),istrip);
        try {
          theDbConditions.noisify((**layerItr).id(), signal, engine);
        } catch (cms::Exception &e) {
          std::cerr << "Bad Noise Matrix for " << (**layerItr).id() << "strip " << istrip << "\n";
          std::cerr << e.what() << std::endl;
        }
      }
    }
  }

  // virtual void analyze(const edm::Event&, const edm::EventSetup&) {}

private:
  CSCDbStripConditions theDbConditions;
  const edm::ESGetToken<CSCGeometry, MuonGeometryRecord> tokGeo_;
};

DEFINE_FWK_MODULE(CSCNoiseMatrixTest);