CrystalBall

#include <RoccoR.h>

Public Member Functions
	CrystalBall ()
void	init ()
double	pdf (double x) const
double	pdf (double x, double ks, double dm) const
double	cdf (double x) const
double	invcdf (double u) const

Public Attributes
double	m
double	s
double	a
double	n
double	B
double	C
double	D
double	N
double	NA
double	Ns
double	NC
double	F
double	G
double	k
double	cdfMa
double	cdfPa

Static Public Attributes
static const double	pi = M_PI
static const double	sqrtPiOver2 = 0.8862269254527579
static const double	sqrt2 = 1.4142135623730951

Constructor & Destructor Documentation

CrystalBall()

CrystalBall ( )

inline

                 :m(0),s(1),a(10),n(10){
    init();
    }

Member Function Documentation

cdf()

double cdf ( double  x ) const

inline

                              {
    double d = (x-m)/s;
    if(d<-a) return NC / pow(F-s*d/G, n-1);
    if(d>a) return NC * (C - pow(F+s*d/G, 1-n) );
    return Ns * (D - sqrtPiOver2 * erf(-d/sqrt2));
    }

init()

void init ( )

inline

               {
    double fa = fabs(a);
    double ex = exp(-fa*fa/2);
    double A  = pow(n/fa, n) * ex;
    double C1 = n/fa/(n-1) * ex; 
    double D1 = 2 * sqrtPiOver2 * erf(fa/sqrt2);
 
    B = n/fa-fa;
    C = (D1+2*C1)/C1;   
    D = (D1+2*C1)/2;   
 
    N = 1.0/s/(D1+2*C1); 
    k = 1.0/(n-1);  
 
    NA = N*A;       
    Ns = N*s;       
    NC = Ns*C1;     
    F = 1-fa*fa/n; 
    G = s*n/fa;    
 
    cdfMa = cdf(m-a*s);
    cdfPa = cdf(m+a*s);
    }

invcdf()

double invcdf ( double  u ) const

inline

                                 {
    if(u<cdfMa) return m + G*(F - pow(NC/u, k));
    if(u>cdfPa) return m - G*(F - pow(C-u/NC, -k) );
    return m - sqrt2 * s * boost::math::erf_inv((D - u/Ns )/sqrtPiOver2);
    }

pdf() [1/2]

double pdf ( double  x ) const

inline

                              { 
    double d=(x-m)/s;
    if(d<-a) return NA*pow(B-d, -n);
    if(d>a) return NA*pow(B+d, -n);
    return N*exp(-d*d/2);
    }

pdf() [2/2]

double pdf ( double  x, double  ks, double  dm  ) const

inline

                                                    { 
    double d=(x-m-dm)/(s*ks);
    if(d<-a) return NA/ks*pow(B-d, -n);
    if(d>a) return NA/ks*pow(B+d, -n);
    return N/ks*exp(-d*d/2);
    }

Member Data Documentation

a

double a

B

double B

C

double C

cdfMa

double cdfMa

cdfPa

double cdfPa

D

double D

F

double F

G

double G

k

double k

m

double m

n

double n

N

double N

NA

double NA

NC

double NC

Ns

double Ns

pi

constexpr double pi = M_PI

staticconstexpr

s

double s

sqrt2

constexpr double sqrt2 = 1.4142135623730951

staticconstexpr

sqrtPiOver2

constexpr double sqrtPiOver2 = 0.8862269254527579

staticconstexpr

---

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

• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Muons/interface/RoccoR.h
• /builds/cms-analysis/general/DasAnalysisSystem/Core/Installer/Core/Muons/src/RoccoR.cc