Changeset 1192 for trunk/source/processes/electromagnetic/lowenergy/src/G4DNAEmfietzoglouExcitationModel.cc

Timestamp:

Nov 19, 2009, 2:53:25 PM (15 years ago)

Author:

garnier

Message:

update par rapport a CVS

File:

• trunk/source/processes/electromagnetic/lowenergy/src/G4DNAEmfietzoglouExcitationModel.cc (modified) (4 diffs)

trunk/source/processes/electromagnetic/lowenergy/src/G4DNAEmfietzoglouExcitationModel.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4DNAEmfietzoglouExcitationModel.cc,v 1.5 2009/02/16 12:46:58 sincerti Exp $
-// GEANT4 tag $Name: geant4-09-03-beta-cand-01 $
+// $Id: G4DNAEmfietzoglouExcitationModel.cc,v 1.8 2009/08/13 11:32:47 sincerti Exp $
+// GEANT4 tag $Name: emlowen-V09-02-64 $
 //
 
@@ -56 +56 @@
   if (verboseLevel > 3)
 
-  G4cout << "Emfietzoglou Excitation model is constructed " << G4endl
-         << "Energy range: "
-         << lowEnergyLimit / eV << " eV - "
-         << highEnergyLimit / MeV << " MeV"
-         << G4endl;
- 
+  if( verboseLevel>0 ) 
+  { 
+    G4cout << "Emfietzoglou Excitation model is constructed " << G4endl
+           << "Energy range: "
+           << lowEnergyLimit / eV << " eV - "
+           << highEnergyLimit / MeV << " MeV"
+           << G4endl;
+  } 
 }
 
@@ -99 +101 @@
 
   //
-  
-  G4cout << "Emfietzoglou Excitation model is initialized " << G4endl
-         << "Energy range: "
-         << LowEnergyLimit() / eV << " eV - "
-         << HighEnergyLimit() / MeV << " MeV"
-         << G4endl;
+  if( verboseLevel>0 ) 
+  { 
+    G4cout << "Emfietzoglou Excitation model is initialized " << G4endl
+           << "Energy range: "
+           << LowEnergyLimit() / eV << " eV - "
+           << HighEnergyLimit() / MeV << " MeV"
+           << G4endl;
+  }
 
   if(!isInitialised) 
@@ -136 +140 @@
             const G4Material* material = couple->GetMaterial();
 
-            size_t j = material->GetNumberOfElements();
-            while (j>0)
+            if (material->GetName() == "G4_WATER") 
             {
-               j--;
-               const G4Element* element(material->GetElement(j));
-               if (element->GetZ() == 8.)
-               {
-                  G4double density = material->GetAtomicNumDensityVector()[j];
-                  if (density > 0.) 
-                  { 
-                    flagMaterialIsWater = true; 
-                    densityWater = density; 
-                    
-                    if (verboseLevel > 3) 
-                    G4cout << "Water material is found with density(cm^-3)=" << density/(cm*cm*cm) << G4endl;
-                  }
-               }
+              G4double density = material->GetAtomicNumDensityVector()[1];
+              flagMaterialIsWater = true; 
+              densityWater = density; 
+              
+              if (verboseLevel > 3) 
+              G4cout << "****** Water material is found with density(cm^-3)=" << density/(cm*cm*cm) << G4endl;
             }
   
           }
-   } // if(numOfCouples>0)
+
+    } // if(numOfCouples>0)
 
   } // if (theCoupleTable)