Geant4 Cross Reference |
1 // 2 // ******************************************************************** 3 // * License and Disclaimer * 4 // * * 5 // * The Geant4 software is copyright of the Copyright Holders of * 6 // * the Geant4 Collaboration. It is provided under the terms and * 7 // * conditions of the Geant4 Software License, included in the file * 8 // * LICENSE and available at http://cern.ch/geant4/license . These * 9 // * include a list of copyright holders. * 10 // * * 11 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file LICENSE and URL above * 16 // * for the full disclaimer and the limitation of liability. * 17 // * * 18 // * This code implementation is the result of the scientific and * 19 // * technical work of the GEANT4 collaboration. * 20 // * By using, copying, modifying or distributing the software (or * 21 // * any work based on the software) you agree to acknowledge its * 22 // * use in resulting scientific publications, and indicate your * 23 // * acceptance of all terms of the Geant4 Software license. * 24 // ******************************************************************** 25 // 26 // 27 // ------------------------------------------------------------------- 28 // 29 // GEANT4 header file 30 // 31 // File name: G4LevelReader 32 // 33 // Author: V.Ivanchenko 34 // 35 // Creation date: 4 January 2012 36 // 37 // Modifications: 38 // 39 // ------------------------------------------------------------------- 40 // 41 // Helper class to read Geant4 nuclear level data 42 // 43 44 #ifndef G4LEVELREADER_HH 45 #define G4LEVELREADER_HH 1 46 47 #include "globals.hh" 48 #include "G4LevelManager.hh" 49 #include <iosfwd> 50 #include <vector> 51 52 class G4NuclearLevelData; 53 54 class G4LevelReader 55 { 56 57 public: 58 59 explicit G4LevelReader(G4NuclearLevelData*); 60 61 // create run manager using G4LEVELGAMMADATA data for Z and A 62 const G4LevelManager* CreateLevelManager(G4int Z, G4int A); 63 64 // create run manager using whatever data 65 const G4LevelManager* MakeLevelManager(G4int Z, G4int A, 66 const G4String& filename); 67 68 inline void SetVerbose(G4int val); 69 70 G4LevelReader(const G4LevelReader & right) = delete; 71 const G4LevelReader& operator=(const G4LevelReader &right) = delete; 72 G4bool operator==(const G4LevelReader &right) const = delete; 73 G4bool operator!=(const G4LevelReader &right) const = delete; 74 75 private: 76 77 G4bool ReadData(std::istringstream& dataFile, G4double& x); 78 79 G4bool ReadDataItem(std::istream& dataFile, G4double& x); 80 81 G4bool ReadDataItem(std::istream& dataFile, G4float& x); 82 83 G4bool ReadDataItem(std::istream& dataFile, G4int& x); 84 85 const std::vector<G4float>* NormalizedICCProbability(G4int Z); 86 87 const G4LevelManager* LevelManager(G4int Z, G4int A, 88 std::ifstream& infile); 89 90 G4NuclearLevelData* fData; 91 92 G4double fEnergy = 0.; 93 G4double fTransEnergy = 0.; 94 G4double fTime = 0.; 95 G4double fTimeFactor; 96 G4double fSpin = 0.; 97 98 G4float fProb = 0.f; 99 G4float fAlpha = 0.f; 100 G4float fAlphaMax; 101 G4float fRatio = 0.f; 102 G4float fNorm1 = 0.f; 103 G4float fICC[10] = {0.f}; 104 105 G4int nbufmax = 20; 106 G4int nbuf1 = 14; 107 G4int nbuf2 = 8; 108 G4int count1 = 0; 109 G4int count2 = 0; 110 111 G4int fVerbose = 1; 112 G4int fLevelMax = 632; 113 G4int fTransMax = 145; 114 G4int ntrans = 0; 115 G4int i1 = 0; 116 G4int i2 = 0; 117 G4int k = 0; 118 G4int kk = 0; 119 G4int tnum = 0; 120 121 char buffer[20] = {' '}; 122 char buff1[14] = {' '}; 123 char buff2[8] = {' '}; 124 125 std::vector<G4double> vEnergy; 126 std::vector<G4int> vSpin; 127 std::vector<const G4NucLevel*> vLevel; 128 129 std::vector<G4int> vTrans; 130 std::vector<G4float> vRatio; 131 std::vector<G4float> vGammaCumProbability; 132 std::vector<G4float> vGammaProbability; 133 std::vector<const std::vector<G4float>*> vShellProbability; 134 135 G4String fPol = " "; 136 G4String fDirectory; 137 }; 138 139 inline void G4LevelReader::SetVerbose(G4int val) 140 { 141 fVerbose = val; 142 } 143 144 #endif 145