Changeset 1005 for trunk/source/processes/electromagnetic/lowenergy/src/G4PenelopeBremsstrahlungContinuous.cc
- Timestamp:
- Apr 20, 2009, 4:53:50 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/electromagnetic/lowenergy/src/G4PenelopeBremsstrahlungContinuous.cc
r991 r1005 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4PenelopeBremsstrahlungContinuous.cc,v 1.1 0 2008/12/09 15:08:13pandola Exp $27 // GEANT4 tag $Name: geant4-09-02 $26 // $Id: G4PenelopeBremsstrahlungContinuous.cc,v 1.11 2008/12/15 09:23:06 pandola Exp $ 27 // GEANT4 tag $Name: geant4-09-02-ref-02 $ 28 28 // 29 29 // -------------------------------------------------------------- … … 97 97 G4String dirFile = pathString + "/penelope/" + name; 98 98 std::ifstream file(dirFile); 99 std::filebuf* lsdp = file.rdbuf(); 100 if (!(lsdp->is_open())) 99 if (!file.is_open()) 101 100 { 102 101 G4String excep = "G4PenelopeBremsstrahlungContinuous - data file " + name + " not found!"; … … 112 111 for (size_t j=0;j<NumberofKPoints;j++){ 113 112 file >> a1; 114 ReducedCS[i][j]=a1*cm2; //coversion present in Penelope source113 ReducedCS[i][j]=a1*cm2; 115 114 } 116 115 //3) read the total cross section, in cm2 117 116 file >> a1; 118 TotalCS[i]=a1*cm2; //conversion present in Penelope source117 TotalCS[i]=a1*cm2; 119 118 // Check closing item 120 119 file >> a1; 121 120 if (a1 != ((G4double) -1)) 122 121 { 123 G4String excep = "G4PenelopeBremsstrahlungContinuous - Check the bremms data file " + name; 122 G4String excep = "G4PenelopeBremsstrahlungContinuous - Check the bremms data file " 123 + name; 124 124 G4Exception(excep); 125 125 } … … 229 229 230 230 //Global x-section factor 231 G4double Fact=Zmat*Zmat*(energy+electron_mass_c2)*(energy+electron_mass_c2)/(energy*(energy+2.0*electron_mass_c2)); 231 G4double Fact=Zmat*Zmat*(energy+electron_mass_c2)*(energy+electron_mass_c2)/ 232 (energy*(energy+2.0*electron_mass_c2)); 232 233 Fact *= PositronCorrection(energy); 233 234
Note: See TracChangeset
for help on using the changeset viewer.