CSCFakeDBNoiseMatrix.h

Go to the documentation of this file.

#ifndef _CSCFAKEDBNOISEMATRIX_H
#define _CSCFAKEDBNOISEMATRIX_H

#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/SourceFactory.h"
#include <memory>

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "CondFormats/CSCObjects/interface/CSCDBNoiseMatrix.h"
#include "CondFormats/DataRecord/interface/CSCDBNoiseMatrixRcd.h"
#include <DataFormats/MuonDetId/interface/CSCDetId.h>

class CSCFakeDBNoiseMatrix : public edm::ESProducer, public edm::EventSetupRecordIntervalFinder {
public:
  CSCFakeDBNoiseMatrix(const edm::ParameterSet &);
  ~CSCFakeDBNoiseMatrix() override;

  inline static CSCDBNoiseMatrix *prefillDBNoiseMatrix();

  typedef std::unique_ptr<CSCDBNoiseMatrix> Pointer;

  Pointer produceDBNoiseMatrix(const CSCDBNoiseMatrixRcd &);

private:
  // ----------member data ---------------------------
  void setIntervalFor(const edm::eventsetup::EventSetupRecordKey &,
                      const edm::IOVSyncValue &,
                      edm::ValidityInterval &) override;
  Pointer cndbNoiseMatrix;
};

#include <fstream>
#include <iostream>
#include <vector>

// to workaround plugin library
inline CSCDBNoiseMatrix *CSCFakeDBNoiseMatrix::prefillDBNoiseMatrix() {
  int seed;
  const int MAX_SIZE = 252288;  // or 252288 for ME4/2 chambers
  const int FACTOR = 1000;

  CSCDBNoiseMatrix *cndbmatrix = new CSCDBNoiseMatrix();
  cndbmatrix->matrix.resize(MAX_SIZE);

  seed = 10000;
  srand(seed);
  cndbmatrix->factor_noise = int(FACTOR);

  for (int i = 0; i < MAX_SIZE; i++) {
    cndbmatrix->matrix[i].elem33 = (short int)(10.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem34 = (short int)(4.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem35 = (short int)(3.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem44 = (short int)(10.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem45 = (short int)(8.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem46 = (short int)(2.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem55 = (short int)(10.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem56 = (short int)(5.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem57 = (short int)(3.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem66 = (short int)(10.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem67 = (short int)(4.0 * FACTOR + 0.5);
    cndbmatrix->matrix[i].elem77 = (short int)(10.0 * FACTOR + 0.5);
  }
  return cndbmatrix;
}

#endif