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::G4ParticleHPIsoData ( )

inline

Definition at line 49 of file G4ParticleHPIsoData.hh.

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

~G4ParticleHPIsoData()

G4ParticleHPIsoData::~G4ParticleHPIsoData ( )

inline

Definition at line 58 of file G4ParticleHPIsoData.hh.

{if(theChannelData!=0) delete 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();
  }

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

GetXsec()

G4double G4ParticleHPIsoData::GetXsec ( G4double  energy )

inline

Definition at line 60 of file G4ParticleHPIsoData.hh.

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

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

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

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

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

MakeCaptureData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeCaptureData ( )

inline

Definition at line 78 of file G4ParticleHPIsoData.hh.

    {return theCaptureData;}

Referenced by G4ParticleHPElementData::UpdateData().

MakeChannelData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeChannelData ( )

inline

Definition at line 82 of file G4ParticleHPIsoData.hh.

    {return theChannelData;}

Referenced by G4ParticleHPChannel::UpdateData().

MakeElasticData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeElasticData ( )

inline

Definition at line 74 of file G4ParticleHPIsoData.hh.

    {return theElasticData;}

Referenced by G4ParticleHPElementData::UpdateData().

MakeFissionData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeFissionData ( )

inline

Definition at line 76 of file G4ParticleHPIsoData.hh.

    {return theFissionData;}

Referenced by G4ParticleHPElementData::UpdateData().

MakeInelasticData()

G4ParticleHPVector * G4ParticleHPIsoData::MakeInelasticData ( )

inline

Definition at line 80 of file G4ParticleHPIsoData.hh.

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

---

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

• G4ParticleHPIsoData.hh
• G4ParticleHPIsoData.cc