include/Geant4/G4PixeShellDataSet.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 as G4EMShellDataSe
0033 //  9 Mar 2008   MGP        Cleaned up unreadable code modified by former developer
0034 //                          (Further clean-up needed)
0035 // 31 Jul 2008   MGP        Revised and renamed to G4PixeShellDataSet
0036 //
0037 // -------------------------------------------------------------------
0038
0039 // Class description:
0040 // Shell data set
0041 // Applies a Composite design pattern for data library management
0042
0043 // -------------------------------------------------------------------
0044
0045 #ifndef  G4PIXESHELLDATASET_HH
0046 #define  G4PIXESHELLDATASET_HH 1
0047
0048 #include <vector>
0049 #include <CLHEP/Units/SystemOfUnits.h>
0050
0051 #include "globals.hh"
0052 #include "G4IDataSet.hh"
0053
0054 class G4IInterpolator;
0055
0056 class G4PixeShellDataSet : public G4IDataSet
0057 {
0058 public:
0059   G4PixeShellDataSet(G4int Z,
0060              G4IInterpolator* algo,
0061              //const G4String& particleTpye="proton",
0062                      const G4String& modelK="ecpssr",
0063              const G4String& modelL="ecpssr",
0064              const G4String& modelM="ecpssr",
0065              G4double eUnit=CLHEP::MeV,
0066              G4double dataUnit=CLHEP::barn);
0067
0068   virtual ~G4PixeShellDataSet();
0069
0070   virtual G4double FindValue(G4double energy, G4int componentId=0) const;
0071
0072   virtual void PrintData(void) const;
0073
0074   virtual const G4IDataSet* GetComponent(G4int componentId) const
0075   { return components[componentId]; }
0076
0077   virtual void AddComponent(G4IDataSet* dataSet)
0078   { components.push_back(dataSet); }
0079
0080   virtual size_t NumberOfComponents(void) const
0081   { return components.size(); }
0082
0083   virtual const G4DataVector& GetEnergies(G4int componentId) const
0084   { return GetComponent(componentId)->GetEnergies(0); }
0085
0086   virtual const G4DataVector& GetData(G4int componentId) const
0087   { return GetComponent(componentId)->GetData(0); }
0088
0089   virtual void SetEnergiesData(G4DataVector* energies,
0090                    G4DataVector* data,
0091                    G4int componentId);
0092
0093   virtual G4bool LoadData(const G4String& fileName);
0094   virtual G4bool SaveData(const G4String& fileName) const;
0095
0096   virtual G4double RandomSelect(G4int /*componentId = 0*/) const {return -1.;};
0097
0098 protected:
0099
0100   G4double GetUnitEnergies() const { return unitEnergies; }
0101   G4double GetUnitData() const { return unitData; }
0102   const G4IInterpolator* GetAlgorithm() const { return algorithm; }
0103
0104   void CleanUpComponents(void);
0105
0106 private:
0107
0108   G4String FullFileName(const G4String& particleType,
0109             const G4String& subShell) const;
0110   G4int TranslateShell(const G4String& subShell) const;
0111
0112   // Hide copy constructor and assignment operator
0113   G4PixeShellDataSet();
0114   G4PixeShellDataSet(const G4PixeShellDataSet& copy);
0115   G4PixeShellDataSet& operator=(const G4PixeShellDataSet& right);
0116
0117   std::vector<G4IDataSet*> components;          // Owned pointers
0118
0119   G4int z;
0120   G4IInterpolator* algorithm;           // Owned pointer
0121   // G4String particle;
0122   // G4String crossModelK;
0123   // G4String crossModelL;
0124   // G4String crossModelM;
0125   std::vector<G4String> crossModel;
0126   G4double unitEnergies;
0127   G4double unitData;
0128   std::vector<G4String> shellName;
0129   std::vector<G4String> subShellName;
0130
0131 };
0132 #endif /* G4PIXESHELLDATASET_HH */