00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026 #ifndef G4NeutronHPJENDLHEData_h
00027 #define G4NeutronHPJENDLHEData_h 1
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038 #include "G4VCrossSectionDataSet.hh"
00039 #include "G4DynamicParticle.hh"
00040 #include "G4Neutron.hh"
00041 #include "G4Element.hh"
00042 #include "G4ParticleDefinition.hh"
00043 #include "G4PhysicsVector.hh"
00044 #include <map>
00045
00046 class G4NeutronHPJENDLHEData : public G4VCrossSectionDataSet
00047 {
00048 public:
00049
00050 G4NeutronHPJENDLHEData();
00051 G4NeutronHPJENDLHEData( G4String , G4ParticleDefinition* );
00052
00053
00054 ~G4NeutronHPJENDLHEData();
00055
00056 G4bool IsApplicable(const G4DynamicParticle*, const G4Element*);
00057
00058 G4bool IsZAApplicable( const G4DynamicParticle* , G4double , G4double )
00059 { return false; }
00060
00061 G4double GetCrossSection(const G4DynamicParticle*, const G4Element*, G4double aT);
00062
00063 void BuildPhysicsTable(const G4ParticleDefinition&);
00064
00065 void DumpPhysicsTable(const G4ParticleDefinition&);
00066
00067 private:
00068
00069 std::vector< G4bool > vElement;
00070
00071 std::map< G4int , std::map< G4int , G4PhysicsVector* >* > mIsotope;
00072
00073 G4bool isThisInMap ( G4int , G4int );
00074 G4bool isThisNewIsotope ( G4int z , G4int a ) { return !( isThisInMap( z , a ) ); };
00075 G4PhysicsVector* readAFile ( std::fstream* );
00076 void registAPhysicsVector ( G4int , G4int , G4PhysicsVector* );
00077
00078 G4double getXSfromThisIsotope ( G4int , G4int , G4double );
00079
00080 G4String reactionName;
00081 G4String particleName;
00082 };
00083
00084 #endif