67 double massThreshold) :
69 _mass(
"mass",
"Mass observable", this, mass),
70 _mu(
"mu",
"Location parameter of the underlying normal distribution.", this, mu),
71 _lambda(
"lambda",
"Width parameter of the underlying normal distribution (=2 lambda)", this, lambda),
72 _gamma(
"gamma",
"Shift of transformation", this,
gamma),
73 _delta(
"delta",
"Scale of transformation", this, delta),
74 _massThreshold(massThreshold)
84 _mass(
"Mass", this, other._mass),
85 _mu(
"mean", this, other._mu),
86 _lambda(
"lambda", this, other._lambda),
87 _gamma(
"gamma", this, other._gamma),
88 _delta(
"delta", this, other._delta),
89 _massThreshold(other._massThreshold)
108 *
exp(-0.5 * expo * expo);
143 const double globalNorm = 1.;
147 double min = -1.E300;
150 double argMin, argMax;
155 }
else if (code ==
kMean) {
166 }
else if (code ==
kGamma) {
170 }
else if (code ==
kDelta) {
187 const double result = 0.5 * (
188 min*max < 0.0 ? 2.0 - (ecmin + ecmax)
189 : max <= 0. ? ecmax - ecmin : ecmin - ecmax
228 throw std::logic_error(
"Generation in other variables not yet implemented.");
Option_t Option_t TPoint TPoint const char GetTextMagnitude GetFillStyle GetLineColor GetLineWidth GetMarkerStyle GetTextAlign GetTextColor GetTextSize void char Point_t Rectangle_t WindowAttributes_t Float_t Float_t Float_t Int_t Int_t UInt_t UInt_t Rectangle_t result
RooAbsReal is the common abstract base class for objects that represent a real value and implements f...
bool matchArgs(const RooArgSet &allDeps, RooArgSet &numDeps, const RooArgProxy &a) const
Utility function for use in getAnalyticalIntegral().
RooArgSet is a container object that can hold multiple RooAbsArg objects.
virtual void compute(cudaStream_t *, Computer, RestrictArr, size_t, const VarVector &, const ArgVector &={})=0
auto & at(RooAbsArg const *arg, RooAbsArg const *=nullptr)
double evaluate() const override
Evaluate this PDF / function / constant. Needs to be overridden by all derived classes.
double analyticalIntegral(Int_t code, const char *rangeName=nullptr) const override
Implements the actual analytical integral(s) advertised by getAnalyticalIntegral.
void generateEvent(Int_t code) override
Generate events based on code obtained by getGenerator().
Int_t getAnalyticalIntegral(RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=nullptr) const override
Interface function getAnalyticalIntergral advertises the analytical integrals that are supported.
void computeBatch(cudaStream_t *, double *output, size_t nEvents, RooFit::Detail::DataMap const &) const override
Compute multiple values of the Johnson distribution.
Int_t getGenerator(const RooArgSet &directVars, RooArgSet &generateVars, bool staticInitOK=true) const override
Advertise which kind of direct event generation is supported.
static TRandom * randomGenerator()
Return a pointer to a singleton random-number generator implementation.
double max(const char *rname=nullptr) const
Query upper limit of range. This requires the payload to be RooAbsRealLValue or derived.
double min(const char *rname=nullptr) const
Query lower limit of range. This requires the payload to be RooAbsRealLValue or derived.
virtual Double_t Gaus(Double_t mean=0, Double_t sigma=1)
Samples a random number from the standard Normal (Gaussian) Distribution with the given mean and sigm...
double erfc(double x)
Complementary error function.
RVec< PromoteType< T > > asinh(const RVec< T > &v)
RVec< PromoteType< T > > abs(const RVec< T > &v)
RVec< PromoteType< T > > sinh(const RVec< T > &v)
RVec< PromoteType< T > > exp(const RVec< T > &v)
VecExpr< UnaryOp< Sqrt< T >, VecExpr< A, T, D >, T >, T, D > sqrt(const VecExpr< A, T, D > &rhs)
R__EXTERN RooBatchComputeInterface * dispatchCUDA
R__EXTERN RooBatchComputeInterface * dispatchCPU
This dispatch pointer points to an implementation of the compute library, provided one has been loade...
void checkRangeOfParameters(const RooAbsReal *callingClass, std::initializer_list< const RooAbsReal * > pars, double min=-std::numeric_limits< double >::max(), double max=std::numeric_limits< double >::max(), bool limitsInAllowedRange=false, std::string const &extraMessage="")
Check if the parameters have a range, and warn if the range extends below / above the set limits.
static constexpr double gauss
constexpr Double_t TwoPi()