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