TotemTimingRecHitProducer.cc

Go to the documentation of this file.

/****************************************************************************
 *
 * This is a part of CTPPS offline software.
 * Authors:
 *   Laurent Forthomme (laurent.forthomme@cern.ch)
 *   Nicola Minafra (nicola.minafra@cern.ch)
 *   Christopher Misan (krzysztof.misan@cern.ch)
 *
 ****************************************************************************/

#include <memory>

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/Framework/interface/ESWatcher.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/StreamID.h"
#include "FWCore/Utilities/interface/ESGetToken.h"

#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/Common/interface/DetSet.h"

#include "DataFormats/CTPPSDetId/interface/TotemTimingDetId.h"
#include "DataFormats/CTPPSDigi/interface/TotemTimingDigi.h"
#include "DataFormats/CTPPSReco/interface/TotemTimingRecHit.h"

#include "RecoPPS/Local/interface/TotemTimingRecHitProducerAlgorithm.h"

#include "Geometry/Records/interface/VeryForwardRealGeometryRecord.h"
#include "Geometry/VeryForwardGeometryBuilder/interface/CTPPSGeometry.h"
#include "CondFormats/DataRecord/interface/PPSTimingCalibrationRcd.h"

class TotemTimingRecHitProducer : public edm::stream::EDProducer<> {
public:
  explicit TotemTimingRecHitProducer(const edm::ParameterSet&);

  static void fillDescriptions(edm::ConfigurationDescriptions&);

private:
  void produce(edm::Event&, const edm::EventSetup&) override;

  edm::EDGetTokenT<edm::DetSetVector<TotemTimingDigi> > digiToken_;
  edm::ESGetToken<PPSTimingCalibration, PPSTimingCalibrationRcd> timingCalibrationToken_;
  edm::ESGetToken<CTPPSGeometry, VeryForwardRealGeometryRecord> geometryToken_;
  TotemTimingRecHitProducerAlgorithm algo_;
  edm::ESWatcher<PPSTimingCalibrationRcd> calibWatcher_;
};

TotemTimingRecHitProducer::TotemTimingRecHitProducer(const edm::ParameterSet& iConfig)
    : digiToken_(consumes<edm::DetSetVector<TotemTimingDigi> >(iConfig.getParameter<edm::InputTag>("digiTag"))),
      timingCalibrationToken_(esConsumes<PPSTimingCalibration, PPSTimingCalibrationRcd>(
          edm::ESInputTag(iConfig.getParameter<std::string>("timingCalibrationTag")))),
      geometryToken_(esConsumes<CTPPSGeometry, VeryForwardRealGeometryRecord>()),
      algo_(iConfig) {
  produces<edm::DetSetVector<TotemTimingRecHit> >();
}

void TotemTimingRecHitProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  std::unique_ptr<edm::DetSetVector<TotemTimingRecHit> > pOut(new edm::DetSetVector<TotemTimingRecHit>);

  // get the digi collection
  edm::Handle<edm::DetSetVector<TotemTimingDigi> > digis;
  iEvent.getByToken(digiToken_, digis);

  // do not retrieve the calibration parameters if no digis were found
  if (!digis->empty()) {
    // check for timing calibration parameters update
    if (calibWatcher_.check(iSetup)) {
      edm::ESHandle<PPSTimingCalibration> hTimingCalib = iSetup.getHandle(timingCalibrationToken_);
      algo_.setCalibration(*hTimingCalib);
    }

    // get the geometry
    edm::ESHandle<CTPPSGeometry> geometry = iSetup.getHandle(geometryToken_);

    // produce the rechits collection
    algo_.build(*geometry, *digis, *pOut);
  }

  iEvent.put(std::move(pOut));
}

void TotemTimingRecHitProducer::fillDescriptions(edm::ConfigurationDescriptions& descr) {
  edm::ParameterSetDescription desc;

  desc.add<edm::InputTag>("digiTag", edm::InputTag("totemTimingRawToDigi", "TotemTiming"))
      ->setComment("input digis collection to retrieve");
  desc.add<std::string>("timingCalibrationTag", "GlobalTag:TotemTimingCalibration")
      ->setComment("input tag for timing calibrations retrieval");
  desc.add<int>("baselinePoints", 8)->setComment("number of points to be used for the baseline");
  desc.add<double>("saturationLimit", 0.1)
      ->setComment(
          "all signals with max > saturationLimit will be considered as saturated for UFSD, similarly with "
          "min<saturationLimit for a Diamond");
  desc.add<double>("cfdFraction", 0.3)->setComment("fraction of the CFD");
  desc.add<int>("smoothingPoints", 20)
      ->setComment("number of points to be used for the smoothing using sinc (lowpass)");
  desc.add<double>("lowPassFrequency", 0.7)
      ->setComment("Frequency (in GHz) for CFD smoothing, 0 for disabling the filter");
  desc.add<double>("hysteresis", 5.e-3)->setComment("hysteresis of the discriminator");
  desc.add<bool>("mergeTimePeaks", true)->setComment("if time peaks schould be merged");
  desc.add<double>("sampicOffset", 1.0)->setComment("offset of the flipped sampic signal");
  desc.add<double>("sampicSamplingPeriodNs", 1. / 7.695)->setComment("sapic sampling period in ns");
  descr.add("totemTimingRecHits", desc);
}

DEFINE_FWK_MODULE(TotemTimingRecHitProducer);