#ifndef ROO_GAUSS_MODEL
#define ROO_GAUSS_MODEL
#include "RooResolutionModel.h"
#include "RooRealProxy.h"
#include "RooComplex.h"
#include "RooMath.h"
class RooGaussModel : public RooResolutionModel {
public:
enum RooGaussBasis { noBasis=0, expBasisMinus= 1, expBasisSum= 2, expBasisPlus= 3,
sinBasisMinus=11, sinBasisSum=12, sinBasisPlus=13,
cosBasisMinus=21, cosBasisSum=22, cosBasisPlus=23,
linBasisPlus=33,
quadBasisPlus=43,
coshBasisMinus=51,coshBasisSum=52,coshBasisPlus=53,
sinhBasisMinus=61,sinhBasisSum=62,sinhBasisPlus=63};
enum BasisType { none=0, expBasis=1, sinBasis=2, cosBasis=3,
linBasis=4, quadBasis=5, coshBasis=6, sinhBasis=7 } ;
enum BasisSign { Both=0, Plus=+1, Minus=-1 } ;
inline RooGaussModel() : _flatSFInt(kFALSE), _asympInt(kFALSE) { }
RooGaussModel(const char *name, const char *title, RooRealVar& x,
RooAbsReal& mean, RooAbsReal& sigma) ;
RooGaussModel(const char *name, const char *title, RooRealVar& x,
RooAbsReal& mean, RooAbsReal& sigma, RooAbsReal& msSF) ;
RooGaussModel(const char *name, const char *title, RooRealVar& x,
RooAbsReal& mean, RooAbsReal& sigma, RooAbsReal& meanSF, RooAbsReal& sigmaSF) ;
RooGaussModel(const RooGaussModel& other, const char* name=0);
virtual TObject* clone(const char* newname) const { return new RooGaussModel(*this,newname) ; }
virtual ~RooGaussModel();
virtual Int_t basisCode(const char* name) const ;
virtual Int_t getAnalyticalIntegral(RooArgSet& allVars, RooArgSet& analVars, const char* rangeName=0) const ;
virtual Double_t analyticalIntegral(Int_t code, const char* rangeName) const ;
Int_t getGenerator(const RooArgSet& directVars, RooArgSet &generateVars, Bool_t staticInitOK=kTRUE) const;
void generateEvent(Int_t code);
void advertiseFlatScaleFactorIntegral(Bool_t flag) { _flatSFInt = flag ; }
void advertiseAymptoticIntegral(Bool_t flag) { _asympInt = flag ; }
protected:
virtual Double_t evaluate() const ;
RooComplex evalCerfApprox(Double_t swt, Double_t u, Double_t c) const ;
inline RooComplex evalCerf(Double_t swt, Double_t u, Double_t c) const {
RooComplex z(swt*c,u+c);
return (z.im()>-4.0) ? RooMath::FastComplexErrFunc(z)*exp(-u*u) : evalCerfApprox(swt,u,c) ;
}
inline Double_t evalCerfRe(Double_t swt, Double_t u, Double_t c) const {
RooComplex z(swt*c,u+c);
return (z.im()>-4.0) ? RooMath::FastComplexErrFuncRe(z)*exp(-u*u) : evalCerfApprox(swt,u,c).re() ;
}
inline Double_t evalCerfIm(Double_t swt, Double_t u, Double_t c) const {
RooComplex z(swt*c,u+c);
return (z.im()>-4.0) ? RooMath::FastComplexErrFuncIm(z)*exp(-u*u) : evalCerfApprox(swt,u,c).im() ;
}
RooComplex evalCerfInt(Double_t sign, Double_t wt, Double_t tau, Double_t umin, Double_t umax, Double_t c) const ;
Bool_t _flatSFInt ;
Bool_t _asympInt ;
RooRealProxy mean ;
RooRealProxy sigma ;
RooRealProxy msf ;
RooRealProxy ssf ;
ClassDef(RooGaussModel,1)
};
#endif