Changeset 963 for trunk/source/processes/optical/src/G4OpRayleigh.cc

Timestamp:

Apr 6, 2009, 12:34:39 PM (15 years ago)

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/optical/src/G4OpRayleigh.cc (modified) (10 diffs)

trunk/source/processes/optical/src/G4OpRayleigh.cc

@@ -25 +25 @@
 //
 //
-// $Id: G4OpRayleigh.cc,v 1.14 2006/06/29 21:08:54 gunter Exp $
-// GEANT4 tag $Name:  $
+// $Id: G4OpRayleigh.cc,v 1.17 2008/10/24 19:51:12 gum Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // 
@@ -57 +57 @@
 
 #include "G4ios.hh"
+#include "G4OpProcessSubType.hh"
+
 #include "G4OpRayleigh.hh"
 
@@ -78 +80 @@
            : G4VDiscreteProcess(processName, type)
 {
+        SetProcessSubType(fOpRayleigh);
 
         thePhysicsTable = 0;
@@ -248 +251 @@
         const G4Material* aMaterial = aTrack.GetMaterial();
 
-        G4double thePhotonMomentum = aParticle->GetTotalMomentum();
+        G4double thePhotonEnergy = aParticle->GetTotalEnergy();
 
         G4double AttenuationLength = DBL_MAX;
@@ -258 +261 @@
            AttenuationLength =
                 (*thePhysicsTable)(aMaterial->GetIndex())->
-                           GetValue(thePhotonMomentum, isOutRange);
+                           GetValue(thePhotonEnergy, isOutRange);
         }
         else {
@@ -270 +273 @@
              if(AttenuationLengthVector){
                AttenuationLength = AttenuationLengthVector ->
-                                    GetProperty(thePhotonMomentum);
+                                    GetProperty(thePhotonEnergy);
              }
              else{
@@ -305 +308 @@
 
         // Retrieve vectors for refraction index
-        // and photon momentum from the material properties table
+        // and photon energy from the material properties table
 
         G4MaterialPropertyVector* Rindex = aMPT->GetProperty("RINDEX");
@@ -325 +328 @@
            while (++(*Rindex)) {
 
-                e = (Rindex->GetPhotonMomentum());
+                e = (Rindex->GetPhotonEnergy());
 
                 refraction_index = Rindex->GetProperty();
@@ -332 +335 @@
 
                 if (verboseLevel>0) {
-                        G4cout << Rindex->GetPhotonMomentum() << " MeV\t";
+                        G4cout << Rindex->GetPhotonEnergy() << " MeV\t";
                         G4cout << xlambda << " mm\t";
                 }
@@ -347 +350 @@
                 }
                 RayleighScatteringLengths->
-                        InsertValues(Rindex->GetPhotonMomentum(), Dist);
+                        InsertValues(Rindex->GetPhotonEnergy(), Dist);
            }