/data/refman/pasoursint/CMSSW_6_1_2_SLHC2_patch1/src/RecoVertex/BeamSpotProducer/interface/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() {} 
  // 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 {return errorDef_;} 
  // -2lnL value based on vector of parameters
  double operator() (const std::vector<double>&) const; 
  // 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