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/G4Scintillation.cc

    r1315 r1337  
    2525//
    2626//
    27 // $Id: G4Scintillation.cc,v 1.31 2010/05/27 20:49:40 gum Exp $
    28 // GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
     27// $Id: G4Scintillation.cc,v 1.32 2010/06/16 15:34:15 gcosmo Exp $
     28// GEANT4 tag $Name: geant4-09-04-beta-01 $
    2929//
    3030////////////////////////////////////////////////////////////////////////
     
    7171#include "G4Scintillation.hh"
    7272
    73 using namespace std;
    74 
    7573/////////////////////////
    7674// Class Implementation 
     
    214212        if (MeanNumberOfPhotons > 10.)
    215213        {
    216           G4double sigma = ResolutionScale * sqrt(MeanNumberOfPhotons);
     214          G4double sigma = ResolutionScale * std::sqrt(MeanNumberOfPhotons);
    217215          NumPhotons = G4int(G4RandGauss::shoot(MeanNumberOfPhotons,sigma)+0.5);
    218216        }
     
    281279                                          GetConstProperty("YIELDRATIO");
    282280                 if ( ExcitationRatio == 1.0 ) {
    283                     Num = G4int (min(YieldRatio,1.0) * NumPhotons);
     281                    Num = G4int (std::min(YieldRatio,1.0) * NumPhotons);
    284282                 }
    285283                 else {
    286                     Num = G4int (min(ExcitationRatio,1.0) * NumPhotons);
     284                    Num = G4int (std::min(ExcitationRatio,1.0) * NumPhotons);
    287285                 }
    288286                 ScintillationTime   = aMaterialPropertiesTable->
     
    330328
    331329                G4double cost = 1. - 2.*G4UniformRand();
    332                 G4double sint = sqrt((1.-cost)*(1.+cost));
     330                G4double sint = std::sqrt((1.-cost)*(1.+cost));
    333331
    334332                G4double phi = twopi*G4UniformRand();
    335                 G4double sinp = sin(phi);
    336                 G4double cosp = cos(phi);
     333                G4double sinp = std::sin(phi);
     334                G4double cosp = std::cos(phi);
    337335
    338336                G4double px = sint*cosp;
     
    355353
    356354                phi = twopi*G4UniformRand();
    357                 sinp = sin(phi);
    358                 cosp = cos(phi);
     355                sinp = std::sin(phi);
     356                cosp = std::cos(phi);
    359357
    360358                photonPolarization = cosp * photonPolarization + sinp * perp;
     
    392390                if (ScintillationRiseTime==0.0) {
    393391                   deltaTime = deltaTime -
    394                           ScintillationTime * log( G4UniformRand() );
     392                          ScintillationTime * std::log( G4UniformRand() );
    395393                } else {
    396394                   deltaTime = deltaTime +
     
    633631          // make sure the envelope function is
    634632          // always larger than the bi-exponential
    635           G4double t = -1.0*tau2*log(1-ran1);
     633          G4double t = -1.0*tau2*std::log(1-ran1);
    636634          G4double g = d*single_exp(t,tau2);
    637635          if (ran2 <= bi_exp(t,tau1,tau2)/g) return t;
Note: See TracChangeset for help on using the changeset viewer.