G4AblaDataFile Class Reference

#include <G4AblaDataFile.hh>

Inheritance diagram for G4AblaDataFile:

Public Member Functions
	G4AblaDataFile ()
virtual	~G4AblaDataFile ()=default
G4bool	readData ()
Public Member Functions inherited from G4AblaVirtualData
G4bool	setAlpha (G4int A, G4int Z, G4double value)
G4bool	setEcnz (G4int A, G4int Z, G4double value)
G4bool	setVgsld (G4int A, G4int Z, G4double value)
G4bool	setRms (G4int A, G4int Z, G4double value)
G4bool	setMexp (G4int A, G4int Z, G4double value)
G4bool	setMexpID (G4int A, G4int Z, G4int value)
G4bool	setBeta2 (G4int A, G4int Z, G4double value)
G4bool	setBeta4 (G4int A, G4int Z, G4double value)
G4double	getAlpha (G4int A, G4int Z)
G4double	getEcnz (G4int A, G4int Z)
G4double	getVgsld (G4int A, G4int Z)
G4double	getRms (G4int A, G4int Z)
G4double	getMexp (G4int A, G4int Z)
G4int	getMexpID (G4int A, G4int Z)
G4double	getBeta2 (G4int A, G4int Z)
G4double	getBeta4 (G4int A, G4int Z)

Additional Inherited Members
Protected Member Functions inherited from G4AblaVirtualData
	G4AblaVirtualData ()
virtual	~G4AblaVirtualData ()=default

Detailed Description

Read ABLA data from files.

Definition at line 42 of file G4AblaDataFile.hh.

Constructor & Destructor Documentation

G4AblaDataFile()

G4AblaDataFile::G4AblaDataFile

(

)

Definition at line 44 of file G4AblaDataFile.cc.

{
  verboseLevel = 0;
}

~G4AblaDataFile()

virtual G4AblaDataFile::~G4AblaDataFile ( )

virtualdefault

Member Function Documentation

readData()

G4bool G4AblaDataFile::readData ( )

virtual

Read all data from files.

Implements G4AblaVirtualData.

Definition at line 52 of file G4AblaDataFile.cc.

{
  if (!G4FindDataDir("G4ABLADATA")) {
    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(G4FindDataDir("G4ABLADATA"));
 
  G4String flAlphaFile(dataPath + "/flalpha.dat");
  G4String frldmFile(dataPath + "/frldm.dat");
  G4String vgsldFile(dataPath + "/vgsld.dat");
  G4String rmsFile(dataPath + "/rms.dat");
  G4String defoFile(dataPath + "/defo.dat");
  G4String massFile(dataPath + "/mass2020.dat");
 
  if (verboseLevel > 1) {
    // G4cout <<"Data path   = " << dataPath    << G4endl;
    // G4cout <<"FlAlphaFile = " << flAlphaFile << G4endl;
    // G4cout <<"FrldmFile   = " << frldmFile   << G4endl;
    // G4cout <<"VgsldFile   = " << vgsldFile   << G4endl;
  }
 
  std::ifstream flalphain(flAlphaFile.c_str());
  std::ifstream frldmin(frldmFile.c_str());
  std::ifstream vgsldin(vgsldFile.c_str());
  std::ifstream rmsin(rmsFile.c_str());
  std::ifstream defoin(defoFile.c_str());
  std::ifstream massin(massFile.c_str());
 
  if (!massin.is_open()) {
    massFile = dataPath + "/mass2016.dat";
    massin.close();
    massin.open(massFile.c_str());
    std::cout << "Mass evaluation file mass2020.dat not found, current file: " << massFile.c_str()
              << std::endl;
 
    if (!massin.is_open()) {
      massFile = dataPath + "/mass2003.dat";
      massin.close();
      massin.open(massFile.c_str());
      std::cout << "Mass evaluation file mass2016.dat not found, current file: " << massFile.c_str()
                << std::endl;
    }
  }
 
  std::filebuf* buf1 = flalphain.rdbuf();
  std::filebuf* buf2 = frldmin.rdbuf();
  std::filebuf* buf3 = vgsldin.rdbuf();
  std::filebuf* buf4 = rmsin.rdbuf();
  std::filebuf* buf5 = defoin.rdbuf();
  std::filebuf* buf6 = massin.rdbuf();
  if (!((buf1->is_open()) && (buf2->is_open()) && (buf3->is_open()) && (buf4->is_open())
        && (buf5->is_open()) && (buf6->is_open())))
  {
    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);
  }
 
  G4double fflalpha, ffrldm, fvgsld, frms;
  G4int fj = 0, fk = 0, a2, a3, a4;
  G4double fbeta2, fbeta4;
  G4double a7;
  const G4int rows = 99;
  const G4int cols = 154;
  const G4int rowsbeta = 137;
  const G4int colsbeta = 251;
 
  for (G4int i = 0; i < zcols; i++) {
    for (G4int j = 0; j < nrows; j++) {
      setAlpha(j, i, 0.0);
      setEcnz(j, i, 0.0);
      setVgsld(j, i, 0.0);
      setRms(j, i, 0.0);
    }
  }
 
  for (G4int i = 0; i < rows; i++) {
    for (G4int 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 (G4int i = 0; i < rowsbeta; i++) {
    for (G4int j = 0; j < colsbeta; j++) {
      setBeta2(j, i, 0.0);
      setBeta4(j, i, 0.0);
    }
  }
 
  defoin >> fj >> fk >> fbeta2 >> fbeta4;
  while (!defoin.eof()) {
    setBeta2(fk, fj, fbeta2);
    setBeta4(fk, fj, fbeta4);
    defoin >> fj >> fk >> fbeta2 >> fbeta4;
  }
 
  for (G4int i = 0; i < zcols; i++) {
    for (G4int j = 0; j < nrows; j++) {
      setMexp(j, i, 0.0);
      setMexpID(j, i, 0);
    }
  }
  massin >> a2 >> a3 >> a4 >> a7;
  while (!massin.eof()) {
    //
    if (a3 < lpcols) {
      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();
 
  return true;
}

---

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

• G4AblaDataFile.hh
• G4AblaDataFile.cc