ROOT logo
/*****************************************************************************
 * 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 "RooResolutionModel.h"
#include "RooRealProxy.h"
#include "RooComplex.h"
#include "RooMath.h"

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, RooRealVar& x, 
	       RooAbsReal& sigma, RooAbsReal& rlife, 
	       Bool_t nlo=kFALSE, Type type=Normal) ; 

  RooGExpModel(const char *name, const char *title, RooRealVar& x, 
	       RooAbsReal& sigma, RooAbsReal& rlife, 
	       RooAbsReal& srSF, 
	       Bool_t nlo=kFALSE, Type type=Normal) ; 

  RooGExpModel(const char *name, const char *title, RooRealVar& x, 
	       RooAbsReal& sigma, RooAbsReal& rlife, 
	       RooAbsReal& sigmaSF, RooAbsReal& rlifeSF,
	       Bool_t nlo=kFALSE, Type type=Normal) ; 

  RooGExpModel(const RooGExpModel& other, const char* name=0);
  virtual TObject* clone(const char* newname) const { return new RooGExpModel(*this,newname) ; }
  virtual ~RooGExpModel();
  
  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=0) 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 advertiseAsymptoticIntegral(Bool_t flag) { _asympInt = flag ; }  // added FMV,07/24/03

protected:

  Double_t logErfC(Double_t x) const ;

  //Double_t calcDecayConv(Double_t sign, Double_t tau, Double_t sig, Double_t rtau) const ;
  Double_t calcDecayConv(Double_t sign, Double_t tau, Double_t sig, Double_t rtau, Double_t fsign) const ;  
	// modified FMV,08/13/03
  RooComplex calcSinConv(Double_t sign, Double_t sig, Double_t tau, Double_t omega, Double_t rtau, Double_t fsign) const ;
  Double_t calcSinConv(Double_t sign, Double_t sig, Double_t tau, Double_t rtau, Double_t fsign) const ;
  RooComplex calcSinConvNorm(Double_t sign, Double_t tau, Double_t omega, 
	                     Double_t sig, Double_t rtau, Double_t fsign, const char* rangeName) const ; // modified FMV,07/24/03
  Double_t calcSinConvNorm(Double_t sign, Double_t tau, 
        Double_t sig, Double_t rtau, Double_t fsign, const char* rangeName) const ; // added FMV,08/18/03
  //Double_t calcSinhConv(Double_t sign, Double_t sign1, Double_t sign2, Double_t tau, Double_t dgamma, Double_t sig, Double_t rtau, Double_t fsign) const ;
  //Double_t calcCoshConv(Double_t sign, Double_t tau, Double_t dgamma, Double_t sig, Double_t rtau, Double_t fsign) const ;
  virtual Double_t evaluate() const ;
  RooComplex evalCerfApprox(Double_t swt, Double_t u, Double_t c) const ;

  // Calculate exp(-u^2) cwerf(swt*c + i(u+c)), taking care of numerical instabilities
  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) ;
  }
    
  // Calculate Re(exp(-u^2) cwerf(swt*c + i(u+c))), taking care of numerical instabilities
  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() ;
  }
  
  // Calculate Im(exp(-u^2) cwerf(swt*c + i(u+c))), taking care of numerical instabilities
  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() ;
  }

  // Calculate Re(exp(-u^2) cwerf(i(u+c)))
  // added FMV, 08/17/03
  inline Double_t evalCerfRe(Double_t u, Double_t c) const {
    Double_t expArg = u*2*c+c*c ;
    if (expArg<300) {
       return exp(expArg) * RooMath::erfc(u+c);
    } else {
       return exp(expArg+logErfC(u+c));
    }
  }

  // Calculate common normalization factors 
  // added FMV,07/24/03
  RooComplex evalCerfInt(Double_t sign, Double_t wt, Double_t tau, Double_t umin, Double_t umax, Double_t c) const ;
  Double_t evalCerfInt(Double_t sign, Double_t tau, Double_t umin, Double_t umax, Double_t c) const ;

  RooRealProxy sigma ;
  RooRealProxy rlife ;
  RooRealProxy ssf ;
  RooRealProxy rsf ;
  Bool_t _flip ;
  Bool_t _nlo ;
  Bool_t _flatSFInt ;
  Bool_t _asympInt ;  // added FMV,07/24/03

  ClassDef(RooGExpModel,1) // Gauss (x) Expontial resolution model
};

#endif





 RooGExpModel.h:1
 RooGExpModel.h:2
 RooGExpModel.h:3
 RooGExpModel.h:4
 RooGExpModel.h:5
 RooGExpModel.h:6
 RooGExpModel.h:7
 RooGExpModel.h:8
 RooGExpModel.h:9
 RooGExpModel.h:10
 RooGExpModel.h:11
 RooGExpModel.h:12
 RooGExpModel.h:13
 RooGExpModel.h:14
 RooGExpModel.h:15
 RooGExpModel.h:16
 RooGExpModel.h:17
 RooGExpModel.h:18
 RooGExpModel.h:19
 RooGExpModel.h:20
 RooGExpModel.h:21
 RooGExpModel.h:22
 RooGExpModel.h:23
 RooGExpModel.h:24
 RooGExpModel.h:25
 RooGExpModel.h:26
 RooGExpModel.h:27
 RooGExpModel.h:28
 RooGExpModel.h:29
 RooGExpModel.h:30
 RooGExpModel.h:31
 RooGExpModel.h:32
 RooGExpModel.h:33
 RooGExpModel.h:34
 RooGExpModel.h:35
 RooGExpModel.h:36
 RooGExpModel.h:37
 RooGExpModel.h:38
 RooGExpModel.h:39
 RooGExpModel.h:40
 RooGExpModel.h:41
 RooGExpModel.h:42
 RooGExpModel.h:43
 RooGExpModel.h:44
 RooGExpModel.h:45
 RooGExpModel.h:46
 RooGExpModel.h:47
 RooGExpModel.h:48
 RooGExpModel.h:49
 RooGExpModel.h:50
 RooGExpModel.h:51
 RooGExpModel.h:52
 RooGExpModel.h:53
 RooGExpModel.h:54
 RooGExpModel.h:55
 RooGExpModel.h:56
 RooGExpModel.h:57
 RooGExpModel.h:58
 RooGExpModel.h:59
 RooGExpModel.h:60
 RooGExpModel.h:61
 RooGExpModel.h:62
 RooGExpModel.h:63
 RooGExpModel.h:64
 RooGExpModel.h:65
 RooGExpModel.h:66
 RooGExpModel.h:67
 RooGExpModel.h:68
 RooGExpModel.h:69
 RooGExpModel.h:70
 RooGExpModel.h:71
 RooGExpModel.h:72
 RooGExpModel.h:73
 RooGExpModel.h:74
 RooGExpModel.h:75
 RooGExpModel.h:76
 RooGExpModel.h:77
 RooGExpModel.h:78
 RooGExpModel.h:79
 RooGExpModel.h:80
 RooGExpModel.h:81
 RooGExpModel.h:82
 RooGExpModel.h:83
 RooGExpModel.h:84
 RooGExpModel.h:85
 RooGExpModel.h:86
 RooGExpModel.h:87
 RooGExpModel.h:88
 RooGExpModel.h:89
 RooGExpModel.h:90
 RooGExpModel.h:91
 RooGExpModel.h:92
 RooGExpModel.h:93
 RooGExpModel.h:94
 RooGExpModel.h:95
 RooGExpModel.h:96
 RooGExpModel.h:97
 RooGExpModel.h:98
 RooGExpModel.h:99
 RooGExpModel.h:100
 RooGExpModel.h:101
 RooGExpModel.h:102
 RooGExpModel.h:103
 RooGExpModel.h:104
 RooGExpModel.h:105
 RooGExpModel.h:106
 RooGExpModel.h:107
 RooGExpModel.h:108
 RooGExpModel.h:109
 RooGExpModel.h:110
 RooGExpModel.h:111
 RooGExpModel.h:112
 RooGExpModel.h:113
 RooGExpModel.h:114
 RooGExpModel.h:115
 RooGExpModel.h:116
 RooGExpModel.h:117
 RooGExpModel.h:118
 RooGExpModel.h:119
 RooGExpModel.h:120
 RooGExpModel.h:121
 RooGExpModel.h:122
 RooGExpModel.h:123
 RooGExpModel.h:124
 RooGExpModel.h:125
 RooGExpModel.h:126
 RooGExpModel.h:127
 RooGExpModel.h:128
 RooGExpModel.h:129
 RooGExpModel.h:130
 RooGExpModel.h:131
 RooGExpModel.h:132
 RooGExpModel.h:133
 RooGExpModel.h:134
 RooGExpModel.h:135
 RooGExpModel.h:136
 RooGExpModel.h:137
 RooGExpModel.h:138
 RooGExpModel.h:139
 RooGExpModel.h:140
 RooGExpModel.h:141