Changeset 1007 for trunk/source/processes/electromagnetic/lowenergy/src/G4PenelopeBremsstrahlungContinuous.cc

Timestamp:

Apr 20, 2009, 5:54:05 PM (15 years ago)

Author:

garnier

Message:

update to geant4.9.2

File:

• trunk/source/processes/electromagnetic/lowenergy/src/G4PenelopeBremsstrahlungContinuous.cc (modified) (4 diffs)

trunk/source/processes/electromagnetic/lowenergy/src/G4PenelopeBremsstrahlungContinuous.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4PenelopeBremsstrahlungContinuous.cc,v 1.11 2008/12/15 09:23:06 pandola Exp $
-// GEANT4 tag $Name: geant4-09-02-ref-02 $
+// $Id: G4PenelopeBremsstrahlungContinuous.cc,v 1.10 2008/12/09 15:08:13 pandola Exp $
+// GEANT4 tag $Name: geant4-09-02 $
 // 
 // --------------------------------------------------------------
@@ -97 +97 @@
  G4String dirFile = pathString + "/penelope/" + name;
  std::ifstream file(dirFile);
- if (!file.is_open())
+ std::filebuf* lsdp = file.rdbuf();
+ if (!(lsdp->is_open()))
      {
       G4String excep = "G4PenelopeBremsstrahlungContinuous - data file " + name + " not found!";
@@ -111 +112 @@
      for (size_t j=0;j<NumberofKPoints;j++){
        file >> a1;
-       ReducedCS[i][j]=a1*cm2; 
+       ReducedCS[i][j]=a1*cm2; //coversion present in Penelope source
      }
      //3) read the total cross section, in cm2
      file >> a1;
-     TotalCS[i]=a1*cm2; 
+     TotalCS[i]=a1*cm2; //conversion present in Penelope source
      // Check closing item
      file >> a1;
      if (a1 != ((G4double) -1))
        {
-         G4String excep = "G4PenelopeBremsstrahlungContinuous - Check the bremms data file " 
-           + name;
+         G4String excep = "G4PenelopeBremsstrahlungContinuous - Check the bremms data file " + name;
          G4Exception(excep);
        }
@@ -229 +229 @@
   
   //Global x-section factor
-  G4double Fact=Zmat*Zmat*(energy+electron_mass_c2)*(energy+electron_mass_c2)/
-    (energy*(energy+2.0*electron_mass_c2));
+  G4double Fact=Zmat*Zmat*(energy+electron_mass_c2)*(energy+electron_mass_c2)/(energy*(energy+2.0*electron_mass_c2));
   Fact *= PositronCorrection(energy);