Changeset 1192 for trunk/source/processes/electromagnetic/lowenergy/src/G4DNADingfelderChargeDecreaseModel.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/G4DNADingfelderChargeDecreaseModel.cc (modified) (4 diffs)

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

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4DNADingfelderChargeDecreaseModel.cc,v 1.3 2009/02/16 11:00:11 sincerti Exp $
-// GEANT4 tag $Name: geant4-09-03-beta-cand-01 $
+// $Id: G4DNADingfelderChargeDecreaseModel.cc,v 1.6 2009/08/13 11:32:47 sincerti Exp $
+// GEANT4 tag $Name: emlowen-V09-02-64 $
 //
 
@@ -49 +49 @@
   // 4 = entering in methods
   
-  G4cout << "Dingfelder charge decrease model is constructed " << G4endl;
-
+  if( verboseLevel>0 ) 
+  { 
+    G4cout << "Dingfelder charge decrease model is constructed " << G4endl;
+  }
 }
 
@@ -185 +187 @@
   //
   
-  G4cout << "Dingfelder charge decrease model is initialized " << G4endl
-         << "Energy range: "
-         << LowEnergyLimit() / keV << " keV - "
-         << HighEnergyLimit() / MeV << " MeV for "
-         << particle->GetParticleName()
-         << G4endl;
-
+  if( verboseLevel>0 ) 
+  { 
+    G4cout << "Dingfelder charge decrease model is initialized " << G4endl
+           << "Energy range: "
+           << LowEnergyLimit() / keV << " keV - "
+           << HighEnergyLimit() / MeV << " MeV for "
+           << particle->GetParticleName()
+           << G4endl;
+  }
+  
   if(!isInitialised) 
   {
@@ -222 +227 @@
             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)