#include <G4AblaDataFile.hh>

Inheritance diagram for G4AblaDataFile:

Public Member Functions
	G4AblaDataFile ()

G4double	getAlpha (G4int A, G4int Z)

G4int	getAlphaCols ()

G4int	getAlphaRows ()

G4double	getBeta2 (G4int A, G4int Z)

G4double	getBeta4 (G4int A, G4int Z)

G4double	getEcnz (G4int A, G4int Z)

G4double	getMexp (G4int A, G4int Z)

G4int	getMexpID (G4int A, G4int Z)

G4double	getPace2 (G4int A, G4int Z)

G4int	getPaceCols ()

G4int	getPaceRows ()

G4double	getRms (G4int A, G4int Z)

G4double	getVgsld (G4int A, G4int Z)

bool	readData ()

G4bool	setAlpha (G4int A, G4int Z, G4double value)

G4bool	setBeta2 (G4int A, G4int Z, G4double value)

G4bool	setBeta4 (G4int A, G4int Z, G4double value)

G4bool	setEcnz (G4int A, G4int Z, G4double value)

G4bool	setMexp (G4int A, G4int Z, G4double value)

G4bool	setMexpID (G4int A, G4int Z, G4int value)

G4bool	setPace2 (G4int A, G4int Z, G4double value)

G4bool	setRms (G4int A, G4int Z, G4double value)

G4bool	setVgsld (G4int A, G4int Z, G4double value)

	~G4AblaDataFile ()

Private Attributes
G4double	alpha [alphaRows][alphaCols]

G4double	beta2 [betaRows][betaCols]

G4double	beta4 [betaRows][betaCols]

G4double	ecnz [alphaRows][alphaCols]

G4double	mexp [massRows][massCols]

G4int	mexpid [massRows][massCols]

G4double	pace2 [paceRows][paceCols]

G4double	rms [rmsRows][rmsCols]

G4int	verboseLevel

G4double	vgsld [alphaRows][alphaCols]

Static Private Attributes
static const G4int	alphaCols = 99

static const G4int	alphaRows = 154

static const G4int	betaCols = 137

static const G4int	betaRows = 251

static const G4int	massCols = 13

static const G4int	massRows = 154

static const G4int	paceCols = 500

static const G4int	paceRows = 500

static const G4int	rmsCols = 99

static const G4int	rmsRows = 154

Detailed Description

Read ABLA data from files.

Definition at line 45 of file G4AblaDataFile.hh.

Constructor & Destructor Documentation

◆ G4AblaDataFile()

G4AblaDataFile::G4AblaDataFile ( )

Definition at line 52 of file G4AblaDataFile.cc.

                               {
#else
G4AblaDataFile::G4AblaDataFile(G4INCL::Config *config)
  : G4AblaVirtualData(config) {
  theConfig = config;
#endif
  verboseLevel = 0;
}

References G4AblaDataFile(), G4AblaVirtualData::G4AblaVirtualData(), and verboseLevel.

Referenced by G4AblaDataFile().

◆ ~G4AblaDataFile()

G4AblaDataFile::~G4AblaDataFile ( )

Definition at line 61 of file G4AblaDataFile.cc.

62{

63}

Member Function Documentation

◆ getAlpha()

double G4AblaVirtualData::getAlpha	(	G4int	A,
		G4int	Z
	)

inherited

Get the value of Alpha.

Definition at line 111 of file G4AblaVirtualData.cc.

{
  return alpha[A][Z];
}

References A, G4AblaVirtualData::alpha, and Z.

Referenced by G4Abla::initEvapora().

◆ getAlphaCols()

int G4AblaVirtualData::getAlphaCols ( )

inherited

Definition at line 161 of file G4AblaVirtualData.cc.

{
  return alphaCols;
}

References G4AblaVirtualData::alphaCols.

◆ getAlphaRows()

int G4AblaVirtualData::getAlphaRows ( )

inherited

Definition at line 156 of file G4AblaVirtualData.cc.

{
  return alphaRows;
}

References G4AblaVirtualData::alphaRows.

◆ getBeta2()

double G4AblaVirtualData::getBeta2	(	G4int	A,
		G4int	Z
	)

inherited

Get the value of beta2 deformation.

Definition at line 146 of file G4AblaVirtualData.cc.

{
  return beta2[A][Z];
}

References A, G4AblaVirtualData::beta2, and Z.

Referenced by G4Abla::initEvapora().

◆ getBeta4()

double G4AblaVirtualData::getBeta4	(	G4int	A,
		G4int	Z
	)

inherited

Get the value of beta4 deformation.

Definition at line 151 of file G4AblaVirtualData.cc.

{
  return beta4[A][Z];
}

References A, G4AblaVirtualData::beta4, and Z.

Referenced by G4Abla::initEvapora().

◆ getEcnz()

double G4AblaVirtualData::getEcnz	(	G4int	A,
		G4int	Z
	)

inherited

Get the value of Ecnz.

Definition at line 116 of file G4AblaVirtualData.cc.

{
  return ecnz[A][Z];
}

References A, G4AblaVirtualData::ecnz, and Z.

Referenced by G4Abla::initEvapora().

◆ getMexp()

double G4AblaVirtualData::getMexp	(	G4int	A,
		G4int	Z
	)

inherited

Get the value of experimental masses.

Definition at line 136 of file G4AblaVirtualData.cc.

{
  return mexp[A][Z];
}

References A, G4AblaVirtualData::mexp, and Z.

Referenced by G4Abla::initEvapora().

◆ getMexpID()

int G4AblaVirtualData::getMexpID	(	G4int	A,
		G4int	Z
	)

inherited

Get the value of experimental masses ID.

Definition at line 141 of file G4AblaVirtualData.cc.

{
  return mexpid[A][Z];
}

References A, G4AblaVirtualData::mexpid, and Z.

Referenced by G4Abla::initEvapora().

◆ getPace2()

double G4AblaVirtualData::getPace2	(	G4int	A,
		G4int	Z
	)

inherited

Get the value of Pace2.

Definition at line 126 of file G4AblaVirtualData.cc.

{
  return pace2[A][Z];
}

References A, G4AblaVirtualData::pace2, and Z.

Referenced by G4Abla::initEvapora(), and readData().

◆ getPaceCols()

int G4AblaVirtualData::getPaceCols ( )

inherited

Definition at line 169 of file G4AblaVirtualData.cc.

{
  return paceCols;
}

References G4AblaVirtualData::paceCols.

Referenced by readData().

◆ getPaceRows()

int G4AblaVirtualData::getPaceRows ( )

inherited

Definition at line 165 of file G4AblaVirtualData.cc.

{
  return paceRows;
}

References G4AblaVirtualData::paceRows.

Referenced by readData().

◆ getRms()

double G4AblaVirtualData::getRms	(	G4int	A,
		G4int	Z
	)

inherited

Get the value of RMS.

Definition at line 131 of file G4AblaVirtualData.cc.

{
  return rms[A][Z];
}

References A, G4AblaVirtualData::rms, and Z.

Referenced by G4Abla::initEvapora().

◆ getVgsld()

double G4AblaVirtualData::getVgsld	(	G4int	A,
		G4int	Z
	)

inherited

Get the value of Vgsld.

Definition at line 121 of file G4AblaVirtualData.cc.

{
  return vgsld[A][Z];
}

References A, G4AblaVirtualData::vgsld, and Z.

Referenced by G4Abla::initEvapora().

◆ readData()

bool G4AblaDataFile::readData ( )

virtual

Read all data from files.

Implements G4AblaVirtualData.

Definition at line 68 of file G4AblaDataFile.cc.

{
#ifdef ABLAXX_IN_GEANT4_MODE
  if(!std::getenv("G4ABLADATA")) {
    //    throw G4HadronicException(__FILE__, __LINE__, "ERROR: Data
    //    missing. Set environment variable G4ABLA3.0 to point to the
    //    directory containing data files needed by INCL and ABLA
    //    models.");
    // G4String errorMessage1 = "ERROR: Data missing. Set environment variable G4ABLADATA\n";
    // G4String errorMessage2 = "\t to point to the directory containing data files needed\n";
    // G4String errorMessage3 = "\t by INCL and ABLA models.\n";
    // G4String errorMessage = errorMessage1 + errorMessage2 + errorMessage3;
    // G4Exception(errorMessage);
    G4ExceptionDescription ed;
    ed << " Data missing: set environment variable G4ABLADATA\n"
       << " to point to the directory containing data files needed\n"
       << " by the ABLA model" << G4endl;
    G4Exception("G4AblaDataFile::readData()","ABLA_001",
                FatalException, ed);
  }
  G4String dataPath(std::getenv("G4ABLADATA"));
#else
  G4String dataPath(theConfig->getABLAXXDataFilePath().c_str());
#endif
  G4String flAlphaFile(dataPath + "/flalpha.dat");
  G4String frldmFile(  dataPath + "/frldm.dat");
  G4String vgsldFile(  dataPath + "/vgsld.dat");
  G4String pace2File(  dataPath + "/pace2.dat");
  G4String rmsFile(  dataPath + "/rms.dat");
  G4String defoFile(  dataPath + "/defo.dat");
  G4String massFile(  dataPath + "/mass2003.dat");
 
  if(verboseLevel > 1) {
    // G4cout <<"Data path   = " << dataPath    << G4endl;
    // G4cout <<"FlAlphaFile = " << flAlphaFile << G4endl;
    // G4cout <<"FrldmFile   = " << frldmFile   << G4endl;
    // G4cout <<"VgsldFile   = " << vgsldFile   << G4endl;
    // G4cout <<"Pace2File   = " << pace2File   << G4endl;
  }
  
  std::ifstream flalphain(flAlphaFile.c_str());
  std::ifstream frldmin(frldmFile.c_str());  
  std::ifstream vgsldin(vgsldFile.c_str());  
  std::ifstream pace2in(pace2File.c_str());
  std::ifstream rmsin(rmsFile.c_str());
  std::ifstream defoin(defoFile.c_str());
  std::ifstream massin(massFile.c_str());
 
  std::filebuf *buf1 = flalphain.rdbuf();
  std::filebuf *buf2 = frldmin.rdbuf();
  std::filebuf *buf3 = vgsldin.rdbuf();
  std::filebuf *buf4 = pace2in.rdbuf();  
  std::filebuf *buf5 = rmsin.rdbuf(); 
  std::filebuf *buf6 = defoin.rdbuf();
  std::filebuf *buf7 = massin.rdbuf();
  if (!((buf1->is_open()) && (buf2->is_open()) && (buf3->is_open()) && (buf4->is_open()) && (buf5->is_open()) && (buf6->is_open()) && (buf7->is_open()))) {
#ifdef ABLAXX_IN_GEANT4_MODE
    G4ExceptionDescription ed;
    ed << "Data missing: could not find ABLA data file in " << dataPath
       << "defined by environment variable G4ABLADATA" << G4endl;
    G4Exception("G4AblaDataFile::readData()", "ABLA", FatalException, ed);
#else
    std::cerr << "Error opening file." << std::endl;
#endif
  }
  
  G4double fflalpha, ffrldm, fvgsld, fpace2, frms;
  int fj,fk,a2,a3,a4;
  G4double fbeta2,fbeta4;
  G4double a7;
  const G4int rows = 99;
  const G4int cols = 154;
  const G4int rowsbeta = 137;
  const G4int colsbeta = 251;
  const G4int rowsmass = 13;
  const G4int colsmass = 154;
  const G4int massnumbers = 263;
  for(int i = 0; i < rows; i++) {
    for(int j = 0; j < cols; j++) {
      setAlpha(j, i, 0.0);
      setEcnz( j, i, 0.0);
      setVgsld(j, i, 0.0);
      setRms(j, i, 0.0);
    }
  }
  
  for(int i = 0; i < rows; i++) {
    for(int j = 0; j < cols; j++) {
      flalphain >> fflalpha;
      frldmin >> ffrldm;
      vgsldin >> fvgsld;  
      rmsin >> frms;    
      setAlpha(j, i, fflalpha);
      setEcnz( j, i, ffrldm);
      setVgsld(j, i, fvgsld);
      setRms(j, i, frms);
    }
  }
 
 
  for(int i = 0; i < rowsbeta; i++) {
    for(int j = 0; j < colsbeta; j++) {
      setBeta2(j, i, 0.0);
      setBeta4(j, i, 0.0);
    }
  }
 
  for(int i = 0; i < 8983; i++) {
  defoin >> fj >> fk >> fbeta2 >> fbeta4;
  setBeta2(fk, fj, fbeta2);
  setBeta4(fk, fj, fbeta4);
  }
 
  for(int i = 0; i < rowsmass; i++) {
    for(int j = 0; j < colsmass; j++) {
      setMexp(j, i, 0.0);
      setMexpID(j,i,0);
    }
  }
   massin >> a2 >> a3 >> a4 >> a7 ;
   while(!massin.eof()){
   //
     if(a3<13.){
     setMexpID(a2,a3,1);
     setMexp(a2,a3,938.7829835*a3+939.5653301*a2-1.*a4*a7/1000.);
     }
     massin >> a2 >> a3 >> a4 >> a7 ;
   }
 
  flalphain.close();
  frldmin.close();  
  vgsldin.close();
  rmsin.close();
  defoin.close();
  massin.close();
 
  G4String str1, str2, str3;
  for(int i = 0; i < 500; i++) {
    for(int j = 0; j < 500; j++) {
      setPace2(i, j, 0.0);
    }
  }
  
  int A = 0, Zbegin = 0, Zend = 0;
  for(int i = 0; i < massnumbers; i++) {
    pace2in >> str1 >> A >> str2 >> Zbegin >> str3 >> Zend;
    if(Zbegin >= 0 && Zbegin < getPaceCols() &&
       A >= 0 && A < getPaceRows()) {
      for(int j = Zbegin; j <= Zend; j++) {
    pace2in >> fpace2;
    setPace2(A, j, fpace2);
      }
    }
  }
  pace2in.close();
  if(std::abs(getPace2(A, Zend) - 114516.10) > 1e-6) {
    std::cerr << "ERROR: Problem in parsing datafile " + pace2File << std::endl;
    return false;
  }
 
  return true;
}

References A, FatalException, G4endl, G4Exception(), G4AblaVirtualData::getPace2(), G4AblaVirtualData::getPaceCols(), G4AblaVirtualData::getPaceRows(), G4AblaVirtualData::setAlpha(), G4AblaVirtualData::setBeta2(), G4AblaVirtualData::setBeta4(), G4AblaVirtualData::setEcnz(), G4AblaVirtualData::setMexp(), G4AblaVirtualData::setMexpID(), G4AblaVirtualData::setPace2(), G4AblaVirtualData::setRms(), G4AblaVirtualData::setVgsld(), and verboseLevel.

Referenced by G4Abla::initEvapora().

◆ setAlpha()

bool G4AblaVirtualData::setAlpha	(	G4int	A,
		G4int	Z,
		G4double	value
	)

inherited

Set the value of Alpha.

Definition at line 47 of file G4AblaVirtualData.cc.

{
  alpha[A][Z] = value;
 
  return true;
}

References A, G4AblaVirtualData::alpha, and Z.

Referenced by readData().

◆ setBeta2()

bool G4AblaVirtualData::setBeta2	(	G4int	A,
		G4int	Z,
		G4double	value
	)

inherited

Set the value of beta2 deformation.

Definition at line 96 of file G4AblaVirtualData.cc.

{
  beta2[A][Z] = value;
 
  return true;
}

References A, G4AblaVirtualData::beta2, and Z.

Referenced by readData().

◆ setBeta4()

bool G4AblaVirtualData::setBeta4	(	G4int	A,
		G4int	Z,
		G4double	value
	)

inherited

Set the value of beta4 deformation.

Definition at line 103 of file G4AblaVirtualData.cc.

{
  beta4[A][Z] = value;
 
  return true;
}

References A, G4AblaVirtualData::beta4, and Z.

Referenced by readData().

◆ setEcnz()

bool G4AblaVirtualData::setEcnz	(	G4int	A,
		G4int	Z,
		G4double	value
	)

inherited

Set the value of Ecnz.

Definition at line 54 of file G4AblaVirtualData.cc.

{
  ecnz[A][Z] = value;
 
  return true;
}

References A, G4AblaVirtualData::ecnz, and Z.

Referenced by readData().

◆ setMexp()

bool G4AblaVirtualData::setMexp	(	G4int	A,
		G4int	Z,
		G4double	value
	)

inherited

Set the value of experimental masses.

Definition at line 82 of file G4AblaVirtualData.cc.

{
  mexp[A][Z] = value;
 
  return true;
}

References A, G4AblaVirtualData::mexp, and Z.

Referenced by readData().

◆ setMexpID()

bool G4AblaVirtualData::setMexpID	(	G4int	A,
		G4int	Z,
		G4int	value
	)

inherited

Set the value of experimental masses ID.

Definition at line 89 of file G4AblaVirtualData.cc.

{
  mexpid[A][Z] = value;
 
  return true;
}

References A, G4AblaVirtualData::mexpid, and Z.

Referenced by readData().

◆ setPace2()

bool G4AblaVirtualData::setPace2	(	G4int	A,
		G4int	Z,
		G4double	value
	)

inherited

Set the value of Pace2.

Definition at line 68 of file G4AblaVirtualData.cc.

{
  pace2[A][Z] = value;
 
  return true;
}

References A, G4AblaVirtualData::pace2, and Z.

Referenced by readData().

◆ setRms()

bool G4AblaVirtualData::setRms	(	G4int	A,
		G4int	Z,
		G4double	value
	)

inherited

Set the value of RMS.

Definition at line 75 of file G4AblaVirtualData.cc.

{
  rms[A][Z] = value;
 
  return true;
}

References A, G4AblaVirtualData::rms, and Z.

Referenced by readData().

◆ setVgsld()

bool G4AblaVirtualData::setVgsld	(	G4int	A,
		G4int	Z,
		G4double	value
	)

inherited

Set the value of Vgsld.

Definition at line 61 of file G4AblaVirtualData.cc.

{
  vgsld[A][Z] = value;
 
  return true;
}

References A, G4AblaVirtualData::vgsld, and Z.

Referenced by readData().

Field Documentation

◆ alpha

G4double G4AblaVirtualData::alpha[alphaRows][alphaCols]

privateinherited

Definition at line 185 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getAlpha(), and G4AblaVirtualData::setAlpha().

◆ alphaCols

const G4int G4AblaVirtualData::alphaCols = 99

staticprivateinherited

Definition at line 171 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getAlphaCols().

◆ alphaRows

const G4int G4AblaVirtualData::alphaRows = 154

staticprivateinherited

Definition at line 170 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getAlphaRows().

◆ beta2

G4double G4AblaVirtualData::beta2[betaRows][betaCols]

privateinherited

Definition at line 192 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getBeta2(), and G4AblaVirtualData::setBeta2().

◆ beta4

G4double G4AblaVirtualData::beta4[betaRows][betaCols]

privateinherited

Definition at line 193 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getBeta4(), and G4AblaVirtualData::setBeta4().

◆ betaCols

const G4int G4AblaVirtualData::betaCols = 137

staticprivateinherited

Definition at line 180 of file G4AblaVirtualData.hh.

◆ betaRows

const G4int G4AblaVirtualData::betaRows = 251

staticprivateinherited

Definition at line 179 of file G4AblaVirtualData.hh.

◆ ecnz

G4double G4AblaVirtualData::ecnz[alphaRows][alphaCols]

privateinherited

Definition at line 186 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getEcnz(), and G4AblaVirtualData::setEcnz().

◆ massCols

const G4int G4AblaVirtualData::massCols = 13

staticprivateinherited

Definition at line 183 of file G4AblaVirtualData.hh.

◆ massRows

const G4int G4AblaVirtualData::massRows = 154

staticprivateinherited

Definition at line 182 of file G4AblaVirtualData.hh.

◆ mexp

G4double G4AblaVirtualData::mexp[massRows][massCols]

privateinherited

Definition at line 190 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getMexp(), and G4AblaVirtualData::setMexp().

◆ mexpid

G4int G4AblaVirtualData::mexpid[massRows][massCols]

privateinherited

Definition at line 191 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getMexpID(), and G4AblaVirtualData::setMexpID().

◆ pace2

G4double G4AblaVirtualData::pace2[paceRows][paceCols]

privateinherited

Definition at line 188 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getPace2(), and G4AblaVirtualData::setPace2().

◆ paceCols

const G4int G4AblaVirtualData::paceCols = 500

staticprivateinherited

Definition at line 174 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getPaceCols().

◆ paceRows

const G4int G4AblaVirtualData::paceRows = 500

staticprivateinherited

Definition at line 173 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getPaceRows().

◆ rms

G4double G4AblaVirtualData::rms[rmsRows][rmsCols]

privateinherited

Definition at line 189 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getRms(), and G4AblaVirtualData::setRms().

◆ rmsCols

const G4int G4AblaVirtualData::rmsCols = 99

staticprivateinherited

Definition at line 177 of file G4AblaVirtualData.hh.

◆ rmsRows

const G4int G4AblaVirtualData::rmsRows = 154

staticprivateinherited

Definition at line 176 of file G4AblaVirtualData.hh.

◆ verboseLevel

G4int G4AblaDataFile::verboseLevel

private

Definition at line 61 of file G4AblaDataFile.hh.

Referenced by G4AblaDataFile(), and readData().

◆ vgsld

G4double G4AblaVirtualData::vgsld[alphaRows][alphaCols]

privateinherited

Definition at line 187 of file G4AblaVirtualData.hh.

Referenced by G4AblaVirtualData::getVgsld(), and G4AblaVirtualData::setVgsld().

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

source/processes/hadronic/models/abla/include/G4AblaDataFile.hh
source/processes/hadronic/models/abla/src/G4AblaDataFile.cc

Public Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4AblaDataFile()

◆ ~G4AblaDataFile()

Member Function Documentation

◆ getAlpha()

◆ getAlphaCols()

◆ getAlphaRows()

◆ getBeta2()

◆ getBeta4()

◆ getEcnz()

◆ getMexp()

◆ getMexpID()

◆ getPace2()

◆ getPaceCols()

◆ getPaceRows()

◆ getRms()

◆ getVgsld()

◆ readData()

◆ setAlpha()

◆ setBeta2()

◆ setBeta4()

◆ setEcnz()

◆ setMexp()

◆ setMexpID()

◆ setPace2()

◆ setRms()

◆ setVgsld()

Field Documentation

◆ alpha

◆ alphaCols

◆ alphaRows

◆ beta2

◆ beta4

◆ betaCols

◆ betaRows

◆ ecnz

◆ massCols

◆ massRows

◆ mexp

◆ mexpid

◆ pace2

◆ paceCols

◆ paceRows

◆ rms

◆ rmsCols

◆ rmsRows

◆ verboseLevel

◆ vgsld