Changeset 19639 for trunk/Mars

Timestamp:

09/15/19 14:23:22 (5 years ago)

Author:

tbretz

Message:

Added a random function to simulate a gaussian shape with long tails which do not extend over pi/2

Location:

trunk/Mars/mbase

Files:

• MMath.cc (modified) (1 diff)
• MMath.h (modified) (1 diff)

trunk/Mars/mbase/MMath.cc

@@ -901 +901 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Returns an exponential distribution
+//
+// -tau*log(x)  with x=[0;1[ and y=[0;1[
+//
 Double_t MMath::RndmExp(Double_t tau)
 {
-    // returns an exponential deviate.
-    //
-    //          exp( -t/tau )
-
     const Double_t x = gRandom->Rndm(); // uniform on ] 0, 1 ]
 
     return -tau * TMath::Log(x);       // convert to exponential distribution
 }
+
+// --------------------------------------------------------------------------
+//
+// Returns a distribution according to
+//
+//  atanh((y+1)/x*sigma) with x and y [-1;1] and x^2+y^2=1
+//
+Double_t MMath::RndmPSF(Double_t sigma)
+{
+    double x, y;
+    gRandom->Circle(x, y, 1);
+    return TMath::ATanH((y+1)/x*sigma);
+}

trunk/Mars/mbase/MMath.h

@@ -91 +91 @@
 
     Double_t RndmExp(Double_t tau);
+    Double_t RndmPSF(Double_t sigma=1);
 }