physics/include/G4RDShellEMDataSet.hh

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 //
0027 //
0028 // Author: Maria Grazia Pia (Maria.Grazia.Pia@cern.ch)
0029 //
0030 // History:
0031 // -----------
0032 // 31 Jul 2001   MGP        Created
0033 //  9 Mar 2008   MGP        Cleaned up unreadable code modified by former developer
0034 //                          (Further clean-up needed)
0035 //
0036 // -------------------------------------------------------------------
0037
0038 // Class description:
0039 // Low Energy Electromagnetic Physics
0040 // Data set for an electromagnetic physics process
0041 // A strategy pattern is used to encapsulate algorithms for data interpolation
0042 // Further documentation available from http://www.ge.infn.it/geant4/lowE
0043
0044 // -------------------------------------------------------------------
0045
0046 #ifndef G4RDSHELLEMDATASET_HH
0047 #define G4RDSHELLEMDATASET_HH 1
0048
0049 #include <vector>
0050 #include <CLHEP/Units/SystemOfUnits.h>
0051
0052 #include "globals.hh"
0053 #include "G4RDVEMDataSet.hh"
0054
0055 class G4RDVDataSetAlgorithm;
0056
0057 class G4RDShellEMDataSet : public G4RDVEMDataSet
0058 {
0059 public:
0060   G4RDShellEMDataSet(G4int Z,
0061            G4RDVDataSetAlgorithm* algo,
0062            G4double eUnit=CLHEP::MeV,
0063            G4double dataUnit=CLHEP::barn);
0064
0065   virtual ~G4RDShellEMDataSet();
0066
0067   virtual G4double FindValue(G4double energy, G4int componentId=0) const;
0068
0069   virtual void PrintData(void) const;
0070
0071   virtual const G4RDVEMDataSet*  GetComponent(G4int componentId) const { return components[componentId]; }
0072   virtual void AddComponent(G4RDVEMDataSet* dataSet) { components.push_back(dataSet); }
0073   virtual size_t NumberOfComponents(void) const { return components.size(); }
0074
0075   virtual const G4DataVector& GetEnergies(G4int componentId) const { return GetComponent(componentId)->GetEnergies(0); }
0076   virtual const G4DataVector& GetData(G4int componentId) const { return GetComponent(componentId)->GetData(0); }
0077   virtual void SetEnergiesData(G4DataVector* energies, G4DataVector* data, G4int componentId);
0078
0079   virtual G4bool LoadData(const G4String& fileName);
0080   virtual G4bool SaveData(const G4String& fileName) const;
0081
0082   virtual G4double RandomSelect(G4int /*componentId = 0*/) const { return -1.; };
0083
0084 protected:
0085
0086   G4double GetUnitEnergies() const { return unitEnergies; }
0087   G4double GetUnitData() const { return unitData; }
0088   const G4RDVDataSetAlgorithm* GetAlgorithm() const { return algorithm; }
0089
0090   void CleanUpComponents(void);
0091
0092 private:
0093
0094   G4String FullFileName(const G4String& fileName) const;
0095
0096   // Hide copy constructor and assignment operator
0097   G4RDShellEMDataSet();
0098   G4RDShellEMDataSet(const G4RDShellEMDataSet& copy);
0099   G4RDShellEMDataSet& operator=(const G4RDShellEMDataSet& right);
0100
0101   std::vector<G4RDVEMDataSet*> components;          // Owned pointers
0102
0103   G4int z;
0104
0105   G4RDVDataSetAlgorithm* algorithm;           // Owned pointer
0106
0107   G4double unitEnergies;
0108   G4double unitData;
0109 };
0110 #endif /* G4SHELLEMDATASET_HH */