69 double massThreshold) :
71 _mass(
"mass",
"Mass observable", this, mass),
72 _mu(
"mu",
"Location parameter of the underlying normal distribution.", this, mu),
73 _lambda(
"lambda",
"Width parameter of the underlying normal distribution (=2 lambda)", this, lambda),
74 _gamma(
"gamma",
"Shift of transformation", this, gamma),
75 _delta(
"delta",
"Scale of transformation", this, delta),
76 _massThreshold(massThreshold)
86 _mass(
"Mass", this, other._mass),
87 _mu(
"mean", this, other._mu),
88 _lambda(
"lambda", this, other._lambda),
89 _gamma(
"gamma", this, other._gamma),
90 _delta(
"delta", this, other._delta),
91 _massThreshold(other._massThreshold)
109 / (
_lambda * sqrt(1. + arg*arg))
110 * exp(-0.5 * expo * expo);
121 {ctx.at(_mass), ctx.at(_mu), ctx.at(_lambda), ctx.at(_gamma), ctx.at(_delta)}, extraArgs);
144 const double globalNorm = 1.;
148 double min = -1.E300;
157 }
else if (code ==
kMean) {
168 }
else if (code ==
kGamma) {
172 }
else if (code ==
kDelta) {
186 const double ecmin = std::erfc(std::abs(min/sqrt(2.)));
187 const double ecmax = std::erfc(std::abs(max/sqrt(2.)));
189 const double result = 0.5 * (
190 min*max < 0.0 ? 2.0 - (ecmin + ecmax)
191 : max <= 0. ? ecmax - ecmin : ecmin - ecmax
230 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
Abstract interface for all probability density functions.
Abstract base class for objects that represent a real value and implements functionality common to al...
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.
std::span< double > output()
RooBatchCompute::Config config(RooAbsArg const *arg) const
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().
void doEval(RooFit::EvalContext &) const override
Compute multiple values of the Johnson distribution.
Int_t getAnalyticalIntegral(RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=nullptr) const override
Interface function getAnalyticalIntergral advertises the analytical integrals that are supported.
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.
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...
void compute(Config cfg, Computer comp, std::span< double > output, VarSpan vars, ArgSpan extraArgs={})
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.
constexpr Double_t TwoPi()