- Timestamp:
- Sep 30, 2010, 2:47:17 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/electromagnetic/xrays/src/G4SynchrotronRadiation.cc
r1228 r1337 25 25 // 26 26 // 27 // $Id: G4SynchrotronRadiation.cc,v 1. 5 2006/06/29 19:56:15 gunterExp $28 // GEANT4 tag $Name: geant4-09-0 3$27 // $Id: G4SynchrotronRadiation.cc,v 1.6 2010/06/16 15:34:15 gcosmo Exp $ 28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 29 29 // 30 30 // -------------------------------------------------------------- … … 44 44 45 45 #include "G4SynchrotronRadiation.hh" 46 // #include "G4Integrator.hh"47 46 #include "G4UnitsTable.hh" 48 49 using namespace std;50 47 51 48 /////////////////////////////////////////////////////////////////////// … … 65 62 fFieldPropagator = transportMgr->GetPropagatorInField(); 66 63 67 fLambdaConst = s qrt(3.0)*electron_mass_c2/64 fLambdaConst = std::sqrt(3.0)*electron_mass_c2/ 68 65 (2.5*fine_structure_const*eplus*c_light) ; 69 66 fEnergyConst = 1.5*c_light*c_light*eplus*hbar_Planck/electron_mass_c2 ; … … 169 166 << " B = " << Btot/tesla << " Tesla" 170 167 << " perpB = " << perpB/tesla << " Tesla" 171 << " Theta = " << Theta << " s in(Theta)=" <<sin(Theta) << '\n'168 << " Theta = " << Theta << " std::sin(Theta)=" << std::sin(Theta) << '\n' 172 169 << " ptot = " << G4BestUnit(ptot,"Energy") 173 170 << " rho = " << G4BestUnit(rho,"Length") … … 262 259 G4double Phi = twopi * G4UniformRand() ; 263 260 264 G4double dirx = s in(Teta)*cos(Phi) ,265 diry = s in(Teta)*sin(Phi) ,266 dirz = cos(Teta) ;261 G4double dirx = std::sin(Teta)*std::cos(Phi) , 262 diry = std::sin(Teta)*std::sin(Phi) , 263 dirz = std::cos(Teta) ; 267 264 268 265 G4ThreeVector gammaDirection ( dirx, diry, dirz); … … 271 268 // polarization of new gamma 272 269 273 // G4double sx = cos(Teta)*cos(Phi);274 // G4double sy = cos(Teta)*sin(Phi);275 // G4double sz = -s in(Teta);270 // G4double sx = std::cos(Teta)*std::cos(Phi); 271 // G4double sy = std::cos(Teta)*std::sin(Phi); 272 // G4double sz = -std::sin(Teta); 276 273 277 274 G4ThreeVector gammaPolarization = FieldValue.cross(gammaDirection); … … 368 365 else if(x<aa3) return Chebyshev(aa2,aa3,cheb2,ncheb2,x); 369 366 else if(x<1-0.0000841363) 370 { G4double y=- log(1-x);367 { G4double y=-std::log(1-x); 371 368 return y*Chebyshev(aa4,aa5,cheb3,ncheb3,y); 372 369 } 373 370 else 374 { G4double y=- log(1-x);371 { G4double y=-std::log(1-x); 375 372 return y*Chebyshev(aa5,aa6,cheb4,ncheb4,y); 376 373 } … … 384 381 static G4bool FirstTime=true; 385 382 if(verboseLevel > 0 && FirstTime) 386 { G4double Emean=8./(15.*s qrt(3.))*Ecr; // mean photon energy387 G4double E_rms=s qrt(211./675.)*Ecr; // rms of photon energy distribution383 { G4double Emean=8./(15.*std::sqrt(3.))*Ecr; // mean photon energy 384 G4double E_rms=std::sqrt(211./675.)*Ecr; // rms of photon energy distribution 388 385 G4cout << "G4SynchrotronRadiation::GetRandomEnergySR :" << '\n' << std::setprecision(4) 389 386 << " Ecr = " << G4BestUnit(Ecr,"Energy") << '\n'
Note: See TracChangeset
for help on using the changeset viewer.