XrayFluoHPGeDetectorType Class Reference

#include <XrayFluoHPGeDetectorType.hh>

Inheritance diagram for XrayFluoHPGeDetectorType:

Public Member Functions
	~XrayFluoHPGeDetectorType ()
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 XrayFluoHPGeDetectorType *	GetInstance ()

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

Detailed Description

Definition at line 50 of file XrayFluoHPGeDetectorType.hh.

Constructor & Destructor Documentation

XrayFluoHPGeDetectorType::~XrayFluoHPGeDetectorType

(

)

Definition at line 60 of file XrayFluoHPGeDetectorType.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 XrayFluoHPGeDetectorType::GetDetectorMaterial ( )

virtual

Implements XrayFluoVDetectorType.

Definition at line 80 of file XrayFluoHPGeDetectorType.cc.

{
  return detectorMaterial;
}

G4double XrayFluoHPGeDetectorType::GetInfData ( G4double  energy, G4double  random, G4int    )

virtual

Implements XrayFluoVDetectorType.

Definition at line 165 of file XrayFluoHPGeDetectorType.cc.

References python.hepunit::keV.

Referenced by ResponseFunction().

{
  G4double value = 0.;
  G4int zMin = 1;
  G4int zMax = 10; 
  
  G4int Z = ((G4int)(energy/keV));
  
  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 partSum = 0;
      G4int index = 0;
      
      while (random> partSum)
        {
          partSum += dataSet[index];
          index++;
        }
      
      
      if (index >= 0 && index < nData)
        {
          value = energySet[index];
          
        }
      
    }
    }
  return value;

}

XrayFluoHPGeDetectorType * XrayFluoHPGeDetectorType::GetInstance ( void  )

static

Definition at line 87 of file XrayFluoHPGeDetectorType.cc.

Referenced by XrayFluoPlaneDetectorConstruction::SetDetectorType(), XrayFluoMercuryDetectorConstruction::SetDetectorType(), and XrayFluoDetectorConstruction::SetDetectorType().

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

G4double XrayFluoHPGeDetectorType::GetSupData ( G4double  energy, G4double  random, G4int    )

virtual

Implements XrayFluoVDetectorType.

Definition at line 209 of file XrayFluoHPGeDetectorType.cc.

References python.hepunit::keV.

Referenced by ResponseFunction().

{
  G4double value = 0.;
  G4int zMin = 1;
  G4int zMax = 10;
  G4int Z = ((G4int)(energy/keV)+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 partSum = 0;
      G4int index = 0;
      
      while (random> partSum)
        {
          partSum += dataSet[index];
          index++;
        }
      
      
      if (index >= 0 && index < nData)
        {
          value = energySet[index];
        }
    }
    }
  return value;

}

void XrayFluoHPGeDetectorType::LoadEfficiencyData ( G4String  fileName )

virtual

Implements XrayFluoVDetectorType.

Definition at line 329 of file XrayFluoHPGeDetectorType.cc.

References python.hepunit::keV.

{
  /*
  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 XrayFluoHPGeDetectorType::LoadResponseData ( G4String  fileName )

virtual

Implements XrayFluoVDetectorType.

Definition at line 248 of file XrayFluoHPGeDetectorType.cc.

References test::a, FatalException, G4Exception(), and python.hepunit::keV.

{
  std::ostringstream ost;
  
  ost << fileName<<".dat";
  
  G4String name = ost.str();
  
  char* path = getenv("XRAYDATA");
  
  G4String pathString(path);
  G4String dirFile = pathString + "/" + name;
  std::ifstream file(dirFile);
  std::filebuf* lsdp = file.rdbuf();
  
  if (! (lsdp->is_open()) )
    {
      G4ExceptionDescription execp;
      execp << "XrayFluoHPGeDetectorType - data file: " + dirFile + " not found";
      G4Exception("XrayFluoHPGeDetectorType::LoadResponseData()","example-xray_fluorescence02",
      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();    
}

G4double XrayFluoHPGeDetectorType::ResponseFunction ( G4double  energy )

virtual

Implements XrayFluoVDetectorType.

Definition at line 98 of file XrayFluoHPGeDetectorType.cc.

References energy(), XrayFluoDataSet::FindValue(), G4UniformRand, GetInfData(), GetSupData(), and python.hepunit::keV.

{
 

  G4double eMin = 1* keV;
  G4double eMax = 10*keV; 
  G4double value = 0.;
  G4double efficiency = 1.;
  
  const XrayFluoDataSet* dataSet = efficiencySet;
  G4int id = 0;
  
  G4double random = G4UniformRand();
 
  if (energy>=eMin && energy <=eMax)
    {
      G4double infEnergy = (G4int)(energy/keV)* keV;
      
      G4double supEnergy = ((G4int)(energy/keV) + 1)*keV;
            
      G4double infData = GetInfData(energy, random, 0);// 0 is not used
      
      G4double supData = GetSupData(energy,random, 0); // 0 is not used
      
      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)
    { 
      G4double infEnergy = eMin;
      G4double supEnergy = eMin/keV +1*keV;
 
      G4double infData = GetInfData(eMin, random, 0);
      G4double supData = GetSupData(eMin,random, 0);
      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)
    { 
      G4double infEnergy = eMax/keV - 1. *keV;
      G4double supEnergy = eMax;
 
      G4double infData = GetInfData(eMax, random, 0);
      G4double supData = GetSupData(eMax,random, 0);
      value = (std::log10(infData)*std::log10(supEnergy/eMax) +
           std::log10(supData)*std::log10(eMax/infEnergy)) / 
    std::log10(supEnergy/infEnergy);
      value = std::pow(10,value);
      value = value+energy- eMax;
    }
  G4double  RandomNum = G4UniformRand(); 
  
  efficiency = dataSet->FindValue(value,id);
  if ( RandomNum>efficiency )
    {
      value = 0.;
    }
 
  return value;
}

---

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

• XrayFluoHPGeDetectorType.hh
• XrayFluoHPGeDetectorType.cc