Main Page
Namespaces
Namespace List
Namespace Members
All
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Functions
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Variables
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Typedefs
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Enumerations
a
b
c
d
e
f
g
h
i
j
k
l
m
o
p
q
r
s
t
u
v
w
z
Enumerator
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Classes
Class List
Class Index
Class Hierarchy
Class Members
All
:
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
~
Functions
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
~
Variables
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Typedefs
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Enumerations
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
Enumerator
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
y
z
Properties
_
a
d
e
f
l
m
o
p
s
t
u
v
Related Functions
:
_
a
b
c
d
e
f
g
h
i
j
k
l
m
n
o
p
q
r
s
t
u
v
w
x
Package Documentation
•
All
Classes
Namespaces
Files
Functions
Variables
Typedefs
Enumerations
Enumerator
Properties
Friends
Macros
Modules
Pages
CommonTools
Statistics
src
IncompleteGammaComplement.cc
Go to the documentation of this file.
1
#include "
CommonTools/Statistics/src/IncompleteGammaComplement.h
"
2
#include "
CommonTools/Statistics/src/GammaContinuedFraction.h
"
3
#include "
CommonTools/Statistics/src/GammaSeries.h
"
4
#include "
CommonTools/Statistics/src/GammaLn.h
"
5
6
#include "
FWCore/MessageLogger/interface/MessageLogger.h
"
7
8
#include <iostream>
9
#include <cmath>
10
11
float
IncompleteGammaComplement::value
(
float
a
,
float
x
) {
12
if
(
x
< 0.0 ||
a
<= 0.0)
13
edm::LogInfo
(
"IncompleteGammaComplement"
) <<
"IncompleteGammaComplement::invalid arguments"
;
14
if
(
x
< (
a
+ 1.0))
15
// take the complement of the series representation
16
return
1. -
GammaSeries
(
a
,
x
) * (
exp
(-
x
+
a
*
log
(
x
) -
GammaLn
(
a
)));
17
else
18
// use the continued fraction representation
19
return
GammaContinuedFraction
(
a
,
x
) * (
exp
(-
x
+
a
*
log
(
x
) -
GammaLn
(
a
)));
20
}
21
22
float
IncompleteGammaComplement::ln
(
float
a
,
float
x
) {
23
if
(
x
< 0.0 ||
a
<= 0.0)
24
edm::LogInfo
(
"IncompleteGammaComplement"
) <<
"IncompleteGammaComplement::invalid arguments"
;
25
if
(
x
< (
a
+ 1.0))
26
// take the complement of the series representation
27
return
log
(1. -
GammaSeries
(
a
,
x
) * (
exp
(-
x
+
a
*
log
(
x
) -
GammaLn
(
a
))));
28
else
29
// use the continued fraction representation
30
return
log
(
GammaContinuedFraction
(
a
,
x
)) -
x
+
a
*
log
(
x
) -
GammaLn
(
a
);
31
}
GammaSeries.h
IncompleteGammaComplement::value
static float value(float a, float x)
Definition:
IncompleteGammaComplement.cc:11
MessageLogger.h
GammaContinuedFraction.h
edm::LogInfo
Definition:
MessageLogger.h:254
GammaSeries
float GammaSeries(float a, float x)
Definition:
GammaSeries.cc:8
DDAxes::x
vertices_cff.x
x
Definition:
vertices_cff.py:29
IncompleteGammaComplement.h
a
double a
Definition:
hdecay.h:119
GammaLn
float GammaLn(float z)
Definition:
GammaLn.cc:4
GammaLn.h
dqm-mbProfile.log
log
Definition:
dqm-mbProfile.py:17
IncompleteGammaComplement::ln
static float ln(float a, float x)
Definition:
IncompleteGammaComplement.cc:22
JetChargeProducer_cfi.exp
exp
Definition:
JetChargeProducer_cfi.py:6
GammaContinuedFraction
float GammaContinuedFraction(float a, float x)
Definition:
GammaContinuedFraction.cc:9
Generated for CMSSW Reference Manual by
1.8.16