1 #ifndef __JMATHSUPPORTKIT__
2 #define __JMATHSUPPORTKIT__
19 namespace JPP {
using namespace JMATH; }
35 const double u = x /
sigma;
52 inline double gauss(
const double x,
const double x0,
const double sigma)
54 return gauss(x - x0, sigma);
79 inline double Gauss(
const double x,
const double x0,
const double sigma)
81 return Gauss(x - x0, sigma);
96 inline double Gamma(
const double a,
const double x)
100 const int max = 1000000;
105 const double gln = lgamma(a);
118 double sum = 1.0 /
a;
121 for (
int i = 0; i != max; ++i) {
128 return sum*
exp(-x + a*
log(x) - gln);
138 double b = x + 1.0 -
a;
139 double c = 1.0 / FPMIN;
143 for (
int i = 1; i != max; ++i) {
145 const double an = -i * (i-
a);
150 if (fabs(d) < FPMIN) {
156 if (fabs(c) < FPMIN) {
162 const double del = d*
c;
167 return 1.0 -
exp(-x + a*
log(x) - gln) * h;
199 for (
size_t i = 2; i <=
n; ++i) {
202 p2 = ((2*i-1) * x*p1 - (i-1) * p0) / i;
220 if (n == 0 || n < k) {
224 if (k == 0 || n == k) {
228 const int k1 = std::min(k, n - k);
229 const int k2 = n - k1;
231 double value = k2 + 1;
233 for (
int i = k1; i != 1; --i) {
234 value *= (double) (k2 + i) / (double) i;
248 inline double poisson(
const size_t n,
const double mu)
255 return exp(n*
log(mu) - lgamma(n+1) - mu);
258 }
else if (mu == 0.0) {
260 return (n == 0 ? 1.0 : 0.0);
274 inline double Poisson(
const size_t n,
const double mu)
276 return 1.0 -
Gamma(n + 1, mu);
then fatal No hydrophone data file $HYDROPHONE_TXT fi sort gr k
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
double Gamma(const double a, const double x)
Incomplete gamma function.
static const double PI
Mathematical constants.
$WORKDIR ev_configure_dqsimulator txt echo process $DQ_SIMULATOR $i $SOURCE_HOST[$index] csh c(setenv ROOTSYS $ROOTSYS &&source $JPP_DIR/setenv.csh $JPP_DIR &&($DQ_SIMULATOR\-u\$NAME\$\-H\$SERVER\$\-M\$LOGGER\$\-d $DEBUG</dev/null > &/dev/null &))'
then JMuonMCEvt f $INPUT_FILE o $INTERMEDIATE_FILE d
double Poisson(const size_t n, const double mu)
Poisson cumulative density distribition.
then set_variable NUMBER_OF_TESTS else set_variable NUMBER_OF_TESTS fi function gauss()
double legendre(const size_t n, const double x)
Legendre polynome.
double poisson(const size_t n, const double mu)
Poisson probability density distribition.
double Gauss(const double x, const double sigma)
Normalised Gauss function.
Exception for accessing a value in a collection that is outside of its range.
then cat $TRIPOD_INITIAL<< EOF1 256877.5 4743716.7-2438.42 256815.5 4743395.0-2435.53 257096.2 4743636.0-2439.5EOFfiif[[!-f $DETECTOR]];then JEditDetector-a $DETECTOR_INITIAL-s"-1 addz -6.9"-o $DETECTOR--!eval`JPrintDetector-a $DETECTOR-O SUMMARY`for STRING in ${STRINGS[*]};do set_variable MODULE`getModule-a $DETECTOR-L"$STRING 0"`JEditDetector-a $DETECTOR-M"$MODULE setz -2.9"-o $DETECTOR--!donefiif[[!-f $TRIPOD]];then cp-p $TRIPOD_INITIAL $TRIPODfiJAcoustics.sh $DETECTOR_IDcat > acoustics_trigger_parameters txt<< EOFQ=0.0;TMax_s=0.020;numberOfHits=90;EOFJAcousticsEventBuilder.sh $DETECTOR $RUNS[*]INPUT_FILES=(`ls KM3NeT_ ${(l:8::0::0:) DETECTOR_ID}_0 *${^RUNS}_event.root`) cd $WORKDIRsource $JPP_DIR/examples/JAcoustics/acoustics-fit-toolkit.shif[!$HOMEDIR-ef $WORKDIR];then cp-p $HOMEDIR/$DETECTOR $WORKDIR cp-p $HOMEDIR/${^ACOUSTICS_KEYS}.txt $WORKDIR cp-p $HOMEDIR/${^INPUT_FILES}$WORKDIRfitimer_startinitialise stage_b 1 0 100.0e-6 10.0 0.002 0.1 0 > &stage log
double binomial(const size_t n, const size_t k)
Binomial function.
then fatal Wrong number of arguments fi set_variable DETECTOR $argv[1] set_variable STRING $argv[2] set_array QUANTILES set_variable FORMULA *[0] exp(-0.5 *(x-[1])*(x-[1])/([2]*[2]))" set_variable MODULE `getModule -a $DETECTOR -L "$STRING 0"` source JAcousticsToolkit.sh typeset -A TRIPODS get_tripods $WORKDIR/tripod.txt TRIPODS XMEAN