|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
|
#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
|