Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-03-17 11:21:39

0001 /****************************************************************************
0002  *
0003  * This is a part of PPS offline software.
0004  * Authors:
0005  *   Laurent Forthomme (laurent.forthomme@cern.ch)
0006  *   Nicola Minafra (nicola.minafra@cern.ch)
0007  *
0008  ****************************************************************************/
0009 
0010 #include <memory>
0011 
0012 #include "FWCore/Framework/interface/Frameworkfwd.h"
0013 #include "FWCore/Framework/interface/stream/EDProducer.h"
0014 #include "FWCore/Framework/interface/Event.h"
0015 #include "FWCore/Framework/interface/MakerMacros.h"
0016 #include "FWCore/Framework/interface/ESHandle.h"
0017 #include "FWCore/Framework/interface/ESWatcher.h"
0018 
0019 #include "FWCore/ParameterSet/interface/ParameterSet.h"
0020 #include "FWCore/Utilities/interface/StreamID.h"
0021 #include "FWCore/Utilities/interface/ESGetToken.h"
0022 
0023 #include "DataFormats/Common/interface/DetSetVector.h"
0024 #include "DataFormats/Common/interface/DetSet.h"
0025 
0026 #include "DataFormats/CTPPSDetId/interface/CTPPSDiamondDetId.h"
0027 #include "DataFormats/CTPPSDigi/interface/CTPPSDiamondDigi.h"
0028 #include "DataFormats/CTPPSReco/interface/CTPPSDiamondRecHit.h"
0029 
0030 #include "RecoPPS/Local/interface/CTPPSDiamondRecHitProducerAlgorithm.h"
0031 
0032 #include "Geometry/Records/interface/VeryForwardRealGeometryRecord.h"
0033 #include "CondFormats/DataRecord/interface/PPSTimingCalibrationRcd.h"
0034 #include "CondFormats/DataRecord/interface/PPSTimingCalibrationLUTRcd.h"
0035 
0036 class CTPPSDiamondRecHitProducer : public edm::stream::EDProducer<> {
0037 public:
0038   explicit CTPPSDiamondRecHitProducer(const edm::ParameterSet&);
0039 
0040   static void fillDescriptions(edm::ConfigurationDescriptions&);
0041 
0042 private:
0043   void produce(edm::Event&, const edm::EventSetup&) override;
0044 
0045   edm::EDGetTokenT<edm::DetSetVector<CTPPSDiamondDigi> > digiToken_;
0046   edm::ESGetToken<PPSTimingCalibration, PPSTimingCalibrationRcd> timingCalibrationToken_;
0047   edm::ESGetToken<PPSTimingCalibrationLUT, PPSTimingCalibrationLUTRcd> timingCalibrationLUTToken_;
0048   edm::ESGetToken<CTPPSGeometry, VeryForwardRealGeometryRecord> geometryToken_;
0049 
0050   /// A watcher to detect timing calibration changes.
0051   edm::ESWatcher<PPSTimingCalibrationRcd> calibWatcher_;
0052 
0053   bool applyCalib_;
0054   CTPPSDiamondRecHitProducerAlgorithm algo_;
0055 };
0056 
0057 CTPPSDiamondRecHitProducer::CTPPSDiamondRecHitProducer(const edm::ParameterSet& iConfig)
0058     : digiToken_(consumes<edm::DetSetVector<CTPPSDiamondDigi> >(iConfig.getParameter<edm::InputTag>("digiTag"))),
0059       geometryToken_(esConsumes<CTPPSGeometry, VeryForwardRealGeometryRecord>()),
0060       applyCalib_(iConfig.getParameter<bool>("applyCalibration")),
0061       algo_(iConfig) {
0062   if (applyCalib_) {
0063     timingCalibrationToken_ = esConsumes<PPSTimingCalibration, PPSTimingCalibrationRcd>(
0064         edm::ESInputTag(iConfig.getParameter<std::string>("timingCalibrationTag")));
0065     timingCalibrationLUTToken_ = esConsumes<PPSTimingCalibrationLUT, PPSTimingCalibrationLUTRcd>();
0066   }
0067   produces<edm::DetSetVector<CTPPSDiamondRecHit> >();
0068 }
0069 
0070 void CTPPSDiamondRecHitProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
0071   auto pOut = std::make_unique<edm::DetSetVector<CTPPSDiamondRecHit> >();
0072 
0073   // get the digi collection
0074   const auto& digis = iEvent.get(digiToken_);
0075 
0076   if (!digis.empty()) {
0077     if (applyCalib_ && calibWatcher_.check(iSetup))
0078       algo_.setCalibration(iSetup.getData(timingCalibrationToken_), iSetup.getData(timingCalibrationLUTToken_));
0079 
0080     // produce the rechits collection
0081     algo_.build(iSetup.getData(geometryToken_), digis, *pOut);
0082   }
0083 
0084   iEvent.put(std::move(pOut));
0085 }
0086 
0087 void CTPPSDiamondRecHitProducer::fillDescriptions(edm::ConfigurationDescriptions& descr) {
0088   edm::ParameterSetDescription desc;
0089 
0090   desc.add<edm::InputTag>("digiTag", edm::InputTag("ctppsDiamondRawToDigi", "TimingDiamond"))
0091       ->setComment("input digis collection to retrieve");
0092   desc.add<std::string>("timingCalibrationTag", "GlobalTag:PPSDiamondTimingCalibration")
0093       ->setComment("input tag for timing calibrations retrieval");
0094   desc.add<double>("timeSliceNs", 25.0 / 1024.0)
0095       ->setComment("conversion constant between HPTDC timing bin size and nanoseconds");
0096   desc.add<bool>("applyCalibration", true)->setComment("switch on/off the timing calibration");
0097 
0098   descr.add("ctppsDiamondRecHits", desc);
0099 }
0100 
0101 DEFINE_FWK_MODULE(CTPPSDiamondRecHitProducer);