G4PWATotalXsecZ Class Reference

#include <G4PWATotalXsecTable.hh>

Collaboration diagram for G4PWATotalXsecZ:

Public Member Functions
G4int	GetLowestEnergyBinIndex () const
G4int	GetHighestEnergyBinIndex () const
G4double	GetLowestEnergy () const
G4double	GetHighestEnergy () const
G4double	GetLowestXsecValue (G4int j) const
G4double	GetHighestXsecValue (G4int j) const
G4int	GetPWATotalXsecEnergyBinIndex (G4double energy) const
G4double	GetInterpXsec (G4double energy, G4int elowindex, G4int j) const
G4double	GetInterpXsec (G4double energy, G4int j) const

Private Member Functions
	G4PWATotalXsecZ (G4int Z)
	~G4PWATotalXsecZ ()
G4PWATotalXsecZ &	operator= (const G4PWATotalXsecZ &right)
	G4PWATotalXsecZ (const G4PWATotalXsecZ &)
void	LoadPWATotalXsecZ (G4int Z)

Private Attributes
G4double	fPWAXsecs [fgNumTotalXsecBins *6]
G4double	fInterpParamA [fgNumTotalXsecBins *6]
G4double	fInterpParamB [fgNumTotalXsecBins *6]

Static Private Attributes
static const G4int	fgNumTotalXsecBins = 106
static const G4double	fgPWATotalXsecEnergyGrid [fgNumTotalXsecBins]

Friends
class	G4PWATotalXsecTable

Detailed Description

Definition at line 64 of file G4PWATotalXsecTable.hh.

Constructor & Destructor Documentation

G4PWATotalXsecZ() [1/2]

G4PWATotalXsecZ::G4PWATotalXsecZ ( G4int  Z )

private

Definition at line 94 of file G4PWATotalXsecTable.cc.

                                       {
  G4int nn = fgNumTotalXsecBins*6;
  for(G4int i=0; i<nn; ++i) {
    fPWAXsecs[i] = 0.0;
    fInterpParamA[i] = 0.0;
    fInterpParamB[i] = 0.0;
  }
  LoadPWATotalXsecZ(Z); 
}

Here is the call graph for this function:

Here is the caller graph for this function:

~G4PWATotalXsecZ()

G4PWATotalXsecZ::~G4PWATotalXsecZ ( )

inlineprivate

Definition at line 104 of file G4PWATotalXsecTable.hh.

{};

Here is the call graph for this function:

G4PWATotalXsecZ() [2/2]

G4PWATotalXsecZ::G4PWATotalXsecZ ( const G4PWATotalXsecZ &  )

private

Member Function Documentation

GetHighestEnergy()

G4double G4PWATotalXsecZ::GetHighestEnergy ( ) const

inline

Definition at line 73 of file G4PWATotalXsecTable.hh.

{return fgPWATotalXsecEnergyGrid[fgNumTotalXsecBins-1];}

Here is the caller graph for this function:

GetHighestEnergyBinIndex()

G4int G4PWATotalXsecZ::GetHighestEnergyBinIndex ( ) const

inline

Definition at line 71 of file G4PWATotalXsecTable.hh.

{return fgNumTotalXsecBins-1;}

GetHighestXsecValue()

G4double G4PWATotalXsecZ::GetHighestXsecValue ( G4int  j ) const

inline

Definition at line 77 of file G4PWATotalXsecTable.hh.

{return fPWAXsecs[(j+1)*fgNumTotalXsecBins-1];}

Here is the call graph for this function:

Here is the caller graph for this function:

GetInterpXsec() [1/2]

G4double G4PWATotalXsecZ::GetInterpXsec	(	G4double	energy,
		G4int	elowindex,
		G4int	j
	)		const

Definition at line 166 of file G4PWATotalXsecTable.cc.

                                                                                       { 
   // protection : out of energy grid range 
   if(energy < GetLowestEnergy())
      return GetLowestXsecValue(j);
   if(energy >= GetHighestEnergy())
      return GetHighestXsecValue(j);     

   // normal case log-log linear intrp.
  G4int k = j*fgNumTotalXsecBins+elowindex; 
  return G4Exp(G4Log(energy)*fInterpParamA[k])*fInterpParamB[k];
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetInterpXsec() [2/2]

G4double G4PWATotalXsecZ::GetInterpXsec	(	G4double	energy,
		G4int	j
	)		const

Definition at line 179 of file G4PWATotalXsecTable.cc.

                                                                      {
   // protection : out of energy grid range 
   if(energy < GetLowestEnergy())
      return GetLowestXsecValue(j);
   if(energy >= GetHighestEnergy())
      return GetHighestXsecValue(j);     

   // normal case log-log linear intrp.
  G4int elowindex = GetPWATotalXsecEnergyBinIndex(energy); 
  G4int k = j*fgNumTotalXsecBins+elowindex; 
  return G4Exp(G4Log(energy)*fInterpParamA[k])*fInterpParamB[k];
}

Here is the call graph for this function:

GetLowestEnergy()

G4double G4PWATotalXsecZ::GetLowestEnergy ( ) const

inline

Definition at line 72 of file G4PWATotalXsecTable.hh.

{return fgPWATotalXsecEnergyGrid[0];}

Here is the caller graph for this function:

GetLowestEnergyBinIndex()

G4int G4PWATotalXsecZ::GetLowestEnergyBinIndex ( ) const

inline

Definition at line 70 of file G4PWATotalXsecTable.hh.

{return 0;}

GetLowestXsecValue()

G4double G4PWATotalXsecZ::GetLowestXsecValue ( G4int  j ) const

inline

Definition at line 76 of file G4PWATotalXsecTable.hh.

{return fPWAXsecs[j*fgNumTotalXsecBins];}

Here is the caller graph for this function:

GetPWATotalXsecEnergyBinIndex()

G4int G4PWATotalXsecZ::GetPWATotalXsecEnergyBinIndex

(

G4double

energy

)

const

Definition at line 154 of file G4PWATotalXsecTable.cc.

                                                                          {
  // log(fgPWATotalXsecEnergyGrid[0]);
  const G4double lne0 = -9.21034037197618e+00; 
  // 1./log(fgPWATotalXsecEnergyGrid[i+1]/fgPWATotalXsecEnergyGrid[i]); 
  const G4double invlnde =  6.51441722854880e+00;   

  return (G4int)((G4Log(energy)-lne0)*invlnde);
}

Here is the call graph for this function:

Here is the caller graph for this function:

LoadPWATotalXsecZ()

void G4PWATotalXsecZ::LoadPWATotalXsecZ ( G4int  Z )

private

Definition at line 105 of file G4PWATotalXsecTable.cc.

                                              {
 G4double dum;
 char fname[512];
 char* path = getenv("G4LEDATA");
 if (!path) {
    G4Exception("G4PWATotalXsecZ::LoadPWATotalXsecZ()","em0006",
        FatalException,
        "Environment variable G4LEDATA not defined");
    return;
 }

 std::string pathString(path);
 sprintf(fname,"%s/msc_GS/xsecs/xsecs_%d",path,Z);
 std::ifstream infile(fname,std::ios::in);
 if(!infile.is_open()){
    char msgc[512];
    sprintf(msgc,"  Total PWA xsection %s not found.",fname);
    G4Exception("G4PWATotalXsecZ::LoadPWATotalXsecZ()","em0006",
        FatalException,
        msgc);
    return;
 }
 G4double dummy;

 for(G4int i=0; i<fgNumTotalXsecBins; ++i)
   for(G4int j=0; j<7; ++j)
     if(j==0) infile >> dum;
     else {
       // load pwa xsection that are stored in cm2 units in file and change to 
       // Geant4 internal length2 units 
       infile >> dummy;
       fPWAXsecs[(j-1)*fgNumTotalXsecBins+i] = dummy*CLHEP::cm2;        
     }
 infile.close();

 // compute log-log linear intrp. parameters
 for(G4int i=0; i<fgNumTotalXsecBins-1; ++i)
    for(G4int k=0; k<6; ++k) {
       G4int j = k*fgNumTotalXsecBins+i; 
       G4double val2 = fPWAXsecs[j+1];
       G4double val1 = fPWAXsecs[j];

       fInterpParamA[j] =  G4Log(val2/val1)/G4Log(fgPWATotalXsecEnergyGrid[i+1]/fgPWATotalXsecEnergyGrid[i]);
       fInterpParamB[j] =  G4Exp(G4Log(val1) - fInterpParamA[j]*G4Log(fgPWATotalXsecEnergyGrid[i]));
 }
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

G4PWATotalXsecZ& G4PWATotalXsecZ::operator= ( const G4PWATotalXsecZ &  right )

private

Here is the caller graph for this function:

Friends And Related Function Documentation

G4PWATotalXsecTable

friend class G4PWATotalXsecTable

friend

Definition at line 66 of file G4PWATotalXsecTable.hh.

Member Data Documentation

fgNumTotalXsecBins

const G4int G4PWATotalXsecZ::fgNumTotalXsecBins = 106

staticprivate

Definition at line 114 of file G4PWATotalXsecTable.hh.

fgPWATotalXsecEnergyGrid

const G4double G4PWATotalXsecZ::fgPWATotalXsecEnergyGrid

staticprivate

Initial value:

={
    
    1.00000000e-04, 1.16591440e-04, 1.35935639e-04, 1.58489319e-04, 1.84784980e-04, 2.15443469e-04, 2.51188643e-04, 2.92864456e-04, 
    3.41454887e-04, 3.98107171e-04, 4.64158883e-04, 5.41169527e-04, 6.30957344e-04, 7.35642254e-04, 8.57695899e-04, 1.00000000e-03, 
    1.16591440e-03, 1.35935639e-03, 1.58489319e-03, 1.84784980e-03, 2.15443469e-03, 2.51188643e-03, 2.92864456e-03, 3.41454887e-03, 
    3.98107171e-03, 4.64158883e-03, 5.41169527e-03, 6.30957344e-03, 7.35642254e-03, 8.57695899e-03, 1.00000000e-02, 1.16591440e-02, 
    1.35935639e-02, 1.58489319e-02, 1.84784980e-02, 2.15443469e-02, 2.51188643e-02, 2.92864456e-02, 3.41454887e-02, 3.98107171e-02, 
    4.64158883e-02, 5.41169527e-02, 6.30957344e-02, 7.35642254e-02, 8.57695899e-02, 1.00000000e-01, 1.16591440e-01, 1.35935639e-01, 
    1.58489319e-01, 1.84784980e-01, 2.15443469e-01, 2.51188643e-01, 2.92864456e-01, 3.41454887e-01, 3.98107171e-01, 4.64158883e-01, 
    5.41169527e-01, 6.30957344e-01, 7.35642254e-01, 8.57695899e-01, 1.00000000e+00, 1.16591440e+00, 1.35935639e+00, 1.58489319e+00, 
    1.84784980e+00, 2.15443469e+00, 2.51188643e+00, 2.92864456e+00, 3.41454887e+00, 3.98107171e+00, 4.64158883e+00, 5.41169527e+00, 
    6.30957344e+00, 7.35642254e+00, 8.57695899e+00, 1.00000000e+01, 1.16591440e+01, 1.35935639e+01, 1.58489319e+01, 1.84784980e+01, 
    2.15443469e+01, 2.51188643e+01, 2.92864456e+01, 3.41454887e+01, 3.98107171e+01, 4.64158883e+01, 5.41169527e+01, 6.30957344e+01, 
    7.35642254e+01, 8.57695899e+01, 1.00000000e+02, 1.16591440e+02, 1.35935639e+02, 1.58489319e+02, 1.84784980e+02, 2.15443469e+02, 
    2.51188643e+02, 2.92864456e+02, 3.41454887e+02, 3.98107171e+02, 4.64158883e+02, 5.41169527e+02, 6.30957344e+02, 7.35642254e+02, 
    8.57695899e+02, 1.00000000e+03
 }

Definition at line 118 of file G4PWATotalXsecTable.hh.

fInterpParamA

G4double G4PWATotalXsecZ::fInterpParamA[fgNumTotalXsecBins *6]

private

Definition at line 124 of file G4PWATotalXsecTable.hh.

fInterpParamB

G4double G4PWATotalXsecZ::fInterpParamB[fgNumTotalXsecBins *6]

private

Definition at line 125 of file G4PWATotalXsecTable.hh.

fPWAXsecs

G4double G4PWATotalXsecZ::fPWAXsecs[fgNumTotalXsecBins *6]

private

Definition at line 122 of file G4PWATotalXsecTable.hh.

---

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

• G4PWATotalXsecTable.hh
• G4PWATotalXsecTable.cc