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