Ignore:
Timestamp:
Sep 30, 2010, 2:47:17 PM (14 years ago)
Author:
garnier
Message:

tag geant4.9.4 beta 1 + modifs locales

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/source/processes/electromagnetic/xrays/src/G4SynchrotronRadiation.cc

    r1228 r1337  
    2525//
    2626//
    27 // $Id: G4SynchrotronRadiation.cc,v 1.5 2006/06/29 19:56:15 gunter Exp $
    28 // GEANT4 tag $Name: geant4-09-03 $
     27// $Id: G4SynchrotronRadiation.cc,v 1.6 2010/06/16 15:34:15 gcosmo Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030// --------------------------------------------------------------
     
    4444
    4545#include "G4SynchrotronRadiation.hh"
    46 // #include "G4Integrator.hh"
    4746#include "G4UnitsTable.hh"
    48 
    49 using namespace std;
    5047
    5148///////////////////////////////////////////////////////////////////////
     
    6562  fFieldPropagator = transportMgr->GetPropagatorInField();
    6663
    67   fLambdaConst = sqrt(3.0)*electron_mass_c2/
     64  fLambdaConst = std::sqrt(3.0)*electron_mass_c2/
    6865                           (2.5*fine_structure_const*eplus*c_light) ;
    6966  fEnergyConst = 1.5*c_light*c_light*eplus*hbar_Planck/electron_mass_c2  ;
     
    169166                << "  B = " << Btot/tesla << " Tesla"
    170167            << "  perpB = " << perpB/tesla << " Tesla"
    171             << "  Theta = " << Theta << " sin(Theta)=" << sin(Theta) << '\n'
     168            << "  Theta = " << Theta << " std::sin(Theta)=" << std::sin(Theta) << '\n'
    172169            << "  ptot  = " << G4BestUnit(ptot,"Energy")
    173170            << "  rho   = " << G4BestUnit(rho,"Length")
     
    262259      G4double Phi  = twopi * G4UniformRand() ;
    263260
    264       G4double dirx = sin(Teta)*cos(Phi) ,
    265                diry = sin(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) ;
    267264
    268265      G4ThreeVector gammaDirection ( dirx, diry, dirz);
     
    271268      // polarization of new gamma
    272269
    273       // G4double sx =  cos(Teta)*cos(Phi);
    274       // G4double sy =  cos(Teta)*sin(Phi);
    275       // G4double sz = -sin(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);
    276273
    277274      G4ThreeVector gammaPolarization = FieldValue.cross(gammaDirection);
     
    368365  else if(x<aa3) return       Chebyshev(aa2,aa3,cheb2,ncheb2,x);
    369366  else if(x<1-0.0000841363)
    370   { G4double y=-log(1-x);
     367  { G4double y=-std::log(1-x);
    371368        return y*Chebyshev(aa4,aa5,cheb3,ncheb3,y);
    372369  }
    373370  else
    374   { G4double y=-log(1-x);
     371  { G4double y=-std::log(1-x);
    375372        return y*Chebyshev(aa5,aa6,cheb4,ncheb4,y);
    376373  }
     
    384381  static G4bool FirstTime=true;
    385382  if(verboseLevel > 0 && FirstTime)
    386   { G4double Emean=8./(15.*sqrt(3.))*Ecr; // mean photon energy
    387     G4double E_rms=sqrt(211./675.)*Ecr; // rms of photon energy distribution
     383  { 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
    388385    G4cout << "G4SynchrotronRadiation::GetRandomEnergySR :" << '\n' << std::setprecision(4)
    389386        << "  Ecr   = "    << G4BestUnit(Ecr,"Energy") << '\n'
Note: See TracChangeset for help on using the changeset viewer.