FcnBeamSpotFitPV.h

Go to the documentation of this file.

#ifndef FcnBeamSpotFitPV_h_
#define FcnBeamSpotFitPV_h_

#include "RecoVertex/BeamSpotProducer/interface/BeamSpotFitPVData.h"
#include "Minuit2/FCNBase.h"

#include <vector>

class FcnBeamSpotFitPV : public ROOT::Minuit2::FCNBase {
public:
  // constructor from vertex data
  FcnBeamSpotFitPV(const std::vector<BeamSpotFitPVData>& data);
  ~FcnBeamSpotFitPV() override {}
  // additional vertex selection using limits in x, y, z
  void setLimits(float xmin, float xmax, float ymin, float ymax, float zmin, float zmax);
  // deltaFcn for definition of the uncertainty
  double Up() const override { return errorDef_; }
  // -2lnL value based on vector of parameters
  double operator()(const std::vector<double>&) const override;
  // vertex count used for the fit (after selection)
  unsigned int nrOfVerticesUsed() const;

private:
  const std::vector<BeamSpotFitPVData>& data_;  //< vertex data
  double errorDef_;                             //< error definition for Minuit

  float lowerLimits_[3];  //< lower limits for x,y,z
  float upperLimits_[3];  //< upper limits for x,y,z
};
#endif