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

Timestamp:

May 28, 2009, 4:26:57 PM (15 years ago)

Author:

garnier

Message:

maj sur la beta de geant 4.9.3

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.10 2008/12/09 15:08:13 pandola Exp $
-// GEANT4 tag $Name: geant4-09-02 $
+// $Id: G4PenelopeBremsstrahlungContinuous.cc,v 1.11 2008/12/15 09:23:06 pandola Exp $
+// GEANT4 tag $Name: geant4-09-03-beta-cand-01 $
 // 
 // --------------------------------------------------------------
@@ -97 +97 @@
  G4String dirFile = pathString + "/penelope/" + name;
  std::ifstream file(dirFile);
- std::filebuf* lsdp = file.rdbuf();
- if (!(lsdp->is_open()))
+ if (!file.is_open())
      {
       G4String excep = "G4PenelopeBremsstrahlungContinuous - data file " + name + " not found!";
@@ -112 +111 @@
      for (size_t j=0;j<NumberofKPoints;j++){
        file >> a1;
-       ReducedCS[i][j]=a1*cm2; //coversion present in Penelope source
+       ReducedCS[i][j]=a1*cm2; 
      }
      //3) read the total cross section, in cm2
      file >> a1;
-     TotalCS[i]=a1*cm2; //conversion present in Penelope source
+     TotalCS[i]=a1*cm2; 
      // 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);