Changeset 9226 for trunk/MagicSoft/Mars/mbase

Timestamp:

01/17/09 14:52:43 (16 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mbase

Files:

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

trunk/MagicSoft/Mars/mbase/MMath.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MMath.cc,v 1.43 2009-01-14 12:31:36 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MMath.cc,v 1.44 2009-01-17 14:52:37 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -20 +20 @@
 !   Author(s): Thomas Bretz  3/2004 <mailto:tbretz@astro.uni-wuerzburg.de>
 !
-!   Copyright: MAGIC Software Development, 2000-2005
+!   Copyright: MAGIC Software Development, 2000-2009
 !
 !
@@ -50 +50 @@
 #ifndef ROOT_TComplex
 #include <TComplex.h>
+#endif
+
+#ifndef ROOT_TRandom
+#include <TRandom.h>  // gRandom in RndmExp
 #endif
 
@@ -875 +879 @@
     e = error.Atof();
 }
+
+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
+}

trunk/MagicSoft/Mars/mbase/MMath.h

@@ -87 +87 @@
 
     void Format(Double_t &v, Double_t &e);
+
+    Double_t RndmExp(Double_t tau);
 }