XrayFluoSiLiDetectorType Class Reference

#include <XrayFluoSiLiDetectorType.hh>

Inheritance diagram for XrayFluoSiLiDetectorType:

Collaboration diagram for XrayFluoSiLiDetectorType:

Public Member Functions
	~XrayFluoSiLiDetectorType ()
G4String	GetDetectorMaterial ()
G4double	ResponseFunction (G4double)
G4double	GetInfData (G4double, G4double, G4int)
G4double	GetSupData (G4double, G4double, G4int)
void	LoadResponseData (G4String)
void	LoadEfficiencyData (G4String)
Public Member Functions inherited from XrayFluoVDetectorType
virtual	~XrayFluoVDetectorType ()

Static Public Member Functions
static XrayFluoSiLiDetectorType *	GetInstance ()

Additional Inherited Members
Protected Member Functions inherited from XrayFluoVDetectorType
	XrayFluoVDetectorType ()

Detailed Description

Definition at line 50 of file XrayFluoSiLiDetectorType.hh.

Constructor & Destructor Documentation

XrayFluoSiLiDetectorType::~XrayFluoSiLiDetectorType

(

)

Definition at line 59 of file XrayFluoSiLiDetectorType.cc.

{
  std::map<G4int,G4DataVector*,std::less<G4int> >::iterator pos;

  for (pos = energyMap.begin(); pos != energyMap.end(); pos++)
    {
      G4DataVector* dataSet = (*pos).second;
      delete dataSet;
      dataSet = 0;
    }
  for (pos = dataMap.begin(); pos != dataMap.end(); pos++)
    {
      G4DataVector* dataSet = (*pos).second;
      delete dataSet;
      dataSet = 0;
    }

  delete interpolation4; 

}

Member Function Documentation

G4String XrayFluoSiLiDetectorType::GetDetectorMaterial ( )

virtual

Implements XrayFluoVDetectorType.

Definition at line 79 of file XrayFluoSiLiDetectorType.cc.

{
  return detectorMaterial;
}

G4double XrayFluoSiLiDetectorType::GetInfData ( G4double  , G4double  random, G4int  posIndex  )

virtual

Implements XrayFluoVDetectorType.

Definition at line 202 of file XrayFluoSiLiDetectorType.cc.

{
  G4double value = 0.;
  G4int zMin = 1;
  G4int zMax = 6; 
  
  G4int Z = posIndex;
  
  if (Z<zMin) {Z=zMin;}
  if (Z>zMax) {Z=zMax;}
  
  if (Z >= zMin && Z <= zMax)
    {
      std::map<G4int,G4DataVector*,std::less<G4int> >::const_iterator pos;
      pos = energyMap.find(Z);
      std::map<G4int,G4DataVector*,std::less<G4int> >::const_iterator posData;
      posData = dataMap.find(Z);


      if (pos!= energyMap.end())
    {
      G4DataVector energySet = *((*pos).second);
      G4DataVector dataSet = *((*posData).second);
      G4int nData = energySet.size();


      G4double dataSum = 0;
      G4double partSum = 0;
      G4int index = 0;

      // if data is not perfectly normalized (it may happen) 
      // rnadom number is renormalized, in case it is higer 
      //than the sum of all energies => segmentation fault.
      
      for  (G4int i = 0; i<nData; i++){
        dataSum += dataSet[i];
      }

      G4double normRandom = random*dataSum;


      while (normRandom> partSum)
        {
          
          partSum += dataSet[index];
          index++;
        }
     
      
      if (index >= 0 && index < nData)
        {
          value = energySet[index];
          
        }
      
    }
    }
  return value;
}

Here is the caller graph for this function:

XrayFluoSiLiDetectorType * XrayFluoSiLiDetectorType::GetInstance ( void  )

static

Definition at line 86 of file XrayFluoSiLiDetectorType.cc.

{
  if (instance == 0)
    {
      instance = new XrayFluoSiLiDetectorType;
     
    }
  return instance;
}

Here is the caller graph for this function:

G4double XrayFluoSiLiDetectorType::GetSupData ( G4double  , G4double  random, G4int  posIndex  )

virtual

Implements XrayFluoVDetectorType.

Definition at line 261 of file XrayFluoSiLiDetectorType.cc.

{
  G4double value = 0.;
  G4int zMin = 1;
  G4int zMax = 6;
  G4int Z = (posIndex+1);
  
  if (Z<zMin) {Z=zMin;}
  if (Z>zMax) {Z=zMax;}
  if (Z >= zMin && Z <= zMax)
    {
      std::map<G4int,G4DataVector*,std::less<G4int> >::const_iterator pos;
      pos = energyMap.find(Z);
      std::map<G4int,G4DataVector*,std::less<G4int> >::const_iterator posData;
      posData = dataMap.find(Z);
      if (pos!= energyMap.end())
    {
      G4DataVector energySet = *((*pos).second);
      G4DataVector dataSet = *((*posData).second);
      G4int nData = energySet.size();
      G4double dataSum = 0;
      G4double partSum = 0;
      G4int index = 0;

      // if data is not perfectly normalized (it may happen) 
      // rnadom number is renormalized, in case it is higer 
      //than the sum of all energies => segmentation fault.

      for  (G4int i = 0; i<nData; i++){
        dataSum += dataSet[i];
      }   

      G4double normRandom = random*dataSum;


      while (normRandom> partSum)
        {
          partSum += dataSet[index];
          index++;
        }
      
      
      if (index >= 0 && index < nData)
        {
          value = energySet[index];
        }
    }
    }
  return value;
}

Here is the caller graph for this function:

void XrayFluoSiLiDetectorType::LoadEfficiencyData ( G4String  fileName )

virtual

Implements XrayFluoVDetectorType.

Definition at line 400 of file XrayFluoSiLiDetectorType.cc.

{
  char* path = getenv("XRAYDATA");
  G4String dirFile;
  if (path) {
    G4String pathString(path);
    dirFile = pathString + "/" + fileName;
  }
  else{ 
    path = getenv("PWD");
    G4String pathString(path);
    dirFile = pathString + "/" + fileName;
  }

  interpolation4 = new G4LogLogInterpolation();
  efficiencySet = new XrayFluoDataSet(1,fileName,interpolation4,keV,1);
}

void XrayFluoSiLiDetectorType::LoadResponseData ( G4String  fileName )

virtual

Implements XrayFluoVDetectorType.

Definition at line 311 of file XrayFluoSiLiDetectorType.cc.

{
  std::ostringstream ost;
  ost << fileName<<".dat";
  G4String name = ost.str();
  
  char* path = getenv("XRAYDATA");
  G4String dirFile;
  if (path) {
    G4String pathString(path);
    pathString += "\0";
    dirFile = pathString + "/" + name;
  }
  else{ 
    path = getenv("PWD");
    G4String pathString(path);
    pathString += "\0";
    dirFile = pathString + "/" + name;
  }
  
  std::ifstream file(dirFile);
  std::filebuf* lsdp = file.rdbuf();
  
  if (! (lsdp->is_open()) )
    {
      G4ExceptionDescription execp;
      execp <<  "XrayFluoSiLiDetectorType - data file: " + dirFile + " not found";
      G4Exception("XrayFluoSiLiDetectorType::LoadResponseData()","example-xray_fluorescence07",
      FatalException, execp);
    }
  G4double a = 0;
  G4int k = 1;
  G4int q = 0;
  
  G4int Z = 1;
  G4DataVector* energies = new G4DataVector;
  G4DataVector* data = new G4DataVector;

  do
    {
      file >> a;
      G4int nColumns = 2;
      if (a == -1)
    {
      if (q == 0)
        {
          // End of a  data set
          energyMap[Z] = energies;
          dataMap[Z] = data;
          // Start of new shell data set
          energies = new G4DataVector;
          data = new G4DataVector;
          Z++;      
        }      
      q++;
      if (q == nColumns)
        {
          q = 0;
        }
    }
      else if (a == -2)
    {
      // End of file; delete the empty vectors 
      //created when encountering the last -1 -1 row
      delete energies;
      delete data;
      
    }
      else
    {
      // 1st column is energy
      if(k%nColumns != 0)
        {   
          G4double e = a * keV;
          energies->push_back(e);
          k++;
        }
      else if (k%nColumns == 0)
        {
          // 2nd column is data
          
          data->push_back(a);
          k = 1;
        }
    }
    } while (a != -2); // end of file
  file.close();    
}

Here is the call graph for this function:

G4double XrayFluoSiLiDetectorType::ResponseFunction ( G4double  energy )

virtual

Implements XrayFluoVDetectorType.

Definition at line 98 of file XrayFluoSiLiDetectorType.cc.

{
 
  G4double eMin = 1.500 *keV;
  G4double eMax = 6.403 *keV; 
  G4double value = 0.;
  G4double efficiency = 1.;
  
  const XrayFluoDataSet* dataSet = efficiencySet;
  G4int id = 0;
  G4DataVector energyVector;
  energyVector.push_back(1.486* keV);
  energyVector.push_back(1.740* keV);
  energyVector.push_back(3.688* keV);
  energyVector.push_back(4.510* keV);
  energyVector.push_back(5.414* keV);
  energyVector.push_back(6.404* keV);

  G4double infEnergy = 0 *keV;
  G4double supEnergy = 10* keV;
  G4int energyNumber = 0;
  G4double random = G4UniformRand();
 
  if (energy>=eMin && energy <=eMax)
    {

      for (G4int i=0; i<(G4int)energyVector.size(); i++){
    if (energyVector[i]/keV < energy/keV){

      infEnergy = energyVector[i];
      supEnergy = energyVector[i+1];

      energyNumber = i+1;

    }
      }
      

      G4double infData = GetInfData(energy, random, energyNumber);

      G4double supData = GetSupData(energy,random, energyNumber);

      value = (std::log10(infData)*std::log10(supEnergy/energy) +
           std::log10(supData)*std::log10(energy/infEnergy)) / 
    std::log10(supEnergy/infEnergy);
      value = std::pow(10,value);


    }
//    else if (energy<eMin || energy>eMax)
//      { 
//        G4double infEnergy = eMin;
//        G4double supEnergy = eMin/keV +1*keV;
 
//        G4double infData = GetInfData(eMin, random);
//        G4double supData = GetSupData(eMin,random);
//        value = (std::log10(infData)*std::log10(supEnergy/eMin) +
//             std::log10(supData)*std::log10(eMin/infEnergy)) / 
//      std::log10(supEnergy/infEnergy);
//        value = std::pow(10,value);
//        value = value-eMin+ energy;


//      }


  else if (energy > eMax)
    {

      energyNumber = 5;

      value = (GetSupData(energy, random, energyNumber))+(energy - 6.404* keV);

    }



  else 
    {
 
      energyNumber = 1;

      value = (GetInfData(energy, random, energyNumber))+(energy - 1.486* keV);


      //      G4double mean = -14.03 * eV + 1.0047*energy/eV;
      //      G4double stdDev = 35.38 * eV + 0.004385*energy/eV;
      
      
      //      value = (G4RandGauss::shoot(mean,stdDev))*eV;

    }
  G4double  RandomNum = G4UniformRand(); 
  
  efficiency = dataSet->FindValue(value,id);
  if ( RandomNum>efficiency )
    {
      value = 0.;
    }

  //  G4cout << value << G4endl; 
  return value;

}

Here is the call graph for this function:

---

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

• source/geant4.10.03.p02/examples/advanced/xray_fluorescence/include/XrayFluoSiLiDetectorType.hh
• source/geant4.10.03.p02/examples/advanced/xray_fluorescence/src/XrayFluoSiLiDetectorType.cc