RooKeysPdf.h

Go to the documentation of this file.

/*****************************************************************************
 * Project: RooFit                                                           *
 * Package: RooFitModels                                                     *
 *    File: $Id: RooKeysPdf.h,v 1.10 2007/05/11 09:13:07 verkerke Exp $
 * Authors:                                                                  *
 *   GR, Gerhard Raven,   UC San Diego,        raven@slac.stanford.edu       *
 *   DK, David Kirkby,    UC Irvine,         dkirkby@uci.edu                 *
 *   WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu       *
 *                                                                           *
 * Copyright (c) 2000-2005, Regents of the University of California          *
 *                          and Stanford University. All rights reserved.    *
 *                                                                           *
 * Redistribution and use in source and binary forms,                        *
 * with or without modification, are permitted according to the terms        *
 * listed in LICENSE (http://roofit.sourceforge.net/license.txt)             *
 *****************************************************************************/
#ifndef ROO_KEYS
#define ROO_KEYS
 
#include "RooAbsPdf.h"
#include "RooRealProxy.h"
 
class RooRealVar;
 
class RooKeysPdf : public RooAbsPdf {
public:
  /// Boundary correction obtained by reflecting the data across the lower
  /// and/or upper edge of the observable range. Symmetric mirroring adds the
  /// reflected events (density flat at the boundary), asymmetric mirroring
  /// subtracts them (density vanishing at the boundary).
  enum Mirror {
     NoMirror,            ///< No boundary correction
     MirrorLeft,          ///< Symmetric mirror at the lower edge
     MirrorRight,         ///< Symmetric mirror at the upper edge
     MirrorBoth,          ///< Symmetric mirror at both edges
     MirrorAsymLeft,      ///< Asymmetric mirror at the lower edge
     MirrorAsymLeftRight, ///< Asymmetric mirror at the lower edge, symmetric mirror at the upper edge
     MirrorAsymRight,     ///< Asymmetric mirror at the upper edge
     MirrorLeftAsymRight, ///< Symmetric mirror at the lower edge, asymmetric mirror at the upper edge
     MirrorAsymBoth       ///< Asymmetric mirror at both edges
  };
  RooKeysPdf() ;
  RooKeysPdf(const char *name, const char *title,
             RooAbsReal& xpdf, RooDataSet& data, Mirror mirror= NoMirror,
        double rho=1);
  RooKeysPdf(const char *name, const char *title,
             RooAbsReal& x, RooRealVar& xdata, RooDataSet& data, Mirror mirror= NoMirror,
        double rho=1);
  RooKeysPdf(const RooKeysPdf& other, const char* name=nullptr);
  TObject* clone(const char* newname=nullptr) const override {return new RooKeysPdf(*this,newname); }
  ~RooKeysPdf() override;
 
  Int_t getAnalyticalIntegral(RooArgSet& allVars, RooArgSet& analVars,
     const char* rangeName = nullptr) const override;
  double analyticalIntegral(Int_t code, const char* rangeName = nullptr) const override;
  Int_t getMaxVal(const RooArgSet& vars) const override;
  double maxVal(Int_t code) const override;
 
  void LoadDataSet( RooDataSet& data);
 
protected:
 
  RooRealProxy _x ;
  double evaluate() const override;
 
private:
  // how far you have to go out in a Gaussian until it is smaller than the
  // machine precision
  static const double _nSigma; ///<!
 
  Int_t _nEvents = 0;
  double *_dataPts = nullptr;  //[_nEvents]
  double *_dataWgts = nullptr; //[_nEvents]
  double *_weights = nullptr;  //[_nEvents]
  double _sumWgt = 0.0;
 
  constexpr static int _nPoints{1000};
  double _lookupTable[_nPoints+1];
 
  double g(double x,double sigma) const;
 
  bool _mirrorLeft = false;
  bool _mirrorRight = false;
  bool _asymLeft = false;
  bool _asymRight = false;
 
  // cached info on variable
  Char_t _varName[128];
  double _lo, _hi, _binWidth;
  double _rho;
 
  ClassDefOverride(RooKeysPdf,2) // One-dimensional non-parametric kernel estimation p.d.f.
};
 
#endif