TotemTimingRecHitProducerAlgorithm.h

Go to the documentation of this file.

/****************************************************************************
 *
 * This is a part of CTPPS offline software.
 * Authors:
 *   Laurent Forthomme (laurent.forthomme@cern.ch)
 *   Nicola Minafra
 *
 ****************************************************************************/
 
#ifndef RecoPPS_Local_TotemTimingRecHitProducerAlgorithm
#define RecoPPS_Local_TotemTimingRecHitProducerAlgorithm
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
 
#include "DataFormats/Common/interface/DetSetVector.h"
 
#include "DataFormats/CTPPSDetId/interface/TotemTimingDetId.h"
#include "DataFormats/CTPPSDigi/interface/TotemTimingDigi.h"
#include "DataFormats/CTPPSReco/interface/TotemTimingRecHit.h"
 
#include "Geometry/VeryForwardGeometryBuilder/interface/CTPPSGeometry.h"
#include "CondFormats/PPSObjects/interface/PPSTimingCalibration.h"
 
#include "RecoPPS/Local/interface/TotemTimingConversions.h"
 
#include <memory>
 
class TotemTimingRecHitProducerAlgorithm {
public:
  TotemTimingRecHitProducerAlgorithm(const edm::ParameterSet& conf);
 
  void setCalibration(const PPSTimingCalibration&);
  void build(const CTPPSGeometry&, const edm::DetSetVector<TotemTimingDigi>&, edm::DetSetVector<TotemTimingRecHit>&);
 
private:
  struct RegressionResults {
    float m, q, rms;
    RegressionResults() : m(0.), q(0.), rms(0.) {}
  };
 
  RegressionResults simplifiedLinearRegression(const std::vector<float>& time,
                                               const std::vector<float>& data,
                                               const unsigned int start_at,
                                               const unsigned int points) const;
 
  int fastDiscriminator(const std::vector<float>& data, float threshold) const;
 
  float constantFractionDiscriminator(const std::vector<float>& time, const std::vector<float>& data);
 
  static constexpr float SINC_COEFFICIENT = M_PI * 2 / 7.8;
 
  std::unique_ptr<TotemTimingConversions> sampicConversions_;
 
  bool mergeTimePeaks_;
  int baselinePoints_;
  double saturationLimit_;
  double cfdFraction_;
  int smoothingPoints_;
  double lowPassFrequency_;
  double hysteresis_;
  TotemTimingRecHit::TimingAlgorithm mode_;
};
 
#endif