G4ParticleHPIsoData Class Reference

#include <G4ParticleHPIsoData.hh>

Public Member Functions
void	FillChannelData (G4ParticleHPVector *aBuffer)
	G4ParticleHPIsoData ()
G4String	GetName (G4int A, G4int Z, G4String base, G4String rest)
G4double	GetXsec (G4double energy)
void	Init (G4int A, G4int Z, G4double abun, G4ParticleDefinition *projectile, const char *dataDirVariable)
G4bool	Init (G4int A, G4int Z, G4double abun, G4String dirName, G4String aFSType)
void	Init (G4int A, G4int Z, G4int M, G4double abun, G4ParticleDefinition *projectile, const char *dataDirVariable)
G4bool	Init (G4int A, G4int Z, G4int M, G4double abun, G4String dirName, G4String aFSType)
G4ParticleHPVector *	MakeCaptureData ()
G4ParticleHPVector *	MakeChannelData ()
G4ParticleHPVector *	MakeElasticData ()
G4ParticleHPVector *	MakeFissionData ()
G4ParticleHPVector *	MakeInelasticData ()
void	ThinOut (G4double precision)
	~G4ParticleHPIsoData ()

Private Attributes
G4ParticleHPVector *	theCaptureData
G4ParticleHPVector *	theChannelData
G4ParticleHPVector *	theElasticData
G4String	theFileName
G4ParticleHPVector *	theFissionData
G4ParticleHPVector *	theInelasticData
G4ParticleHPNames	theNames

Detailed Description

Definition at line 45 of file G4ParticleHPIsoData.hh.

Constructor & Destructor Documentation

G4ParticleHPIsoData()

G4ParticleHPIsoData::G4ParticleHPIsoData ( )

inline

Definition at line 49 of file G4ParticleHPIsoData.hh.

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

References theCaptureData, theChannelData, theElasticData, theFissionData, and theInelasticData.

~G4ParticleHPIsoData()

G4ParticleHPIsoData::~G4ParticleHPIsoData ( )

inline

Definition at line 58 of file G4ParticleHPIsoData.hh.

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

References theChannelData.

Member Function Documentation

FillChannelData()

void G4ParticleHPIsoData::FillChannelData ( G4ParticleHPVector *  aBuffer )

inline

Definition at line 87 of file G4ParticleHPIsoData.hh.

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

References G4ParticleHPVector::GetPoint(), G4ParticleHPVector::GetVectorLength(), G4ParticleHPVector::Hash(), G4ParticleHPVector::SetPoint(), and theChannelData.

Referenced by G4ParticleHPChannel::UpdateData().

GetName()

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

Definition at line 170 of file G4ParticleHPIsoData.cc.

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

References A, G4ParticleHPNames::GetName(), theNames, and Z.

GetXsec()

G4double G4ParticleHPIsoData::GetXsec ( G4double  energy )

inline

Definition at line 60 of file G4ParticleHPIsoData.hh.

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

References G4INCL::KinematicsUtils::energy(), G4ParticleHPVector::GetXsec(), G4INCL::Math::max(), and theChannelData.

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

Init() [1/4]

void G4ParticleHPIsoData::Init ( G4int  A, G4int  Z, G4double  abun, G4ParticleDefinition *  projectile, const char *  dataDirVariable  )

inline

Definition at line 70 of file G4ParticleHPIsoData.hh.

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

References A, Init(), M, and Z.

Init() [2/4]

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

inline

Definition at line 66 of file G4ParticleHPIsoData.hh.

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

References A, Init(), M, and Z.

Referenced by Init(), G4ParticleHPChannel::UpdateData(), and G4ParticleHPElementData::UpdateData().

Init() [3/4]

void G4ParticleHPIsoData::Init	(	G4int	A,
		G4int	Z,
		G4int	M,
		G4double	abun,
		G4ParticleDefinition *	projectile,
		const char *	dataDirVariable
	)

Definition at line 90 of file G4ParticleHPIsoData.cc.

  {
 
  G4String particleName;
  if ( projectile == G4Neutron::Neutron() ) {
     ;
  } else if ( projectile == G4Proton::Proton() ) {
     particleName = "Proton";
  } else if ( projectile == G4Deuteron::Deuteron() ) {
     particleName = "Deuteron";
  } else if ( projectile == G4Triton::Triton() ) {
     particleName = "Triton";
  } else if ( projectile == G4He3::He3() ) {
     particleName = "He3";
  } else if ( projectile == G4Alpha::Alpha() ) {
     particleName = "Alpha";
  } else {
     G4String message("G4ParticleHPInelastic may only be called for neutron, proton, deuteron, triton, He3 or alpha, while it is called for " + projectile->GetParticleName());
     throw G4HadronicException(__FILE__, __LINE__,message.c_str());
  }
 
  G4String baseName;
  if ( std::getenv( dataDirVariable ) ) {
     baseName = std::getenv( dataDirVariable );
  } else {
     baseName = std::getenv( "G4PARTICLEHPDATA" );
     baseName += "/" + particleName;
  }
    
   // G4String baseName = getenv(dataDirVariable);
    G4String dirName;
    if( projectile == G4Neutron::Neutron() ){
      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 G4ParticleHPVector;
    }
      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();
 
  }

References A, G4Alpha::Alpha(), G4Deuteron::Deuteron(), G4ParticleDefinition::GetParticleName(), G4He3::He3(), Init(), M, G4Neutron::Neutron(), G4Proton::Proton(), theCaptureData, theChannelData, theElasticData, theFissionData, theInelasticData, G4Triton::Triton(), and Z.

Init() [4/4]

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

Definition at line 39 of file G4ParticleHPIsoData.cc.

  {
    theChannelData = 0;
    G4double abundance = abun/100.;
    G4String filename;
    G4bool result = true;
    //G4ParticleHPDataUsed aFile = theNames.GetName(A, Z, dirName, aFSType, result);
    G4ParticleHPDataUsed 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);
    G4ParticleHPManager::GetInstance()->GetDataStream(filename,theChannel);
 
#ifdef G4PHPDEBUG
    if(std::getenv("G4ParticleHPDebug")) G4cout << "G4ParticleHPIsoData::Init = "<< filename <<" "<< A << " " << Z <<G4endl;
#endif
    
    if(Z==1 && (aFile.GetZ()!=Z || std::abs(aFile.GetA()-A)>0.0001) )
    {
      if(std::getenv("G4ParticleHPDebug")) G4cout << "Skipped = "<< filename <<" "<<A<<" "<<Z<<G4endl;
      //080901 TKDB No more necessary below protection, cross sections set to 0 in G4ParticleHPNames
      //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 G4ParticleHPVector;
    G4int nData;
    theChannel >> nData;
    theChannelData->Init(theChannel, nData, CLHEP::eV, abundance*CLHEP::barn);
//    G4cout << "Channel Data Statistics: "<<theChannelData->GetVectorLength()<<G4endl;
//    G4cout << "Channel data"<<G4endl;
//     G4int hpw;
//     G4cin >> hpw;
//    theChannelData->Dump();
//    theChannel.close();
    return result;
  }

References A, CLHEP::barn, CLHEP::eV, G4cout, G4endl, G4ParticleHPDataUsed::GetA(), G4ParticleHPManager::GetDataStream(), G4ParticleHPManager::GetInstance(), G4ParticleHPDataUsed::GetName(), G4ParticleHPNames::GetName(), G4ParticleHPDataUsed::GetZ(), G4ParticleHPVector::Init(), M, theChannelData, theNames, and Z.

MakeCaptureData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeCaptureData ( )

inline

Definition at line 78 of file G4ParticleHPIsoData.hh.

    {return theCaptureData;}

References theCaptureData.

Referenced by G4ParticleHPElementData::UpdateData().

MakeChannelData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeChannelData ( )

inline

Definition at line 82 of file G4ParticleHPIsoData.hh.

    {return theChannelData;}

References theChannelData.

Referenced by G4ParticleHPChannel::UpdateData().

MakeElasticData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeElasticData ( )

inline

Definition at line 74 of file G4ParticleHPIsoData.hh.

    {return theElasticData;}

References theElasticData.

Referenced by G4ParticleHPElementData::UpdateData().

MakeFissionData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeFissionData ( )

inline

Definition at line 76 of file G4ParticleHPIsoData.hh.

    {return theFissionData;}

References theFissionData.

Referenced by G4ParticleHPElementData::UpdateData().

MakeInelasticData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeInelasticData ( )

inline

Definition at line 80 of file G4ParticleHPIsoData.hh.

    {return theInelasticData;}

References theInelasticData.

Referenced by G4ParticleHPElementData::UpdateData().

ThinOut()

void G4ParticleHPIsoData::ThinOut ( G4double  precision )

inline

Definition at line 98 of file G4ParticleHPIsoData.hh.

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

References theCaptureData, theElasticData, theFissionData, theInelasticData, and G4ParticleHPVector::ThinOut().

Field Documentation

theCaptureData

G4ParticleHPVector* G4ParticleHPIsoData::theCaptureData

private

Definition at line 109 of file G4ParticleHPIsoData.hh.

Referenced by G4ParticleHPIsoData(), Init(), MakeCaptureData(), and ThinOut().

theChannelData

G4ParticleHPVector* G4ParticleHPIsoData::theChannelData

private

Definition at line 112 of file G4ParticleHPIsoData.hh.

Referenced by FillChannelData(), G4ParticleHPIsoData(), GetXsec(), Init(), MakeChannelData(), and ~G4ParticleHPIsoData().

theElasticData

G4ParticleHPVector* G4ParticleHPIsoData::theElasticData

private

Definition at line 110 of file G4ParticleHPIsoData.hh.

Referenced by G4ParticleHPIsoData(), Init(), MakeElasticData(), and ThinOut().

theFileName

G4String G4ParticleHPIsoData::theFileName

private

Definition at line 114 of file G4ParticleHPIsoData.hh.

theFissionData

G4ParticleHPVector* G4ParticleHPIsoData::theFissionData

private

Definition at line 108 of file G4ParticleHPIsoData.hh.

Referenced by G4ParticleHPIsoData(), Init(), MakeFissionData(), and ThinOut().

theInelasticData

G4ParticleHPVector* G4ParticleHPIsoData::theInelasticData

private

Definition at line 111 of file G4ParticleHPIsoData.hh.

Referenced by G4ParticleHPIsoData(), Init(), MakeInelasticData(), and ThinOut().

theNames

G4ParticleHPNames G4ParticleHPIsoData::theNames

private

Definition at line 115 of file G4ParticleHPIsoData.hh.

Referenced by GetName(), and Init().

---

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

• source/processes/hadronic/models/particle_hp/include/G4ParticleHPIsoData.hh
• source/processes/hadronic/models/particle_hp/src/G4ParticleHPIsoData.cc