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 43 of file G4NeutronHPIsoData.hh.

Constructor & Destructor Documentation

G4NeutronHPIsoData::G4NeutronHPIsoData ( )

inline

Definition at line 47 of file G4NeutronHPIsoData.hh.

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

G4NeutronHPIsoData::~G4NeutronHPIsoData ( )

inline

Definition at line 56 of file G4NeutronHPIsoData.hh.

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

Member Function Documentation

void G4NeutronHPIsoData::FillChannelData ( G4NeutronHPVector *  aBuffer )

inline

Definition at line 85 of file G4NeutronHPIsoData.hh.

References G4NeutronHPVector::GetPoint(), G4NeutronHPVector::GetVectorLength(), and G4NeutronHPVector::SetPoint().

Referenced by G4NeutronHPChannel::UpdateData().

  {
    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));
    }
  }

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

Definition at line 130 of file G4NeutronHPIsoData.cc.

References G4NeutronHPNames::GetName().

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

G4double G4NeutronHPIsoData::GetXsec ( G4double  energy )

inline

Definition at line 58 of file G4NeutronHPIsoData.hh.

References G4NeutronHPVector::GetXsec(), and G4INCL::Math::max().

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

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

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

inline

Definition at line 64 of file G4NeutronHPIsoData.hh.

References Init().

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

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

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

Definition at line 38 of file G4NeutronHPIsoData.cc.

References python.hepunit::barn, python.hepunit::eV, G4cout, G4endl, G4NeutronHPDataUsed::GetA(), G4NeutronHPManager::GetDataStream(), G4NeutronHPManager::GetInstance(), G4NeutronHPNames::GetName(), G4NeutronHPDataUsed::GetName(), G4NeutronHPDataUsed::GetZ(), and G4NeutronHPVector::Init().

  {
    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);
   std::istringstream theChannel(filename,std::ios::in);
   G4NeutronHPManager::GetInstance()->GetDataStream(filename,theChannel);
    
    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;
  }

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

inline

Definition at line 68 of file G4NeutronHPIsoData.hh.

References Init().

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

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

Definition at line 80 of file G4NeutronHPIsoData.cc.

References Init().

  {
    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();

  }

G4NeutronHPVector* G4NeutronHPIsoData::MakeCaptureData ( )

inline

Definition at line 76 of file G4NeutronHPIsoData.hh.

Referenced by G4NeutronHPElementData::UpdateData().

    {return theCaptureData;}

G4NeutronHPVector* G4NeutronHPIsoData::MakeChannelData ( )

inline

Definition at line 80 of file G4NeutronHPIsoData.hh.

Referenced by G4NeutronHPChannel::UpdateData().

    {return theChannelData;}

G4NeutronHPVector* G4NeutronHPIsoData::MakeElasticData ( )

inline

Definition at line 72 of file G4NeutronHPIsoData.hh.

Referenced by G4NeutronHPElementData::UpdateData().

    {return theElasticData;}

G4NeutronHPVector* G4NeutronHPIsoData::MakeFissionData ( )

inline

Definition at line 74 of file G4NeutronHPIsoData.hh.

Referenced by G4NeutronHPElementData::UpdateData().

    {return theFissionData;}

G4NeutronHPVector* G4NeutronHPIsoData::MakeInelasticData ( )

inline

Definition at line 78 of file G4NeutronHPIsoData.hh.

Referenced by G4NeutronHPElementData::UpdateData().

    {return theInelasticData;}

void G4NeutronHPIsoData::ThinOut ( G4double  precision )

inline

Definition at line 95 of file G4NeutronHPIsoData.hh.

References G4NeutronHPVector::ThinOut().

  {
    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