RPCFakeCalibration.cc

Go to the documentation of this file.

#include "CalibMuon/RPCCalibration/interface/RPCCalibSetUp.h"
#include "CalibMuon/RPCCalibration/interface/RPCFakeCalibration.h"
#include "CondFormats/DataRecord/interface/RPCStripNoisesRcd.h"
#include "CondFormats/RPCObjects/interface/RPCStripNoises.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/MuonDetId/interface/RPCDetId.h"
#include "FWCore/Framework/interface/EDAnalyzer.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
 
#include <cmath>
#include <fstream>
#include <iostream>
#include <memory>
 
RPCFakeCalibration::RPCFakeCalibration(const edm::ParameterSet &pset) : RPCPerformanceESSource(pset) {
  edm::LogInfo("RPCFakeCalibration::RPCFakeCalibration");
  theRPCCalibSetUp = new RPCCalibSetUp(pset);
}
 
RPCStripNoises *RPCFakeCalibration::makeNoise() {
  RPCStripNoises *obj = new RPCStripNoises();
 
  std::map<int, std::vector<double>>::iterator itc;
  for (itc = (theRPCCalibSetUp->_clsMap).begin(); itc != (theRPCCalibSetUp->_clsMap).end(); ++itc) {
    for (unsigned int n = 0; n < (itc->second).size(); ++n) {
      (obj->v_cls).push_back((itc->second)[n]);
    }
  }
 
  RPCStripNoises::NoiseItem tipoprova;
  for (std::map<uint32_t, std::vector<float>>::iterator it = (theRPCCalibSetUp->_mapDetIdNoise).begin();
       it != (theRPCCalibSetUp->_mapDetIdNoise).end();
       ++it) {
    tipoprova.dpid = it->first;
    tipoprova.time = theRPCCalibSetUp->getTime(it->first);
 
    for (unsigned int k = 0; k < 96; ++k) {
      tipoprova.noise = ((it->second))[k];
      tipoprova.eff = (theRPCCalibSetUp->getEff(it->first))[k];
      (obj->v_noises).push_back(tipoprova);
    }
  }
 
  return obj;
}
 
RPCClusterSize *RPCFakeCalibration::makeCls() {
  RPCClusterSize *obj = new RPCClusterSize();
  RPCClusterSize::ClusterSizeItem rpcClsItem;
 
  for (std::map<uint32_t, std::vector<double>>::iterator it = (theRPCCalibSetUp->_mapDetClsMap).begin();
       it != (theRPCCalibSetUp->_mapDetClsMap).end();
       ++it) {
    rpcClsItem.dpid = it->first;
 
    for (unsigned int k = 0; k < 100; k++) {
      rpcClsItem.clusterSize = (theRPCCalibSetUp->getCls(it->first))[k];
      (obj->v_cls).push_back(rpcClsItem);
    }
  }
  return obj;
}