00001 #ifndef _CSCFAKEDBNOISEMATRIX_H 00002 #define _CSCFAKEDBNOISEMATRIX_H 00003 00004 #include <memory> 00005 #include "FWCore/Framework/interface/SourceFactory.h" 00006 #include "FWCore/Framework/interface/Frameworkfwd.h" 00007 #include "FWCore/Framework/interface/ESProducer.h" 00008 #include "FWCore/Framework/interface/Event.h" 00009 #include "FWCore/Framework/interface/MakerMacros.h" 00010 #include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h" 00011 #include "FWCore/Framework/interface/ESHandle.h" 00012 #include "FWCore/Framework/interface/EventSetup.h" 00013 00014 #include "FWCore/ParameterSet/interface/ParameterSet.h" 00015 00016 #include "CondFormats/CSCObjects/interface/CSCDBNoiseMatrix.h" 00017 #include "CondFormats/DataRecord/interface/CSCDBNoiseMatrixRcd.h" 00018 #include <DataFormats/MuonDetId/interface/CSCDetId.h> 00019 #include <boost/shared_ptr.hpp> 00020 00021 class CSCFakeDBNoiseMatrix: public edm::ESProducer, public edm::EventSetupRecordIntervalFinder { 00022 public: 00023 CSCFakeDBNoiseMatrix(const edm::ParameterSet&); 00024 ~CSCFakeDBNoiseMatrix(); 00025 00026 inline static CSCDBNoiseMatrix * prefillDBNoiseMatrix(); 00027 00028 typedef boost::shared_ptr<CSCDBNoiseMatrix> Pointer; 00029 00030 Pointer produceDBNoiseMatrix(const CSCDBNoiseMatrixRcd&); 00031 00032 private: 00033 // ----------member data --------------------------- 00034 void setIntervalFor(const edm::eventsetup::EventSetupRecordKey &, const edm::IOVSyncValue&, edm::ValidityInterval & ); 00035 Pointer cndbNoiseMatrix ; 00036 00037 }; 00038 00039 #include<fstream> 00040 #include<vector> 00041 #include<iostream> 00042 00043 // to workaround plugin library 00044 inline CSCDBNoiseMatrix * CSCFakeDBNoiseMatrix::prefillDBNoiseMatrix() 00045 { 00046 int seed; 00047 long int M; 00048 float mean,min, minchi; 00049 const int MAX_SIZE = 217728; //or 252288 for ME4/2 chambers 00050 const int FACTOR=1000; 00051 00052 CSCDBNoiseMatrix * cndbmatrix = new CSCDBNoiseMatrix(); 00053 cndbmatrix->matrix.resize(MAX_SIZE); 00054 00055 seed = 10000; 00056 srand(seed); 00057 mean=6.8, min=-10.0, minchi=1.0, M=1000; 00058 cndbmatrix->factor_noise = int (FACTOR); 00059 00060 for(int i=0; i<MAX_SIZE;i++){ 00061 cndbmatrix->matrix[i].elem33 = (short int) (10.0*FACTOR+0.5); 00062 cndbmatrix->matrix[i].elem34 = (short int) (4.0*FACTOR+0.5); 00063 cndbmatrix->matrix[i].elem44 = (short int) (10.0*FACTOR+0.5); 00064 cndbmatrix->matrix[i].elem35 = (short int) (3.0*FACTOR+0.5); 00065 cndbmatrix->matrix[i].elem45 = (short int) (8.0*FACTOR+0.5); 00066 cndbmatrix->matrix[i].elem55 = (short int) (10.0*FACTOR+0.5); 00067 cndbmatrix->matrix[i].elem46 = (short int) (2.0*FACTOR+0.5); 00068 cndbmatrix->matrix[i].elem56 = (short int) (5.0*FACTOR+0.5); 00069 cndbmatrix->matrix[i].elem66 = (short int) (10.0*FACTOR+0.5); 00070 cndbmatrix->matrix[i].elem57 = (short int) (3.0*FACTOR+0.5); 00071 cndbmatrix->matrix[i].elem67 = (short int) (4.0*FACTOR+0.5); 00072 cndbmatrix->matrix[i].elem77 = (short int) (10.0*FACTOR+0.5); 00073 } 00074 return cndbmatrix; 00075 } 00076 00077 #endif 00078