/*****************************************************************************
 * Project: RooFit                                                           *
 * Package: RooFitModels                                                     *
 *    File: $Id: RooGExpModel.h,v 1.16 2007/05/11 09:13:07 verkerke Exp $
 * Authors:                                                                  *
 *   WV, Wouter Verkerke, UC Santa Barbara, verkerke@slac.stanford.edu       *
 *   DK, David Kirkby,    UC Irvine,         dkirkby@uci.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_GEXP_MODEL
#define ROO_GEXP_MODEL

#include "Rtypes.h"
#include "RooResolutionModel.h"
#include "RooRealProxy.h"

#include <complex>

class RooGExpModel : public RooResolutionModel {
public:

  enum RooGExpBasis { noBasis=0, expBasisMinus= 1, expBasisSum= 2, expBasisPlus= 3,
                       sinBasisMinus=11, sinBasisSum=12, sinBasisPlus=13,
                                 cosBasisMinus=21, cosBasisSum=22, cosBasisPlus=23,
             sinhBasisMinus=31,sinhBasisSum=32,sinhBasisPlus=33,
             coshBasisMinus=41,coshBasisSum=42,coshBasisPlus=43} ;



  enum BasisType { none=0, expBasis=1, sinBasis=2, cosBasis=3, sinhBasis=4, coshBasis=5 } ;
  enum BasisSign { Both=0, Plus=+1, Minus=-1 } ;
  enum Type { Normal, Flipped };

  // Constructors, assignment etc
  inline RooGExpModel() {
    // coverity[UNINIT_CTOR]
  }

  RooGExpModel(const char *name, const char *title, RooAbsRealLValue& x,
          RooAbsReal& mean, RooAbsReal& sigma, RooAbsReal& rlife,
          RooAbsReal& meanSF, RooAbsReal& sigmaSF, RooAbsReal& rlifeSF,
          bool nlo=false, Type type=Normal) ;

  RooGExpModel(const char *name, const char *title, RooAbsRealLValue& x,
          RooAbsReal& sigma, RooAbsReal& rlife,
          bool nlo=false, Type type=Normal) ;

  RooGExpModel(const char *name, const char *title, RooAbsRealLValue& x,
          RooAbsReal& sigma, RooAbsReal& rlife,
          RooAbsReal& srSF,
          bool nlo=false, Type type=Normal) ;

  RooGExpModel(const char *name, const char *title, RooAbsRealLValue& x,
          RooAbsReal& sigma, RooAbsReal& rlife,
          RooAbsReal& sigmaSF, RooAbsReal& rlifeSF,
          bool nlo=false, Type type=Normal) ;



  RooGExpModel(const RooGExpModel& other, const char* name=nullptr);
  TObject* clone(const char* newname) const override { return new RooGExpModel(*this,newname) ; }
  ~RooGExpModel() override;

  Int_t basisCode(const char* name) const 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 getGenerator(const RooArgSet& directVars, RooArgSet &generateVars, bool staticInitOK=true) const override;
  void generateEvent(Int_t code) override;

  void advertiseFlatScaleFactorIntegral(bool flag) { _flatSFInt = flag ; }

  void advertiseAsymptoticIntegral(bool flag) { _asympInt = flag ; }  // added FMV,07/24/03

protected:
  double evaluate() const override ;

private:
  //double calcDecayConv(double sign, double tau, double sig, double rtau) const ;
  double calcDecayConv(double sign, double tau, double sig, double rtau, double fsign) const ;
   // modified FMV,08/13/03
  std::complex<double> calcSinConv(double sign, double sig, double tau, double omega, double rtau, double fsign) const ;
  double calcSinConv(double sign, double sig, double tau, double rtau, double fsign) const ;
  std::complex<double> calcSinConvNorm(double sign, double tau, double omega,
                        double sig, double rtau, double fsign, const char* rangeName) const ; // modified FMV,07/24/03
  double calcSinConvNorm(double sign, double tau,
        double sig, double rtau, double fsign, const char* rangeName) const ; // added FMV,08/18/03
  //double calcSinhConv(double sign, double sign1, double sign2, double tau, double dgamma, double sig, double rtau, double fsign) const ;
  //double calcCoshConv(double sign, double tau, double dgamma, double sig, double rtau, double fsign) const ;

  static std::complex<double> evalCerfApprox(double swt, double u, double c);

  // Calculate common normalization factors
  // added FMV,07/24/03
  std::complex<double> evalCerfInt(double sign, double wt, double tau, double umin, double umax, double c) const ;
  double evalCerfInt(double sign, double tau, double umin, double umax, double c) const ;

  RooRealProxy _mean;
  RooRealProxy sigma ;
  RooRealProxy rlife ;
  RooRealProxy _meanSF;
  RooRealProxy ssf ;
  RooRealProxy rsf ;

  bool _flip ;
  bool _nlo ;
  bool _flatSFInt ;
  bool _asympInt ;  // added FMV,07/24/03

  ClassDefOverride(RooGExpModel,2) // Gauss (x) Exponential resolution model
};

#endif