RooJohnsonSU Class Reference

#include <RooJohnsonSU.h>

Inheritance diagram for RooJohnsonSU:

Public Member Functions
	RooJohnsonSU ()
	RooJohnsonSU (const char *name, const char *title, RooAbsReal &_x, RooAbsReal &_mean, RooAbsReal &_width, RooAbsReal &_nu, RooAbsReal &_tau)
	RooJohnsonSU (const RooJohnsonSU &other, const char *name=0)
TObject *	clone (const char *newname) const
virtual	~RooJohnsonSU ()

Protected Member Functions
Double_t	evaluate () const
Int_t	getAnalyticalIntegral (RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=0) const
Double_t	analyticalIntegral (Int_t code, const char *rangeName=0) const

Protected Attributes
RooRealProxy	x
RooRealProxy	mean
RooRealProxy	width
RooRealProxy	nu
RooRealProxy	tau

Detailed Description

Definition at line 31 of file RooJohnsonSU.h.

Constructor & Destructor Documentation

RooJohnsonSU() [1/3]

RooJohnsonSU::RooJohnsonSU ( )

inline

Definition at line 33 of file RooJohnsonSU.h.

{} ;

RooJohnsonSU() [2/3]

RooJohnsonSU::RooJohnsonSU	(	const char *	name,
		const char *	title,
		RooAbsReal &	_x,
		RooAbsReal &	_mean,
		RooAbsReal &	_width,
		RooAbsReal &	_nu,
		RooAbsReal &	_tau
	)

Definition at line 39 of file RooJohnsonSU.cpp.

                                             :
RooAbsPdf(name,title), 
x("x","x",this,_x),
mean("mean","mean",this,_mean),
width("width","width",this,_width),
nu("nu","nu",this,_nu),
tau("tau","tau",this,_tau)
{ 
} 

RooJohnsonSU() [3/3]

RooJohnsonSU::RooJohnsonSU	(	const RooJohnsonSU &	other,
		const char *	name = 0
	)

Definition at line 55 of file RooJohnsonSU.cpp.

                                                                      :  
RooAbsPdf(other,name), 
x("x",this,other.x),
mean("mean",this,other.mean),
width("width",this,other.width),
nu("nu",this,other.nu),
tau("tau",this,other.tau)
{ 
} 

~RooJohnsonSU()

virtual RooJohnsonSU::~RooJohnsonSU ( )

inlinevirtual

Definition at line 42 of file RooJohnsonSU.h.

{ }

Member Function Documentation

analyticalIntegral()

Double_t RooJohnsonSU::analyticalIntegral ( Int_t  code, const char *  rangeName = 0  ) const

protected

Definition at line 94 of file RooJohnsonSU.cpp.

{
    assert(code==1) ;
    
    // calculate a few variables
    double w = exp( tau * tau);
    double omega = - nu * tau;
    double c = 0.5 * (w-1) * (w * TMath::CosH(2 * omega) + 1);
    c = pow(c, -0.5);
    double zmax = (- x.max(rangeName) + (mean + c * width * sqrt(w) * TMath::SinH(omega) )) / c / width;
    double zmin = (- x.min(rangeName) + (mean + c * width * sqrt(w) * TMath::SinH(omega) )) / c / width;
    static const Double_t pi = atan2(0.0,-1.0);
    static const Double_t PiBy2 = pi/2.0;
    static const Double_t rootPiBy2 = sqrt(PiBy2);
    
    // the integral calculation
    Double_t ret = 0;
    if(code==1){  
        //    ret = rootPiBy2 * sigx * (RooMath::erf((x.max(rangeName)-mean)/xscale)-RooMath::erf((x.min(rangeName)-mean)/xscale));
        ret =  -0.25/rootPiBy2* ( RooMath::erf( (nu*tau + TMath::ASinH( zmax ) )/(sqrt(2)*tau) )-
                                 RooMath::erf( (nu*tau + TMath::ASinH( zmin ) )/(sqrt(2)*tau) ) );
        
        /*-rootPiBy2 * ( RooMath::erf( PiBy2 * nu - ( PiBy2 * TMath::ASinH( zmax ) / tau ) ) - 
         RooMath::erf( PiBy2 * nu - ( PiBy2 * TMath::ASinH( zmin ) / tau ) ) );
         ret *= c * width;
         ret *= 0.5 / c / width / pi;*/
        //     if (gDebug>2) {
        //       cout << "Int_gauss_dx(mean=" << mean << ",sigma=" << sigma << ", xmin=" << x.min(rangeName) << ", xmax=" << x.max(rangeName) << ")=" << ret << endl ;
        //     }
    } else{
        cout << "Error in RooJohnsonSU::analyticalIntegral" << endl;
    }
    return ret ;
    
}

clone()

TObject* RooJohnsonSU::clone ( const char *  newname ) const

inline

Definition at line 41 of file RooJohnsonSU.h.

{ return new RooJohnsonSU(*this,newname); }

evaluate()

Double_t RooJohnsonSU::evaluate ( ) const

protected

Definition at line 67 of file RooJohnsonSU.cpp.

{ 
    // ENTER EXPRESSION IN TERMS OF VARIABLE ARGUMENTS HERE 
    // calculate a few variables
    double w = exp( tau * tau);
    double omega = - nu * tau;
    double c = 0.5 * (w-1) * (w * TMath::CosH(2 * omega) + 1);
    c = pow(c, -0.5);
    double z = (x - (mean + c * width * sqrt(w) * TMath::SinH(omega) )) / c / width;
    double r = -nu + TMath::ASinH(z) / tau;
    // the actual value
    double val = 1. / (c * width * 2 * TMath::Pi());
    val *= 1. / (tau * sqrt(z*z+1));
    val *= exp(-0.5 * r * r);
    // return
    return val ; 
} 

getAnalyticalIntegral()

Int_t RooJohnsonSU::getAnalyticalIntegral ( RooArgSet &  allVars, RooArgSet &  analVars, const char *  rangeName = 0  ) const

protected

Definition at line 87 of file RooJohnsonSU.cpp.

{
    if (matchArgs(allVars,analVars,x)) return 1 ;
    return 0 ;
}

Member Data Documentation

mean

RooRealProxy RooJohnsonSU::mean

protected

Definition at line 47 of file RooJohnsonSU.h.

nu

RooRealProxy RooJohnsonSU::nu

protected

Definition at line 49 of file RooJohnsonSU.h.

tau

RooRealProxy RooJohnsonSU::tau

protected

Definition at line 50 of file RooJohnsonSU.h.

width

RooRealProxy RooJohnsonSU::width

protected

Definition at line 48 of file RooJohnsonSU.h.

x

RooRealProxy RooJohnsonSU::x

protected

Definition at line 46 of file RooJohnsonSU.h.

---

The documentation for this class was generated from the following files:

• Mint/RooJohnsonSU.h
• src/Mojito/TimeDependent/RooJohnsonSU.cpp