Changeset 961 for trunk/source/processes/electromagnetic/lowenergy/src/G4CompositeEMDataSet.cc

Timestamp:

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

Author:

garnier

Message:

update processes

File:

• trunk/source/processes/electromagnetic/lowenergy/src/G4CompositeEMDataSet.cc (modified) (2 diffs)

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

@@ -25 +25 @@
 //
 //
-// $Id: G4CompositeEMDataSet.cc,v 1.9 2006/06/29 19:38:46 gunter Exp $
-// GEANT4 tag $Name: geant4-09-01-patch-02 $
+// $Id: G4CompositeEMDataSet.cc,v 1.13 2008/03/17 13:40:53 pia Exp $
+// GEANT4 tag $Name: geant4-09-02-ref-02 $
 //
 // Author: Maria Grazia Pia (Maria.Grazia.Pia@cern.ch)
@@ -42 +42 @@
 #include <sstream>
 
-                                                G4CompositeEMDataSet :: G4CompositeEMDataSet(G4VDataSetAlgorithm* argAlgorithm, G4double argUnitEnergies, G4double argUnitData, G4int argMinZ, G4int argMaxZ)
-:
- algorithm(argAlgorithm),
- unitEnergies(argUnitEnergies),
- unitData(argUnitData),
- minZ(argMinZ),
- maxZ(argMaxZ)
+G4CompositeEMDataSet::G4CompositeEMDataSet(G4VDataSetAlgorithm* argAlgorithm, 
+                                           G4double argUnitEnergies, 
+                                           G4double argUnitData, 
+                                           G4int argMinZ, 
+                                           G4int argMaxZ)
+  :
+  algorithm(argAlgorithm),
+  unitEnergies(argUnitEnergies),
+  unitData(argUnitData),
+  minZ(argMinZ),
+  maxZ(argMaxZ)
 {
- if (algorithm == 0) 
-  G4Exception("G4CompositeEMDataSet::G4CompositeEMDataSet - interpolation == 0");
+  if (algorithm == 0) 
+    G4Exception("G4CompositeEMDataSet::G4CompositeEMDataSet - interpolation == 0");
 }
 
 
 
-                                                G4CompositeEMDataSet :: ~G4CompositeEMDataSet()
+G4CompositeEMDataSet::~G4CompositeEMDataSet()
 {
- CleanUpComponents();
+  CleanUpComponents();
+  if (algorithm) delete algorithm;
+}
 
- if (algorithm)
-  delete algorithm;
+
+G4double G4CompositeEMDataSet::FindValue(G4double argEnergy, G4int argComponentId) const
+{
+  const G4VEMDataSet* component(GetComponent(argComponentId));
+ 
+  if (component) return component->FindValue(argEnergy);
+
+  std::ostringstream message;
+  message << "G4CompositeEMDataSet::FindValue - component " << argComponentId << " not found";
+ 
+  G4Exception(message.str().c_str());
+ 
+  return 0.;
+}
+
+void G4CompositeEMDataSet::PrintData(void) const
+{
+  const size_t n(NumberOfComponents());
+
+  G4cout << "The data set has " << n << " components" << G4endl;
+  G4cout << G4endl;
+ 
+  size_t i(0);
+ 
+  while (i<n)
+    {
+      G4cout << "--- Component " << i << " ---" << G4endl;
+      GetComponent(i)->PrintData();
+      i++;
+    }
+}
+
+void G4CompositeEMDataSet::SetEnergiesData(G4DataVector* argEnergies, G4DataVector* argData, G4int argComponentId)
+{
+  G4VEMDataSet * component(components[argComponentId]);
+ 
+  if (component)
+    {
+      component->SetEnergiesData(argEnergies, argData, 0);
+      return;
+    }
+
+  std::ostringstream message;
+  message << "G4CompositeEMDataSet::SetEnergiesData - component " << argComponentId << " not found";
+ 
+  G4Exception(message.str().c_str());
+}
+
+G4bool G4CompositeEMDataSet::LoadData(const G4String& argFileName)
+{
+  CleanUpComponents(); 
+
+  for (G4int z(minZ); z<maxZ; z++)
+    {
+      G4VEMDataSet* component = new G4EMDataSet(z, algorithm->Clone(), unitEnergies, unitData);
+      if (!component->LoadData(argFileName))
+        {
+          delete component;
+          return false;
+        }
+      AddComponent(component);
+    }
+  return true;
 }
 
 
 
+G4bool G4CompositeEMDataSet::SaveData(const G4String& argFileName) const
+{
+  for (G4int z=minZ; z<maxZ; z++)
+    {
+      const G4VEMDataSet* component(GetComponent(z-minZ));
+  
+      if (!component)
+        {
+          std::ostringstream message;
+          message << "G4CompositeEMDataSet::SaveData - component " << (z-minZ) << " not found";
+          G4Exception(message.str().c_str());
+        }
 
+      if (!component->SaveData(argFileName))
+        return false;
+    }
+ 
+  return true;
+}
 
-
-G4double                                        G4CompositeEMDataSet :: FindValue(G4double argEnergy, G4int argComponentId) const
+void G4CompositeEMDataSet::CleanUpComponents(void)
 {
- const G4VEMDataSet * component(GetComponent(argComponentId));
- 
- if (component)
-  return component->FindValue(argEnergy);
-
- std::ostringstream message;
- message << "G4CompositeEMDataSet::FindValue - component " << argComponentId << " not found";
- 
- G4Exception(message.str().c_str());
- 
- return 0.;
+  while (!components.empty())
+    {
+      if (components.back())
+        delete components.back();
+      components.pop_back();
+    }
 }
 
 
-
-
-
-void                                            G4CompositeEMDataSet :: PrintData(void) const
+G4double G4CompositeEMDataSet::RandomSelect(G4int componentId) const
 {
- const size_t n(NumberOfComponents());
-
- G4cout << "The data set has " << n << " components" << G4endl;
- G4cout << G4endl;
- 
- size_t i(0);
- 
- while (i<n)
- {
-  G4cout << "--- Component " << i << " ---" << G4endl;
-  GetComponent(i)->PrintData();
-  i++;
- }
+  G4double value = 0.;
+  if (componentId >= 0 && componentId < (G4int)components.size())
+    {
+      const G4VEMDataSet* dataSet = GetComponent(componentId);
+      value = dataSet->RandomSelect();
+    }
+  return value;
 }
-
-
-
-
-
-void                                            G4CompositeEMDataSet :: SetEnergiesData(G4DataVector * argEnergies, G4DataVector * argData, G4int argComponentId)
-{
- G4VEMDataSet * component(components[argComponentId]);
- 
- if (component)
- {
-  component->SetEnergiesData(argEnergies, argData, 0);
-  return;
- }
-
- std::ostringstream message;
- message << "G4CompositeEMDataSet::SetEnergiesData - component " << argComponentId << " not found";
- 
- G4Exception(message.str().c_str());
-}
-
-
-
-
-
-G4bool                                          G4CompositeEMDataSet :: LoadData(const G4String & argFileName)
-{
- CleanUpComponents(); 
-
- for (G4int z(minZ); z<maxZ; z++)
- {
-  G4VEMDataSet * component=new G4EMDataSet(z, algorithm->Clone(), unitEnergies, unitData);
-  if (!component->LoadData(argFileName))
-  {
-   delete component;
-   return false;
-  }
-  
-  AddComponent(component);
- }
- 
- return true;
-}
-
-
-
-G4bool                                          G4CompositeEMDataSet :: SaveData(const G4String & argFileName) const
-{
- for (G4int z(minZ); z<maxZ; z++)
- {
-  const G4VEMDataSet * component(GetComponent(z-minZ));
-  
-  if (!component)
-  {
-   std::ostringstream message;
-   message << "G4CompositeEMDataSet::SaveData - component " << (z-minZ) << " not found";
- 
-   G4Exception(message.str().c_str());
-  }
-
-  if (!component->SaveData(argFileName))
-   return false;
- }
- 
- return true;
-}
-
-
-
-
-
-void                                            G4CompositeEMDataSet :: CleanUpComponents(void)
-{
- while (!components.empty())
- {
-  if (components.back())
-   delete components.back();
-
-  components.pop_back();
- }
-}
-