File indexing completed on 2024-10-17 04:57:59
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020 #include <memory>
0021
0022
0023 #include "CalibTracker/SiPixelESProducers/interface/SiPixelDetInfoFileReader.h"
0024 #include "CondFormats/DataRecord/interface/SiPixelGainCalibrationOfflineRcd.h"
0025 #include "CondFormats/SiPixelObjects/interface/SiPixelGainCalibrationOffline.h"
0026 #include "FWCore/Framework/interface/ESProducer.h"
0027 #include "FWCore/Framework/interface/EventSetupRecordIntervalFinder.h"
0028 #include "FWCore/Framework/interface/MakerMacros.h"
0029 #include "FWCore/Framework/interface/ModuleFactory.h"
0030 #include "FWCore/Framework/interface/SourceFactory.h"
0031 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0032 #include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
0033 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0034 #include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
0035 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
0036 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
0037
0038
0039
0040
0041
0042 class SiPixelFakeGainOfflineESSource : public edm::ESProducer, public edm::EventSetupRecordIntervalFinder {
0043 public:
0044 SiPixelFakeGainOfflineESSource(const edm::ParameterSet&);
0045 ~SiPixelFakeGainOfflineESSource() override = default;
0046
0047 static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
0048
0049 virtual std::unique_ptr<SiPixelGainCalibrationOffline> produce(const SiPixelGainCalibrationOfflineRcd&);
0050
0051 protected:
0052 void setIntervalFor(const edm::eventsetup::EventSetupRecordKey&,
0053 const edm::IOVSyncValue&,
0054 edm::ValidityInterval&) override;
0055
0056 private:
0057 edm::FileInPath fp_;
0058 };
0059
0060
0061
0062
0063 SiPixelFakeGainOfflineESSource::SiPixelFakeGainOfflineESSource(const edm::ParameterSet& conf_)
0064 : fp_(conf_.getParameter<edm::FileInPath>("file")) {
0065 edm::LogInfo("SiPixelFakeGainOfflineESSource::SiPixelFakeGainOfflineESSource");
0066
0067
0068 setWhatProduced(this);
0069 findingRecord<SiPixelGainCalibrationOfflineRcd>();
0070 }
0071
0072 std::unique_ptr<SiPixelGainCalibrationOffline> SiPixelFakeGainOfflineESSource::produce(
0073 const SiPixelGainCalibrationOfflineRcd&) {
0074 using namespace edm::es;
0075 SiPixelGainCalibrationOffline* obj = new SiPixelGainCalibrationOffline(25., 30., 2., 3.);
0076 SiPixelDetInfoFileReader reader(fp_.fullPath());
0077 const std::vector<uint32_t>& DetIds = reader.getAllDetIds();
0078
0079
0080 for (std::vector<uint32_t>::const_iterator detit = DetIds.begin(); detit != DetIds.end(); detit++) {
0081 std::vector<char> theSiPixelGainCalibrationOffline;
0082 const std::pair<int, int>& detUnitDimensions = reader.getDetUnitDimensions(*detit);
0083
0084
0085 for (int i = 0; i < detUnitDimensions.first; i++) {
0086 float totalGain = 0.0;
0087 float totalEntries = 0.0;
0088 for (int j = 0; j < detUnitDimensions.second; j++) {
0089 totalGain += 2.8;
0090 float ped = 28.2;
0091 totalEntries += 1.0;
0092 obj->setDataPedestal(ped, theSiPixelGainCalibrationOffline);
0093 if ((j + 1) % 80 == 0)
0094 {
0095 float gain = totalGain / totalEntries;
0096 obj->setDataGain(gain, 80, theSiPixelGainCalibrationOffline);
0097 totalGain = 0;
0098 totalEntries = 0.0;
0099 }
0100 }
0101 }
0102
0103
0104
0105 SiPixelGainCalibrationOffline::Range range(theSiPixelGainCalibrationOffline.begin(),
0106 theSiPixelGainCalibrationOffline.end());
0107
0108 if (!obj->put(*detit, range, detUnitDimensions.first))
0109 edm::LogError("SiPixelFakeGainOfflineESSource")
0110 << "[SiPixelFakeGainOfflineESSource::produce] detid already exists" << std::endl;
0111 }
0112
0113
0114 return std::unique_ptr<SiPixelGainCalibrationOffline>(obj);
0115 }
0116
0117 void SiPixelFakeGainOfflineESSource::setIntervalFor(const edm::eventsetup::EventSetupRecordKey&,
0118 const edm::IOVSyncValue& iosv,
0119 edm::ValidityInterval& oValidity) {
0120 edm::ValidityInterval infinity(iosv.beginOfTime(), iosv.endOfTime());
0121 oValidity = infinity;
0122 }
0123
0124 void SiPixelFakeGainOfflineESSource::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
0125 edm::ParameterSetDescription desc;
0126 desc.add<edm::FileInPath>("file", edm::FileInPath("CalibTracker/SiPixelESProducers/data/PixelSkimmedGeometry.txt"));
0127 descriptions.addWithDefaultLabel(desc);
0128 }
0129
0130 DEFINE_FWK_EVENTSETUP_SOURCE(SiPixelFakeGainOfflineESSource);