G4ComponentSAIDTotalXS Class Reference

#include <G4ComponentSAIDTotalXS.hh>

Inheritance diagram for G4ComponentSAIDTotalXS:

Collaboration diagram for G4ComponentSAIDTotalXS:

Public Member Functions
	G4ComponentSAIDTotalXS ()
virtual	~G4ComponentSAIDTotalXS ()
virtual G4double	GetTotalElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4double)
virtual G4double	GetTotalIsotopeCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)
virtual G4double	GetInelasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4double)
virtual G4double	GetInelasticIsotopeCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)
virtual G4double	GetElasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4double)
virtual G4double	GetElasticIsotopeCrossSection (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)
G4double	GetChargeExchangeCrossSection (const G4ParticleDefinition *prim, const G4ParticleDefinition *sec, G4double kinEnergy, G4int, G4int)
virtual void	Description () const
Public Member Functions inherited from G4VComponentCrossSection
	G4VComponentCrossSection (const G4String &nam="")
virtual	~G4VComponentCrossSection ()
G4double	GetTotalElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, const G4Element *)
G4double	GetInelasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, const G4Element *)
G4double	GetElasticElementCrossSection (const G4ParticleDefinition *, G4double kinEnergy, const G4Element *)
virtual G4double	ComputeQuasiElasticRatio (const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
const G4String &	GetName () const

Private Member Functions
G4SAIDCrossSectionType	GetType (const G4ParticleDefinition *prim, const G4ParticleDefinition *sec, G4int Z, G4int N)
void	Initialise (G4SAIDCrossSectionType tp)
void	ReadData (G4int index, G4PhysicsVector *, const G4String &, const G4String &)
void	PrintWarning (const G4ParticleDefinition *prim, const G4ParticleDefinition *sec, G4int, G4int, const G4String &, const G4String &)
G4ComponentSAIDTotalXS &	operator= (const G4ComponentSAIDTotalXS &right)
	G4ComponentSAIDTotalXS (const G4ComponentSAIDTotalXS &)

Private Attributes
G4PhysicsVector *	elastdata [numberOfSaidXS]
G4PhysicsVector *	inelastdata [numberOfSaidXS]

Static Private Attributes
static const G4String	fnames [numberOfSaidXS]

Detailed Description

Definition at line 74 of file G4ComponentSAIDTotalXS.hh.

Constructor & Destructor Documentation

G4ComponentSAIDTotalXS() [1/2]

G4ComponentSAIDTotalXS::G4ComponentSAIDTotalXS

(

)

Definition at line 50 of file G4ComponentSAIDTotalXS.cc.

  : G4VComponentCrossSection("xsSAID")
{
  for(G4int i=0; i<numberOfSaidXS; ++i) {
    elastdata[i] = 0;
    inelastdata[i] = 0;
  }
}

~G4ComponentSAIDTotalXS()

G4ComponentSAIDTotalXS::~G4ComponentSAIDTotalXS ( )

virtual

Definition at line 59 of file G4ComponentSAIDTotalXS.cc.

{
  for(G4int i=0; i<numberOfSaidXS; ++i) {
    if(elastdata[i]) {
      delete elastdata[i];
      elastdata[i] = 0;
    }
    if(inelastdata[i]) {
      delete inelastdata[i];
      inelastdata[i] = 0;
    }
  }
}

G4ComponentSAIDTotalXS() [2/2]

G4ComponentSAIDTotalXS::G4ComponentSAIDTotalXS ( const G4ComponentSAIDTotalXS &  )

private

Member Function Documentation

Description()

void G4ComponentSAIDTotalXS::Description ( ) const

virtual

Reimplemented from G4VComponentCrossSection.

Definition at line 168 of file G4ComponentSAIDTotalXS.cc.

{
}

GetChargeExchangeCrossSection()

G4double G4ComponentSAIDTotalXS::GetChargeExchangeCrossSection	(	const G4ParticleDefinition *	prim,
		const G4ParticleDefinition *	sec,
		G4double	kinEnergy,
		G4int	Z,
		G4int	N
	)

Definition at line 150 of file G4ComponentSAIDTotalXS.cc.

{
  G4double cross = 0.0;
  G4SAIDCrossSectionType tp = GetType(prim,sec,Z,N);
  if(saidUnknown != tp) {
    G4int idx = G4int(tp);
    if(!inelastdata[idx]) { Initialise(tp); }
    if(inelastdata[idx]) { 
      cross = (inelastdata[idx])->Value(kinEnergy); 
    }
  }
  return cross;
}

Here is the call graph for this function:

GetElasticElementCrossSection()

G4double G4ComponentSAIDTotalXS::GetElasticElementCrossSection ( const G4ParticleDefinition *  part, G4double  kinEnergy, G4int  Z, G4double  N  )

virtual

Implements G4VComponentCrossSection.

Definition at line 122 of file G4ComponentSAIDTotalXS.cc.

{
  PrintWarning(part,0,Z,G4lrint(N),
           "G4ComponentSAIDTotalXS::GetTotalElementCrossSection",
           "Method is not implemented");
  return 0.0; 
}

Here is the call graph for this function:

GetElasticIsotopeCrossSection()

G4double G4ComponentSAIDTotalXS::GetElasticIsotopeCrossSection ( const G4ParticleDefinition *  part, G4double  kinEnergy, G4int  Z, G4int  N  )

virtual

Implements G4VComponentCrossSection.

Definition at line 133 of file G4ComponentSAIDTotalXS.cc.

{
  G4double cross = 0.0;
  G4SAIDCrossSectionType tp = GetType(part,0,Z,N);
  if(saidUnknown != tp) {
    G4int idx = G4int(tp);
    if(!elastdata[idx]) { Initialise(tp); }
    if(elastdata[idx]) { 
      cross = (elastdata[idx])->Value(kinEnergy); 
    }
  }
  return cross;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetInelasticElementCrossSection()

G4double G4ComponentSAIDTotalXS::GetInelasticElementCrossSection ( const G4ParticleDefinition *  part, G4double  kinEnergy, G4int  Z, G4double  N  )

virtual

Implements G4VComponentCrossSection.

Definition at line 94 of file G4ComponentSAIDTotalXS.cc.

{
  PrintWarning(part,0,Z,G4lrint(N),
           "G4ComponentSAIDTotalXS::GetTotalElementCrossSection",
           "Method is not implemented");
  return 0.0; 
}

Here is the call graph for this function:

GetInelasticIsotopeCrossSection()

G4double G4ComponentSAIDTotalXS::GetInelasticIsotopeCrossSection ( const G4ParticleDefinition *  part, G4double  kinEnergy, G4int  Z, G4int  N  )

virtual

Implements G4VComponentCrossSection.

Definition at line 105 of file G4ComponentSAIDTotalXS.cc.

{
  G4double cross = 0.0;
  G4SAIDCrossSectionType tp = GetType(part,0,Z,N);
  if(saidUnknown != tp) {
    G4int idx = G4int(tp);
    if(!inelastdata[idx]) { Initialise(tp); }
    if(inelastdata[idx]) { 
      cross = (inelastdata[idx])->Value(kinEnergy);
    }
  }
  return cross;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetTotalElementCrossSection()

G4double G4ComponentSAIDTotalXS::GetTotalElementCrossSection ( const G4ParticleDefinition *  part, G4double  kinEnergy, G4int  Z, G4double  N  )

virtual

Implements G4VComponentCrossSection.

Definition at line 74 of file G4ComponentSAIDTotalXS.cc.

{
  PrintWarning(part,0,Z,G4lrint(N),
           "G4ComponentSAIDTotalXS::GetTotalElementCrossSection",
           "Method is not implemented");
  return 0.0; 
}

Here is the call graph for this function:

GetTotalIsotopeCrossSection()

G4double G4ComponentSAIDTotalXS::GetTotalIsotopeCrossSection ( const G4ParticleDefinition *  part, G4double  kinEnergy, G4int  Z, G4int  N  )

virtual

Implements G4VComponentCrossSection.

Definition at line 85 of file G4ComponentSAIDTotalXS.cc.

{
  return GetInelasticIsotopeCrossSection(part,kinEnergy,Z,N)
    + GetElasticIsotopeCrossSection(part,kinEnergy,Z,N);
}

Here is the call graph for this function:

GetType()

G4SAIDCrossSectionType G4ComponentSAIDTotalXS::GetType ( const G4ParticleDefinition *  prim, const G4ParticleDefinition *  sec, G4int  Z, G4int  N  )

private

Definition at line 173 of file G4ComponentSAIDTotalXS.cc.

{
  G4SAIDCrossSectionType type = saidUnknown;
  if(1 == N && prim) {
    G4int code = prim->GetPDGEncoding();

    // only gamma + N x-sections available
    if(0 == Z && sec && 22 == code) {
      G4int code1 = sec->GetPDGEncoding();
      if(-211 == code1)     { type = saidGN_PINP; }
      else if(111 == code1) { type = saidGN_PI0N; }

      // x-sections off proton
    } else if(1 == Z) {
      if(sec) {
    G4int code1 = sec->GetPDGEncoding();
        if(-211 == code) {
          if(111 == code1)      { type = saidPINP_PI0N; }
          else if(221 == code1) { type = saidPINP_ETAN; }

    } else if(22 == code) {
          if(111 == code1)      { type = saidGP_PI0P; }
          else if(211 == code1) { type = saidGP_PIPN; }
          else if(221 == code1) { type = saidGP_ETAP; }
          else if(331 == code1) { type = saidGP_ETAPP; }
    }
      } else {
        if(2212 == code)        { type = saidPP; } 
        else if(2112 == code)   { type = saidNP; } 
        else if(211 == code)    { type = saidPIPP; } 
        else if(-211 == code)   { type = saidPINP; } 
      }
    }
  }
  //G4cout << "G4ComponentSAIDTotalXS::Type= " << type << G4endl;
  return type;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Initialise()

void G4ComponentSAIDTotalXS::Initialise ( G4SAIDCrossSectionType  tp )

private

Definition at line 213 of file G4ComponentSAIDTotalXS.cc.

{
  G4int idx = G4int(tp);
  // check environment variable 
  // Build the complete string identifying the file with the data set
  char* path = getenv("G4SAIDXSDATA");
  if (!path){
    G4Exception("G4ComponentSAIDTotalXS::Initialise(..)","had013",
        FatalException,
        "Environment variable G4SAIDXSDATA is not defined");
    return;
  }
  if(idx <= 4) {
    elastdata[idx] = new G4LPhysicsFreeVector();
    inelastdata[idx] = new G4LPhysicsFreeVector();
    ReadData(idx,elastdata[idx],path,"_el.dat");
    ReadData(idx,inelastdata[idx],path,"_in.dat");
  } else {
    inelastdata[idx] = new G4LPhysicsFreeVector();
    ReadData(idx,inelastdata[idx],path,".dat");
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

operator=()

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

private

PrintWarning()

void G4ComponentSAIDTotalXS::PrintWarning ( const G4ParticleDefinition *  prim, const G4ParticleDefinition *  sec, G4int  Z, G4int  N, const G4String &  ss1, const G4String &  ss2  )

private

Definition at line 263 of file G4ComponentSAIDTotalXS.cc.

{
  G4cout << ss1 << ": " << ss2 << G4endl;
  G4cout << "For Z= " << Z << " N= " << N << " of ";
  if(prim) { G4cout << prim->GetParticleName() << " "; }
  if(sec) { G4cout << " x-section to " << sec->GetParticleName(); }
  G4cout << G4endl;
}

Here is the call graph for this function:

Here is the caller graph for this function:

ReadData()

void G4ComponentSAIDTotalXS::ReadData ( G4int  index, G4PhysicsVector *  v, const G4String &  ss1, const G4String &  ss2  )

private

Definition at line 236 of file G4ComponentSAIDTotalXS.cc.

{
  std::ostringstream ost;
  ost << ss1 << "/" << fnames[index] << ss2;
  std::ifstream filein(ost.str().c_str());
  if (!(filein)) {
    G4ExceptionDescription ed;
    ed << "Data file <" << ost.str().c_str()
       << "> is not opened!";
    G4Exception("G4ComponentSAIDTotalXS::ReadData(..)","had014",
                FatalException, ed, "Check G4SAIDXSDATA");
  } else {
    if(GetVerboseLevel() > 1) {
      G4cout << "File " << ost.str() 
             << " is opened by G4ComponentSAIDTotalXS" << G4endl;
    }
    // retrieve data from DB
    v->Retrieve(filein, true);
    v->ScaleVector(CLHEP::MeV,CLHEP::millibarn);
    v->SetSpline(true);
  } 
}

Here is the call graph for this function:

Here is the caller graph for this function:

Member Data Documentation

elastdata

G4PhysicsVector* G4ComponentSAIDTotalXS::elastdata[numberOfSaidXS]

private

Definition at line 140 of file G4ComponentSAIDTotalXS.hh.

fnames

const G4String G4ComponentSAIDTotalXS::fnames

staticprivate

Initial value:

= {
  "","pp","np","pip","pim",
  "pin","pie",
  "gp_pi0p","gp_pi+n","gn_pi-p","gn_pi0n","gp_etap","gp_etapp"
}

Definition at line 139 of file G4ComponentSAIDTotalXS.hh.

inelastdata

G4PhysicsVector* G4ComponentSAIDTotalXS::inelastdata[numberOfSaidXS]

private

Definition at line 141 of file G4ComponentSAIDTotalXS.hh.

---

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

• G4ComponentSAIDTotalXS.hh
• G4ComponentSAIDTotalXS.cc