From 67e3efa5ba72f20b6744b994f3c8e8e4f96de3b3 Mon Sep 17 00:00:00 2001
From: nemu <nemu>
Date: Tue, 25 Mar 2008 12:51:19 +0000
Subject: [PATCH] minor change for the case of identical sigmas

---
 src/classes/PTheory.cpp | 2 ++
 1 file changed, 2 insertions(+)

diff --git a/src/classes/PTheory.cpp b/src/classes/PTheory.cpp
index 41cc9997..3251e8ca 100644
--- a/src/classes/PTheory.cpp
+++ b/src/classes/PTheory.cpp
@@ -1071,6 +1071,8 @@ double PTheory::SkewedGauss(register double t, const vector<double>& paramValues
 
   if ((zp >= 25.0) || (zm >= 25.0)) // needed to prevent crash of 1F1
     skg = 2.0e300;
+  else if (val[2] == val[3]) // sigma+ == sigma- -> Gaussian
+    skg = TMath::Cos(phase+freq*t) * gp;
   else
     skg = TMath::Cos(phase+freq*t) * (wm*gm + wp*gp) +
           TMath::Sin(phase+freq*t) * (wm*gm*2.0*zm/SQRT_PI*gsl_sf_hyperg_1F1(0.5,1.5,zm*zm) -