XrayFluoDataSet Class Reference

#include <XrayFluoDataSet.hh>

Inheritance diagram for XrayFluoDataSet:

Collaboration diagram for XrayFluoDataSet:

Public Member Functions
	XrayFluoDataSet (G4int Z, G4DataVector *points, G4DataVector *values, const G4VDataSetAlgorithm *interpolation, G4double unitE=CLHEP::MeV, G4double unitData=CLHEP::barn)
	XrayFluoDataSet (G4int Z, const G4String &dataFile, const G4VDataSetAlgorithm *interpolation, G4double unitE=CLHEP::MeV, G4double unitData=CLHEP::barn)
	~XrayFluoDataSet ()
G4double	FindValue (G4double e, G4int) const
virtual const G4VEMDataSet *	GetComponent (G4int) const
virtual void	AddComponent (G4VEMDataSet *)
virtual size_t	NumberOfComponents () const
virtual void	SetEnergiesData (G4DataVector *, G4DataVector *, int=0)
virtual G4bool	SaveData (const G4String &) const
virtual G4double	RandomSelect (int=0) const
void	PrintData () const
const G4DataVector &	GetEnergies (G4int) const
const G4DataVector &	GetData (G4int) const
virtual const G4DataVector &	GetLogEnergies (G4int) const
virtual const G4DataVector &	GetLogData (G4int) const
virtual void	SetLogEnergiesData (G4DataVector *, G4DataVector *, G4DataVector *, G4DataVector *, G4int)
virtual G4bool	LoadNonLogData (const G4String &)
Public Member Functions inherited from G4VEMDataSet
	G4VEMDataSet ()
virtual	~G4VEMDataSet ()

Private Member Functions
G4bool	LoadData (const G4String &dataFile)
G4int	FindBinLocation (G4double energy) const

Private Attributes
G4DataVector *	energies
G4DataVector *	data
const G4VDataSetAlgorithm *	algorithm
G4double	unit1
G4double	unit2
size_t	numberOfBins

Detailed Description

Definition at line 51 of file XrayFluoDataSet.hh.

Constructor & Destructor Documentation

XrayFluoDataSet() [1/2]

XrayFluoDataSet::XrayFluoDataSet	(	G4int	Z,
		G4DataVector *	points,
		G4DataVector *	values,
		const G4VDataSetAlgorithm *	interpolation,
		G4double	unitE = CLHEP::MeV,
		G4double	unitData = CLHEP::barn
	)

Definition at line 42 of file XrayFluoDataSet.cc.

  :energies(points), data(values), algorithm(interpolation)
{
  numberOfBins = energies->size();
  unit1 = unitE;
  unit2 = unitData;
}

XrayFluoDataSet() [2/2]

XrayFluoDataSet::XrayFluoDataSet	(	G4int	Z,
		const G4String &	dataFile,
		const G4VDataSetAlgorithm *	interpolation,
		G4double	unitE = CLHEP::MeV,
		G4double	unitData = CLHEP::barn
	)

Definition at line 56 of file XrayFluoDataSet.cc.

  : algorithm(interpolation)
{
  energies = new G4DataVector;
  data = new G4DataVector;
  unit1 = unitE;
  unit2 = unitData;  
  LoadData(dataFile);
  numberOfBins = energies->size();
}

Here is the call graph for this function:

~XrayFluoDataSet()

XrayFluoDataSet::~XrayFluoDataSet

(

)

Definition at line 73 of file XrayFluoDataSet.cc.

{ 
  delete energies;
  delete data;
}

Member Function Documentation

AddComponent()

virtual void XrayFluoDataSet::AddComponent ( G4VEMDataSet *  )

inlinevirtual

Implements G4VEMDataSet.

Definition at line 74 of file XrayFluoDataSet.hh.

{ }

FindBinLocation()

G4int XrayFluoDataSet::FindBinLocation ( G4double  energy ) const

private

Definition at line 107 of file XrayFluoDataSet.cc.

{
  // Protection against call outside allowed range
  G4double e0 = (*energies)[0];
  if (energy < e0)
    {
    
      energy = e0;
    }

  size_t lowerBound = 0;
  size_t upperBound = numberOfBins - 1;
  
  // Binary search
  while (lowerBound <= upperBound) 
    {
      size_t midBin = (lowerBound + upperBound)/2;
      if ( energy < (*energies)[midBin] ) upperBound = midBin-1;
      else lowerBound = midBin+1;
  }
 
  return upperBound;
}

Here is the caller graph for this function:

FindValue()

G4double XrayFluoDataSet::FindValue ( G4double  e, G4int    ) const

virtual

Implements G4VEMDataSet.

Definition at line 81 of file XrayFluoDataSet.cc.

{
  G4double value;
  G4double e0 = (*energies)[0];
  // Protections
  size_t bin = FindBinLocation(e);
  if (bin == numberOfBins)
    {
     
      value = (*data)[bin];
    }
  else if (e <= e0)
    {
    
      value = (*data)[0];
    }
  else
    {
      value = algorithm->Calculate(e,bin,*energies,*data);
    }
  
  return value;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetComponent()

virtual const G4VEMDataSet* XrayFluoDataSet::GetComponent ( G4int  ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 72 of file XrayFluoDataSet.hh.

{ return 0;} 

GetData()

const G4DataVector& XrayFluoDataSet::GetData ( G4int  ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 88 of file XrayFluoDataSet.hh.

{ return *data; }

GetEnergies()

const G4DataVector& XrayFluoDataSet::GetEnergies ( G4int  ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 87 of file XrayFluoDataSet.hh.

{ return *energies; }

GetLogData()

virtual const G4DataVector& XrayFluoDataSet::GetLogData ( G4int  ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 91 of file XrayFluoDataSet.hh.

{ return *data; }

GetLogEnergies()

virtual const G4DataVector& XrayFluoDataSet::GetLogEnergies ( G4int  ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 90 of file XrayFluoDataSet.hh.

{ return *energies; }

LoadData()

G4bool XrayFluoDataSet::LoadData ( const G4String &  dataFile )

privatevirtual

Implements G4VEMDataSet.

Definition at line 133 of file XrayFluoDataSet.cc.

{
  // Build the complete string identifying the file with the data set
  G4String dirFile = "";
  
  char* path;

  path = getenv("XRAYDATA");
  if (!path)
    path = getenv("PWD");

  //G4cout << path << G4endl;
  //G4cout << fileName << G4endl;


  G4String pathString(path);
  pathString += "\0";
  dirFile = pathString + "/" + fileName + ".dat";
    
  std::ifstream file(dirFile);
  std::filebuf* lsdp = file.rdbuf();
  
  if (! (lsdp->is_open()) )
    {
      G4ExceptionDescription execp;
      execp << "XrayFluoDataSet - data file: " + dirFile + " not found"<<G4endl;
      G4Exception("XrayFluoDataSet::LoadData()","example-xray_fluorescence01",
      FatalException, execp);
    }
  G4double a = 0;
  G4int k = 1;

  do
    {
      file >> a;
      G4int nColumns = 2;
      // The file is organized into two columns:
      // 1st column is the energy
      // 2nd column is the corresponding value
      // The file terminates with the pattern: -1   -1
      //                                       -2   -2
      if (a == -1 || a == -2)
    {
     
    }
      else
    {
      if (k%nColumns != 0)
        {   
          G4double e = a * unit1;
          energies->push_back(e);
        
          k++;

        }
      else if (k%nColumns == 0)
        {
          G4double value = a * unit2;
          data->push_back(value);
         
          k = 1;
        }
    }
      
    } while (a != -2); // end of file
  
  file.close();
  return true;

}

Here is the call graph for this function:

Here is the caller graph for this function:

LoadNonLogData()

virtual G4bool XrayFluoDataSet::LoadNonLogData ( const G4String &  )

inlinevirtual

Implements G4VEMDataSet.

Definition at line 100 of file XrayFluoDataSet.hh.

{return false;}

Here is the call graph for this function:

NumberOfComponents()

virtual size_t XrayFluoDataSet::NumberOfComponents ( void  ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 76 of file XrayFluoDataSet.hh.

{ return 0; }

PrintData()

void XrayFluoDataSet::PrintData ( void  ) const

virtual

Implements G4VEMDataSet.

Definition at line 206 of file XrayFluoDataSet.cc.

{
  size_t size = numberOfBins;
  for (size_t i=0; i<size; i++)
    {
      G4double e = (*energies)[i]  / unit1;
      G4double sigma = (*data)[i] / unit2 ;
      G4cout << "Point: "
         << e
         << " - Data value : "
         << sigma 
         << G4endl; 
    }
}

Here is the caller graph for this function:

RandomSelect()

virtual G4double XrayFluoDataSet::RandomSelect ( int  = 0 ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 83 of file XrayFluoDataSet.hh.

{return 0;}

Here is the call graph for this function:

SaveData()

virtual G4bool XrayFluoDataSet::SaveData ( const G4String &  ) const

inlinevirtual

Implements G4VEMDataSet.

Definition at line 81 of file XrayFluoDataSet.hh.

{return true;}

SetEnergiesData()

virtual void XrayFluoDataSet::SetEnergiesData ( G4DataVector *  , G4DataVector *  , int  = 0  )

inlinevirtual

Implements G4VEMDataSet.

Definition at line 78 of file XrayFluoDataSet.hh.

         {}

SetLogEnergiesData()

virtual void XrayFluoDataSet::SetLogEnergiesData ( G4DataVector *  , G4DataVector *  , G4DataVector *  , G4DataVector *  , G4int    )

inlinevirtual

Implements G4VEMDataSet.

Definition at line 94 of file XrayFluoDataSet.hh.

                                          {;}

Member Data Documentation

algorithm

const G4VDataSetAlgorithm* XrayFluoDataSet::algorithm

private

Definition at line 116 of file XrayFluoDataSet.hh.

data

G4DataVector* XrayFluoDataSet::data

private

Definition at line 114 of file XrayFluoDataSet.hh.

energies

G4DataVector* XrayFluoDataSet::energies

private

Definition at line 113 of file XrayFluoDataSet.hh.

numberOfBins

size_t XrayFluoDataSet::numberOfBins

private

Definition at line 121 of file XrayFluoDataSet.hh.

unit1

G4double XrayFluoDataSet::unit1

private

Definition at line 118 of file XrayFluoDataSet.hh.

unit2

G4double XrayFluoDataSet::unit2

private

Definition at line 119 of file XrayFluoDataSet.hh.

---

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

• XrayFluoDataSet.hh
• XrayFluoDataSet.cc