SiPixelESProducers/plugins/SiPixelFakeGainOfflineESSource.cc

0001 // -*- C++ -*-
0002 //
0003 // Package:    CalibTracker/SiPixelGainESProducers
0004 // Class:      SiPixelFakeGainOfflineESSource
0005 //
0006 /**\class SiPixelFakeGainOfflineESSource SiPixelFakeGainOfflineESSource.cc CalibTracker/SiPixelGainESProducers/plugins/SiPixelFakeGainOfflineESSource.cc
0007
0008  Description: <one line class summary>
0009
0010  Implementation:
0011      <Notes on implementation>
0012 */
0013 //
0014 // Original Author:  Vincenzo Chiochia
0015 //         Created:  Tue 8 12:31:25 CEST 2007
0016 //
0017 //
0018
0019 // system include files
0020 #include <memory>
0021
0022 // user include files
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 // class decleration
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 // constructors and destructor
0062 //
0063 SiPixelFakeGainOfflineESSource::SiPixelFakeGainOfflineESSource(const edm::ParameterSet& conf_)
0064     : fp_(conf_.getParameter<edm::FileInPath>("file")) {
0065   edm::LogInfo("SiPixelFakeGainOfflineESSource::SiPixelFakeGainOfflineESSource");
0066   //the following line is needed to tell the framework what
0067   // data is being produced
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   // Loop over detectors
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     // Loop over columns and rows
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)  //compute the gain average after each ROC
0094         {
0095           float gain = totalGain / totalEntries;
0096           obj->setDataGain(gain, 80, theSiPixelGainCalibrationOffline);
0097           totalGain = 0;
0098           totalEntries = 0.0;
0099         }
0100       }
0101     }
0102
0103     //std::cout << "detid " << (*detit) << std::endl;
0104
0105     SiPixelGainCalibrationOffline::Range range(theSiPixelGainCalibrationOffline.begin(),
0106                                                theSiPixelGainCalibrationOffline.end());
0107     // the 80 in the line below represents the number of columns averaged over.
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);