G4NeutronHPIsoData Class Reference

#include <G4NeutronHPIsoData.hh>

Public Member Functions
	G4NeutronHPIsoData ()
	~G4NeutronHPIsoData ()
G4double	GetXsec (G4double energy)
G4bool	Init (G4int A, G4int Z, G4double abun, G4String dirName, G4String aFSType)
G4bool	Init (G4int A, G4int Z, G4int M, G4double abun, G4String dirName, G4String aFSType)
void	Init (G4int A, G4int Z, G4double abun)
void	Init (G4int A, G4int Z, G4int M, G4double abun)
G4NeutronHPVector *	MakeElasticData ()
G4NeutronHPVector *	MakeFissionData ()
G4NeutronHPVector *	MakeCaptureData ()
G4NeutronHPVector *	MakeInelasticData ()
G4NeutronHPVector *	MakeChannelData ()
G4String	GetName (G4int A, G4int Z, G4String base, G4String rest)
void	FillChannelData (G4NeutronHPVector *aBuffer)
void	ThinOut (G4double precision)

Detailed Description

Definition at line 44 of file G4NeutronHPIsoData.hh.

Constructor & Destructor Documentation

G4NeutronHPIsoData()

G4NeutronHPIsoData::G4NeutronHPIsoData ( )

inline

Definition at line 48 of file G4NeutronHPIsoData.hh.

  {
    theChannelData = 0;
    theFissionData = 0;
    theCaptureData = 0;
    theElasticData = 0;
    theInelasticData = 0;
  }

~G4NeutronHPIsoData()

G4NeutronHPIsoData::~G4NeutronHPIsoData ( )

inline

Definition at line 57 of file G4NeutronHPIsoData.hh.

{if(theChannelData!=0) delete theChannelData;}

Member Function Documentation

FillChannelData()

void G4NeutronHPIsoData::FillChannelData ( G4NeutronHPVector *  aBuffer )

inline

Definition at line 86 of file G4NeutronHPIsoData.hh.

  {
    if(theChannelData!=0) throw G4HadronicException(__FILE__, __LINE__, "IsoData has channel full already!!!");
    theChannelData = new G4NeutronHPVector;
    for(G4int i=0; i<aBuffer->GetVectorLength(); i++)
    {
      theChannelData->SetPoint(i, aBuffer->GetPoint(i));
    }
  }

Referenced by G4NeutronHPChannel::UpdateData().

GetName()

G4String G4NeutronHPIsoData::GetName	(	G4int	A,
		G4int	Z,
		G4String	base,
		G4String	rest
	)

Definition at line 127 of file G4NeutronHPIsoData.cc.

  {
    G4bool dbool;
    return (theNames.GetName(A, Z, base, rest, dbool)).GetName();
  }

GetXsec()

G4double G4NeutronHPIsoData::GetXsec ( G4double  energy )

inline

Definition at line 59 of file G4NeutronHPIsoData.hh.

  {
    return std::max(0., theChannelData->GetXsec(energy));
  }

Referenced by G4NeutronHPChannel::ApplyYourself(), and G4NeutronHPChannel::GetWeightedXsec().

Init() [1/4]

void G4NeutronHPIsoData::Init ( G4int  A, G4int  Z, G4double  abun  )

inline

Definition at line 69 of file G4NeutronHPIsoData.hh.

                                              { G4int M =0;
  Init( A, Z, M, abun); }; 

Init() [2/4]

G4bool G4NeutronHPIsoData::Init ( G4int  A, G4int  Z, G4double  abun, G4String  dirName, G4String  aFSType  )

inline

Definition at line 65 of file G4NeutronHPIsoData.hh.

{ G4int M = 0 ; return Init( A, Z, M, abun, dirName, aFSType); };

Referenced by Init(), G4NeutronHPChannel::UpdateData(), and G4NeutronHPElementData::UpdateData().

Init() [3/4]

void G4NeutronHPIsoData::Init	(	G4int	A,
		G4int	Z,
		G4int	M,
		G4double	abun
	)

Definition at line 77 of file G4NeutronHPIsoData.cc.

  {
    G4String dirName;
    if(!getenv("G4NEUTRONHPDATA")) 
       throw G4HadronicException(__FILE__, __LINE__, "Please setenv G4NEUTRONHPDATA to point to the neutron cross-section files.");
    G4String baseName = getenv("G4NEUTRONHPDATA");
    dirName = baseName+"/Fission";
    //if(Z>89) 
    if(Z>87) //TK Modifed for ENDF VII.0 
    {
      //Init(A, Z, abun, dirName, "/CrossSection/");
      Init(A, Z, M, abun, dirName, "/CrossSection");
    }
    else
    {
       theChannelData = new G4NeutronHPVector;
    }
    theFissionData = theChannelData;
    theChannelData = 0; // fast fix for double delete; revisit later. @@@@@@@
    dirName = baseName+"/Capture";
    //Init(A, Z, abun, dirName, "/CrossSection/");
    Init(A, Z, M, abun, dirName, "/CrossSection");
    theCaptureData = theChannelData;
    theChannelData = 0;
    dirName = baseName+"/Elastic";
    //Init(A, Z, abun, dirName, "/CrossSection/");
    Init(A, Z, M, abun, dirName, "/CrossSection");
    theElasticData = theChannelData;
    theChannelData = 0;
    dirName = baseName+"/Inelastic";
    //Init(A, Z, abun, dirName, "/CrossSection/");
    Init(A, Z, M, abun, dirName, "/CrossSection");
    theInelasticData = theChannelData;
    theChannelData = 0;
    
//    if(theInelasticData!=0) G4cout << "Inelastic Data Statistics: "<<theInelasticData->GetVectorLength()<<G4endl;
//    if(theElasticData!=0) G4cout << "Elastic Data Statistics: "<<theElasticData->GetVectorLength()<<G4endl;
//    if(theCaptureData!=0) G4cout << "Capture Data Statistics: "<<theCaptureData->GetVectorLength()<<G4endl;
//    if(theFissionData!=0) G4cout << "Fission Data Statistics: "<<theFissionData->GetVectorLength()<<G4endl;
//  G4cout << "Inelastic data"<<G4endl;
//  if(theInelasticData!=0) theInelasticData->Dump();
//  G4cout << "Elastic data"<<G4endl;
//  if(theElasticData!=0) theElasticData->Dump();
//  G4cout << "Capture data"<<G4endl;
//  if(theCaptureData!=0) theCaptureData->Dump();
//  G4cout << "Fission data"<<G4endl;
//  if(theFissionData!=0) theFissionData->Dump();
 
  }

Init() [4/4]

G4bool G4NeutronHPIsoData::Init	(	G4int	A,
		G4int	Z,
		G4int	M,
		G4double	abun,
		G4String	dirName,
		G4String	aFSType
	)

Definition at line 37 of file G4NeutronHPIsoData.cc.

  {
    theChannelData = 0;
    G4double abundance = abun/100.;
    G4String filename;
    G4bool result = true;
    //G4NeutronHPDataUsed aFile = theNames.GetName(A, Z, dirName, aFSType, result);
    G4NeutronHPDataUsed aFile = theNames.GetName(A, Z, M, dirName, aFSType, result);
    filename = aFile.GetName();
//    if(filename=="") return false;
    std::ifstream theChannel(filename);
    
    if(Z==1 && (aFile.GetZ()!=Z || std::abs(aFile.GetA()-A)>0.0001) )
    {
      if(getenv("NeutronHPNamesLogging")) G4cout << "Skipped = "<< filename <<" "<<A<<" "<<Z<<G4endl;
      //080901 TKDB No more necessary below protection, cross sections set to 0 in G4NeutronHPNames
      //And below two lines causes trouble with G4PhysicsVector
      //theChannel.close();
      //return false;
    }
    if(!theChannel) {theChannel.close(); return false;}
    // accommodating deficiencie of some compilers
    if(theChannel.eof()) {theChannel.close(); return false;} 
    if(!theChannel) {theChannel.close(); return false;}
    G4int dummy; 
    theChannel >> dummy >> dummy;
    theChannelData = new G4NeutronHPVector;
    G4int nData;
    theChannel >> nData;
    theChannelData->Init(theChannel, nData, eV, abundance*barn);
//    G4cout << "Channel Data Statistics: "<<theChannelData->GetVectorLength()<<G4endl;
//    G4cout << "Channel data"<<G4endl;
//     G4int hpw;
//     G4cin >> hpw;
//    theChannelData->Dump();
    theChannel.close();
    return result;
  }

MakeCaptureData()

G4NeutronHPVector * G4NeutronHPIsoData::MakeCaptureData ( )

inline

Definition at line 77 of file G4NeutronHPIsoData.hh.

    {return theCaptureData;}

Referenced by G4NeutronHPElementData::UpdateData().

MakeChannelData()

G4NeutronHPVector * G4NeutronHPIsoData::MakeChannelData ( )

inline

Definition at line 81 of file G4NeutronHPIsoData.hh.

    {return theChannelData;}

Referenced by G4NeutronHPChannel::UpdateData().

MakeElasticData()

G4NeutronHPVector * G4NeutronHPIsoData::MakeElasticData ( )

inline

Definition at line 73 of file G4NeutronHPIsoData.hh.

    {return theElasticData;}

Referenced by G4NeutronHPElementData::UpdateData().

MakeFissionData()

G4NeutronHPVector * G4NeutronHPIsoData::MakeFissionData ( )

inline

Definition at line 75 of file G4NeutronHPIsoData.hh.

    {return theFissionData;}

Referenced by G4NeutronHPElementData::UpdateData().

MakeInelasticData()

G4NeutronHPVector * G4NeutronHPIsoData::MakeInelasticData ( )

inline

Definition at line 79 of file G4NeutronHPIsoData.hh.

    {return theInelasticData;}

Referenced by G4NeutronHPElementData::UpdateData().

ThinOut()

void G4NeutronHPIsoData::ThinOut ( G4double  precision )

inline

Definition at line 96 of file G4NeutronHPIsoData.hh.

  {
    if(theFissionData) theFissionData->ThinOut(precision);
    if(theCaptureData) theCaptureData->ThinOut(precision);
    if(theElasticData) theElasticData->ThinOut(precision);
    if(theInelasticData) theInelasticData->ThinOut(precision);
  }

---

The documentation for this class was generated from the following files:

• G4NeutronHPIsoData.hh
• G4NeutronHPIsoData.cc