Changeset 961 for trunk/source/processes/electromagnetic/standard/src/G4eBremsstrahlungModel.cc

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/standard/src/G4eBremsstrahlungModel.cc (modified) (8 diffs)

trunk/source/processes/electromagnetic/standard/src/G4eBremsstrahlungModel.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4eBremsstrahlungModel.cc,v 1.39 2007/05/23 08:47:35 vnivanch Exp $
-// GEANT4 tag $Name:  $
+// $Id: G4eBremsstrahlungModel.cc,v 1.43 2008/11/13 19:28:58 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // -------------------------------------------------------------------
@@ -54 +54 @@
 // 27-03-06  Fix calculation of fl parameter at low energy (energy loss) (VI)
 // 15-02-07  correct LPMconstant by a factor 2, thanks to G. Depaola (mma)
+// 09-09-08  MigdalConstant increased in (2pi)^2 times (A.Schaelicke) 
 //
 // Class Description:
@@ -85 +86 @@
     isElectron(true),
     probsup(1.0),
-    MigdalConstant(classic_electr_radius*electron_Compton_length*electron_Compton_length/pi),
+    MigdalConstant(classic_electr_radius*electron_Compton_length*electron_Compton_length*4.0*pi),
     LPMconstant(fine_structure_const*electron_mass_c2*electron_mass_c2/(4.*pi*hbarc)),
-    theLPMflag(true),
     isInitialised(false)
 {
@@ -93 +93 @@
   theGamma = G4Gamma::Gamma();
   minThreshold = 1.0*keV;
-  highKinEnergy= 100.*TeV;
-  lowKinEnergy = 1.0*keV;
-  highEnergyTh = DBL_MAX;
 }
 
@@ -164 +161 @@
   if(isInitialised) return;
 
-  if(pParticleChange)
+  if(pParticleChange) {
     fParticleChange = reinterpret_cast<G4ParticleChangeForLoss*>(pParticleChange);
-  else
+  } else {
     fParticleChange = new G4ParticleChangeForLoss();
-
+  }
   isInitialised = true;
 }
@@ -834 +831 @@
     gammaEnergy = x*kineticEnergy; 
 
-    if (theLPMflag) {
+    if (LPMFlag()) {
      // take into account the supression due to the LPM effect
       if (G4UniformRand() <= SupressionFunction(material,kineticEnergy,
@@ -879 +876 @@
 
   // stop tracking and create new secondary instead of primary
-  if(gammaEnergy > highEnergyTh) {
+  if(gammaEnergy > SecondaryThreshold()) {
     fParticleChange->ProposeTrackStatus(fStopAndKill);
     fParticleChange->SetProposedKineticEnergy(0.0);
@@ -951 +948 @@
   G4double supr = 1.0;
 
-  if (theLPMflag) {
+  if (LPMFlag()) {
 
     G4double s2lpm = LPMEnergy*gammaEnergy/totEnergySquare;