Updated STARLIGHT from r176 ro r188 (http://starlight.hepforge.org/svn/)

[u/mrichter/AliRoot.git] / STARLIGHT / starlight / src / gammaavm.cpp

diff --git a/STARLIGHT/starlight/src/gammaavm.cpp b/STARLIGHT/starlight/src/gammaavm.cpp
index 55cda2399649eb476c195d295e0aa7bd521bb82e..b95284d6499c675ca8d7ee455a8987a9a8d39725 100644 (file)
--- a/STARLIGHT/starlight/src/gammaavm.cpp
+++ b/STARLIGHT/starlight/src/gammaavm.cpp
@@ -20,9 +20,9 @@
 ///////////////////////////////////////////////////////////////////////////
 //
 // File and Version Information:
-// $Rev:: 164                         $: revision of last commit
-// $Author:: odjuvsla                 $: author of last commit
-// $Date:: 2013-10-06 16:18:08 +0200 #$: date of last commit
+// $Rev:: 179                         $: revision of last commit
+// $Author:: jnystrand                $: author of last commit
+// $Date:: 2014-09-11 20:58:05 +0200 #$: date of last commit
 //
 // Description:
 //    Added incoherent t2-> pt2 selection.  Following pp selection scheme
@@ -572,7 +572,11 @@ double Gammaavectormeson::pTgamma(double E)
     ereds = (E/_VMgamma_em)*(E/_VMgamma_em);
     //sqrt(3)*E/gamma_em is p_t where the distribution is a maximum
     Cm = sqrt(3.)*E/_VMgamma_em;
-
+    // If E is very small, the drawing of a pT below is extremely slow. 
+    // ==> Set pT = sqrt(3.)*E/_VMgamma_em for very small E. 
+    // Should have no observable consequences (JN, SRK 11-Sep-2014)
+    if( E < 0.0005 )return Cm; 
+ 
     //the amplitude of the p_t spectrum at the maximum
 
     if( _bbs.beam1().A()==1 && _bbs.beam2().A() != 1){