Changeset 961 for trunk/source/processes/electromagnetic/highenergy/src/G4GammaConversionToMuons.cc

Timestamp:

Apr 6, 2009, 12:21:12 PM (15 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/highenergy/src/G4GammaConversionToMuons.cc (modified) (6 diffs)

trunk/source/processes/electromagnetic/highenergy/src/G4GammaConversionToMuons.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4GammaConversionToMuons.cc,v 1.4 2006/06/29 19:32:40 gunter Exp $
-// GEANT4 tag $Name: geant4-09-01-patch-02 $
+// $Id: G4GammaConversionToMuons.cc,v 1.7 2008/10/16 14:29:48 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 //         ------------ G4GammaConversionToMuons physics process ------
@@ -51 +51 @@
     HighestEnergyLimit(1e21*eV), // ok to 1e21eV=1e12GeV, then LPM suppression
     CrossSecFactor(1.)
-{ }
+{ 
+  SetProcessSubType(15);
+}
 
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.....
@@ -261 +263 @@
     G4double xxp=1.-4./3.*xPM; // the main xPlus dependence
     result=xxp*log(W)*LogWmaxInv;
-    if(result>1.)
-    { G4cout << "error in dSigxPlusGen, result=" << result << " is >1" << '\n';
-      exit(10);
+    if(result>1.) {
+      G4cout << "G4GammaConversionToMuons::PostStepDoIt WARNING:"
+             << " in dSigxPlusGen, result=" << result << " > 1" << G4endl;
     }
   }
@@ -285 +287 @@
       f1=(1.-2.*xPM+4.*xPM*t*(1.-t)) / (1.+C1/(t*t));
       if(f1<0 || f1> f1_max) // should never happend
-      { G4cout << "outside allowed range f1=" << f1 << G4endl;
-        exit(1);
-      }
+        {
+          G4cout << "G4GammaConversionToMuons::PostStepDoIt WARNING:"
+                 << "outside allowed range f1=" << f1 << " is set to zero"
+                 << G4endl;
+          f1 = 0.0;
+        }
     }
     while ( G4UniformRand()*f1_max > f1);
@@ -299 +304 @@
       f2=1.-2.*xPM+4.*xPM*t*(1.-t)*(1.+cos(2.*psi));
       if(f2<0 || f2> f2_max) // should never happend
-      { G4cout << "outside allowed range f2=" << f2 << G4endl;
-        exit(1);
-      }
+        {
+          G4cout << "G4GammaConversionToMuons::PostStepDoIt WARNING:"
+                 << "outside allowed range f2=" << f2 << " is set to zero"
+                 << G4endl;
+          f2 = 0.0;
+        }
     }
     while ( G4UniformRand()*f2_max > f2);
@@ -387 +395 @@
 void G4GammaConversionToMuons::PrintInfoDefinition()
 {
-  G4String comments ="gamma->mu+mu- Bethe Heitler process.\n";
+  G4String comments ="gamma->mu+mu- Bethe Heitler process, SubType= ";
   G4cout << G4endl << GetProcessName() << ":  " << comments
-         << "        good cross section parametrization from "
+         << GetProcessSubType() << G4endl;
+  G4cout << "        good cross section parametrization from "
          << G4BestUnit(LowestEnergyLimit,"Energy")
          << " to " << HighestEnergyLimit/GeV << " GeV for all Z." << G4endl;