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 RecoCTPPS_TotemRPLocal_TotemTimingRecHitProducerAlgorithm
#define RecoCTPPS_TotemRPLocal_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/CTPPSReadoutObjects/interface/PPSTimingCalibration.h"

#include "RecoCTPPS/TotemRPLocal/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