G4ParticleHPIsoData Class Reference

#include <G4ParticleHPIsoData.hh>

Public Member Functions
	G4ParticleHPIsoData ()
	~G4ParticleHPIsoData ()
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, G4ParticleDefinition *projectile, const char *dataDirVariable)
void	Init (G4int A, G4int Z, G4int M, G4double abun, G4ParticleDefinition *projectile, const char *dataDirVariable)
G4ParticleHPVector *	MakeElasticData ()
G4ParticleHPVector *	MakeFissionData ()
G4ParticleHPVector *	MakeCaptureData ()
G4ParticleHPVector *	MakeInelasticData ()
G4ParticleHPVector *	MakeChannelData ()
G4String	GetName (G4int A, G4int Z, G4String base, G4String rest)
void	FillChannelData (G4ParticleHPVector *aBuffer)
void	ThinOut (G4double precision)

Detailed Description

Definition at line 45 of file G4ParticleHPIsoData.hh.

Constructor & Destructor Documentation

G4ParticleHPIsoData::G4ParticleHPIsoData ( )

inline

Definition at line 49 of file G4ParticleHPIsoData.hh.

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

G4ParticleHPIsoData::~G4ParticleHPIsoData ( )

inline

Definition at line 58 of file G4ParticleHPIsoData.hh.

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

Member Function Documentation

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

Here is the call graph for this function:

Here is the caller graph for this function:

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

Here is the call graph for this function:

G4double G4ParticleHPIsoData::GetXsec ( G4double  energy )

inline

Definition at line 60 of file G4ParticleHPIsoData.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

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

Here is the call graph for this function:

Here is the caller graph for this function:

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(getenv("G4ParticleHPDebug")) G4cout << "G4ParticleHPIsoData::Init = "<< filename <<" "<< A << " " << Z <<G4endl;
#endif
    
    if(Z==1 && (aFile.GetZ()!=Z || std::abs(aFile.GetA()-A)>0.0001) )
    {
      if(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;
  }

Here is the call graph for this function:

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

Here is the call graph for this function:

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 ( getenv( dataDirVariable ) ) {
     baseName = getenv( dataDirVariable );
  } else {
     baseName = 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();

  }

Here is the call graph for this function:

G4ParticleHPVector* G4ParticleHPIsoData::MakeCaptureData ( )

inline

Definition at line 78 of file G4ParticleHPIsoData.hh.

    {return theCaptureData;}

Here is the caller graph for this function:

G4ParticleHPVector* G4ParticleHPIsoData::MakeChannelData ( )

inline

Definition at line 82 of file G4ParticleHPIsoData.hh.

    {return theChannelData;}

Here is the caller graph for this function:

G4ParticleHPVector* G4ParticleHPIsoData::MakeElasticData ( )

inline

Definition at line 74 of file G4ParticleHPIsoData.hh.

    {return theElasticData;}

Here is the caller graph for this function:

G4ParticleHPVector* G4ParticleHPIsoData::MakeFissionData ( )

inline

Definition at line 76 of file G4ParticleHPIsoData.hh.

    {return theFissionData;}

Here is the caller graph for this function:

G4ParticleHPVector* G4ParticleHPIsoData::MakeInelasticData ( )

inline

Definition at line 80 of file G4ParticleHPIsoData.hh.

    {return theInelasticData;}

Here is the caller graph for this function:

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

Here is the call graph for this function:

---

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

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