G4ChipsHyperonElasticXS Class Reference

#include <G4ChipsHyperonElasticXS.hh>

Inheritance diagram for G4ChipsHyperonElasticXS:

Collaboration diagram for G4ChipsHyperonElasticXS:

Public Member Functions
	G4ChipsHyperonElasticXS ()
	~G4ChipsHyperonElasticXS ()
virtual void	CrossSectionDescription (std::ostream &) const
virtual G4bool	IsIsoApplicable (const G4DynamicParticle *Pt, G4int Z, G4int A, const G4Element *elm, const G4Material *mat)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int tgZ, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
virtual G4double	GetChipsCrossSection (G4double momentum, G4int Z, G4int N, G4int pdg)
G4double	GetExchangeT (G4int tZ, G4int tN, G4int pPDG)
Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")
virtual	~G4VCrossSectionDataSet ()
virtual G4bool	IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0)
virtual G4double	GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=0)
virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual G4int	GetVerboseLevel () const
virtual void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
const G4String &	GetName () const

Static Public Member Functions
static const char *	Default_Name ()

Private Member Functions
G4double	CalculateCrossSection (G4bool CS, G4int F, G4int I, G4int pPDG, G4int Z, G4int N, G4double pP)
G4double	GetSlope (G4int tZ, G4int tN, G4int pPDG)
G4double	GetHMaxT ()
G4double	GetPTables (G4double lpP, G4double lPm, G4int PDG, G4int tZ, G4int tN)
G4double	GetTabValues (G4double lp, G4int pPDG, G4int tgZ, G4int tgN)
G4double	GetQ2max (G4int pPDG, G4int tgZ, G4int tgN, G4double pP)

Private Attributes
const G4int	nPoints
const G4int	nLast
G4double	lPMin
G4double	lPMax
G4double	dlnP
G4bool	onlyCS
G4double	lastSIG
G4double	lastLP
G4double	lastTM
G4int	lastN
G4int	lastZ
G4double	lastP
G4double	lastTH
G4double	lastCS
G4int	lastI
G4double	theSS
G4double	theS1
G4double	theB1
G4double	theS2
G4double	theB2
G4double	theS3
G4double	theB3
G4double	theS4
G4double	theB4
G4int	lastTZ
G4int	lastTN
G4double	lastPIN
G4double *	lastCST
G4double *	lastPAR
G4double *	lastSST
G4double *	lastS1T
G4double *	lastB1T
G4double *	lastS2T
G4double *	lastB2T
G4double *	lastS3T
G4double *	lastB3T
G4double *	lastS4T
G4double *	lastB4T
std::vector< G4double * >	PAR
std::vector< G4double * >	CST
std::vector< G4double * >	SST
std::vector< G4double * >	S1T
std::vector< G4double * >	B1T
std::vector< G4double * >	S2T
std::vector< G4double * >	B2T
std::vector< G4double * >	S3T
std::vector< G4double * >	B3T
std::vector< G4double * >	S4T
std::vector< G4double * >	B4T
std::vector< G4int >	colN
std::vector< G4int >	colZ
std::vector< G4double >	colP
std::vector< G4double >	colTH
std::vector< G4double >	colCS
std::vector< G4double >	PIN

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 46 of file G4ChipsHyperonElasticXS.hh.

Constructor & Destructor Documentation

G4ChipsHyperonElasticXS()

G4ChipsHyperonElasticXS::G4ChipsHyperonElasticXS

(

)

Definition at line 64 of file G4ChipsHyperonElasticXS.cc.

                                                :G4VCrossSectionDataSet(Default_Name()), nPoints(128), nLast(nPoints-1)
{
  lPMin=-8.;  //Min tabulatedLogarithmMomentum(D)
  lPMax= 8.;  //Max tabulatedLogarithmMomentum(D)
  dlnP=(lPMax-lPMin)/nLast;// LogStep inTable (D)
  onlyCS=true;//Flag toCalculOnlyCS(not Si/Bi)(L)
  lastSIG=0.; //Last calculated cross section (L)
  lastLP=-10.;//LastLog(mom_of IncidentHadron)(L)
  lastTM=0.;  //Last t_maximum                (L)
  theSS=0.;   //TheLastSqSlope of 1st difr.Max(L)
  theS1=0.;   //TheLastMantissa of 1st difrMax(L)
  theB1=0.;   //TheLastSlope of 1st difructMax(L)
  theS2=0.;   //TheLastMantissa of 2nd difrMax(L)
  theB2=0.;   //TheLastSlope of 2nd difructMax(L)
  theS3=0.;   //TheLastMantissa of 3d difr.Max(L)
  theB3=0.;   //TheLastSlope of 3d difruct.Max(L)
  theS4=0.;   //TheLastMantissa of 4th difrMax(L)
  theB4=0.;   //TheLastSlope of 4th difructMax(L)
  lastTZ=0;   // Last atomic number of the target
  lastTN=0;   // Last # of neutrons in the target
  lastPIN=0.; // Last initialized max momentum
  lastCST=0;  // Elastic cross-section table
  lastPAR=0;  // ParametersForFunctionCalculation
  lastSST=0;  // E-dep ofSqardSlope of 1st difMax
  lastS1T=0;  // E-dep of mantissa of 1st dif.Max
  lastB1T=0;  // E-dep of the slope of 1st difMax
  lastS2T=0;  // E-dep of mantissa of 2nd difrMax
  lastB2T=0;  // E-dep of the slope of 2nd difMax
  lastS3T=0;  // E-dep of mantissa of 3d difr.Max
  lastB3T=0;  // E-dep of the slope of 3d difrMax
  lastS4T=0;  // E-dep of mantissa of 4th difrMax
  lastB4T=0;  // E-dep of the slope of 4th difMax
  lastN=0;    // The last N of calculated nucleus
  lastZ=0;    // The last Z of calculated nucleus
  lastP=0.;   // LastUsed inCrossSection Momentum
  lastTH=0.;  // Last threshold momentum
  lastCS=0.;  // Last value of the Cross Section
  lastI=0;    // The last position in the DAMDB
}

~G4ChipsHyperonElasticXS()

G4ChipsHyperonElasticXS::~G4ChipsHyperonElasticXS

(

)

Definition at line 104 of file G4ChipsHyperonElasticXS.cc.

{
  std::vector<G4double*>::iterator pos;
  for (pos=CST.begin(); pos<CST.end(); pos++)
  { delete [] *pos; }
  CST.clear();
  for (pos=PAR.begin(); pos<PAR.end(); pos++)
  { delete [] *pos; }
  PAR.clear();
  for (pos=SST.begin(); pos<SST.end(); pos++)
  { delete [] *pos; }
  SST.clear();
  for (pos=S1T.begin(); pos<S1T.end(); pos++)
  { delete [] *pos; }
  S1T.clear();
  for (pos=B1T.begin(); pos<B1T.end(); pos++)
  { delete [] *pos; }
  B1T.clear();
  for (pos=S2T.begin(); pos<S2T.end(); pos++)
  { delete [] *pos; }
  S2T.clear();
  for (pos=B2T.begin(); pos<B2T.end(); pos++)
  { delete [] *pos; }
  B2T.clear();
  for (pos=S3T.begin(); pos<S3T.end(); pos++)
  { delete [] *pos; }
  S3T.clear();
  for (pos=B3T.begin(); pos<B3T.end(); pos++)
  { delete [] *pos; }
  B3T.clear();
  for (pos=S4T.begin(); pos<S4T.end(); pos++)
  { delete [] *pos; }
  S4T.clear();
  for (pos=B4T.begin(); pos<B4T.end(); pos++)
  { delete [] *pos; }
  B4T.clear();
}

Member Function Documentation

CalculateCrossSection()

G4double G4ChipsHyperonElasticXS::CalculateCrossSection ( G4bool  CS, G4int  F, G4int  I, G4int  pPDG, G4int  Z, G4int  N, G4double  pP  )

private

Definition at line 242 of file G4ChipsHyperonElasticXS.cc.

{
  G4double pMom=pIU/GeV;                // All calculations are in GeV
  onlyCS=CS;                            // Flag to calculate only CS (not Si/Bi)
  lastLP=G4Log(pMom);                // Make a logarithm of the momentum for calculation
  if(F)                                 // This isotope was found in AMDB =>RETRIEVE/UPDATE
  {
    if(F<0)                             // the AMDB must be loded
    {
      lastPIN = PIN[I];                 // Max log(P) initialised for this table set
      lastPAR = PAR[I];                 // Pointer to the parameter set

      lastCST = CST[I];                 // Pointer to the total sross-section table
      lastSST = SST[I];                 // Pointer to the first squared slope
      lastS1T = S1T[I];                 // Pointer to the first mantissa
      lastB1T = B1T[I];                 // Pointer to the first slope
      lastS2T = S2T[I];                 // Pointer to the second mantissa
      lastB2T = B2T[I];                 // Pointer to the second slope
      lastS3T = S3T[I];                 // Pointer to the third mantissa
      lastB3T = B3T[I];                 // Pointer to the rhird slope
      lastS4T = S4T[I];                 // Pointer to the 4-th mantissa
      lastB4T = B4T[I];                 // Pointer to the 4-th slope
    }
    if(lastLP>lastPIN && lastLP<lPMax)
    {
      lastPIN=GetPTables(lastLP,lastPIN,PDG,tgZ,tgN);// Can update upper logP-Limit in tabs
      PIN[I]=lastPIN;                   // Remember the new P-Limit of the tables
    }
  }
  else                                  // This isotope wasn't initialized => CREATE
  {
    lastPAR = new G4double[nPoints];    // Allocate memory for parameters of CS function
    lastPAR[nLast]=0;                   // Initialization for VALGRIND
    lastCST = new G4double[nPoints];    // Allocate memory for Tabulated CS function    
    lastSST = new G4double[nPoints];    // Allocate memory for Tabulated first sqaredSlope 
    lastS1T = new G4double[nPoints];    // Allocate memory for Tabulated first mantissa 
    lastB1T = new G4double[nPoints];    // Allocate memory for Tabulated first slope    
    lastS2T = new G4double[nPoints];    // Allocate memory for Tabulated second mantissa
    lastB2T = new G4double[nPoints];    // Allocate memory for Tabulated second slope   
    lastS3T = new G4double[nPoints];    // Allocate memory for Tabulated third mantissa 
    lastB3T = new G4double[nPoints];    // Allocate memory for Tabulated third slope    
    lastS4T = new G4double[nPoints];    // Allocate memory for Tabulated 4-th mantissa 
    lastB4T = new G4double[nPoints];    // Allocate memory for Tabulated 4-th slope    
    lastPIN = GetPTables(lastLP,lPMin,PDG,tgZ,tgN); // Returns the new P-limit for tables
    PIN.push_back(lastPIN);             // Fill parameters of CS function to AMDB
    PAR.push_back(lastPAR);             // Fill parameters of CS function to AMDB
    CST.push_back(lastCST);             // Fill Tabulated CS function to AMDB    
    SST.push_back(lastSST);             // Fill Tabulated first sq.slope to AMDB 
    S1T.push_back(lastS1T);             // Fill Tabulated first mantissa to AMDB 
    B1T.push_back(lastB1T);             // Fill Tabulated first slope to AMDB    
    S2T.push_back(lastS2T);             // Fill Tabulated second mantissa to AMDB 
    B2T.push_back(lastB2T);             // Fill Tabulated second slope to AMDB    
    S3T.push_back(lastS3T);             // Fill Tabulated third mantissa to AMDB 
    B3T.push_back(lastB3T);             // Fill Tabulated third slope to AMDB    
    S4T.push_back(lastS4T);             // Fill Tabulated 4-th mantissa to AMDB 
    B4T.push_back(lastB4T);             // Fill Tabulated 4-th slope to AMDB    
  } // End of creation/update of the new set of parameters and tables
  // =-----------= NOW Update (if necessary) and Calculate the Cross Section =-----------=
  if(lastLP>lastPIN && lastLP<lPMax)
  {
    lastPIN = GetPTables(lastLP,lastPIN,PDG,tgZ,tgN);
  }
  if(!onlyCS) lastTM=GetQ2max(PDG, tgZ, tgN, pMom); // Calculate (-t)_max=Q2_max (GeV2)
  if(lastLP>lPMin && lastLP<=lastPIN)   // Linear fit is made using precalculated tables
  {
    if(lastLP==lastPIN)
    {
      G4double shift=(lastLP-lPMin)/dlnP+.000001; // Log distance from lPMin
      G4int    blast=static_cast<int>(shift); // this is a bin number of the lower edge (0)
      if(blast<0 || blast>=nLast)G4cout<<"G4QHyperElCS::CCS:b="<<blast<<","<<nLast<<G4endl;
      lastSIG = lastCST[blast];
      if(!onlyCS)                       // Skip the differential cross-section parameters
      {
        theSS  = lastSST[blast];
        theS1  = lastS1T[blast];
        theB1  = lastB1T[blast];
        theS2  = lastS2T[blast];
        theB2  = lastB2T[blast];
        theS3  = lastS3T[blast];
        theB3  = lastB3T[blast];
        theS4  = lastS4T[blast];
        theB4  = lastB4T[blast];
      }
    }
    else
    {
      G4double shift=(lastLP-lPMin)/dlnP;        // a shift from the beginning of the table
      G4int    blast=static_cast<int>(shift);    // the lower bin number
      if(blast<0)   blast=0;
      if(blast>=nLast) blast=nLast-1;            // low edge of the last bin
      shift-=blast;                              // step inside the unit bin
      G4int lastL=blast+1;                       // the upper bin number
      G4double SIGL=lastCST[blast];              // the basic value of the cross-section
      lastSIG= SIGL+shift*(lastCST[lastL]-SIGL); // calculated total elastic cross-section
      if(!onlyCS)                       // Skip the differential cross-section parameters
      {
        G4double SSTL=lastSST[blast];           // the low bin of the first squared slope
        theSS=SSTL+shift*(lastSST[lastL]-SSTL); // the basic value of the first sq.slope
        G4double S1TL=lastS1T[blast];           // the low bin of the first mantissa
        theS1=S1TL+shift*(lastS1T[lastL]-S1TL); // the basic value of the first mantissa
        G4double B1TL=lastB1T[blast];           // the low bin of the first slope
        theB1=B1TL+shift*(lastB1T[lastL]-B1TL); // the basic value of the first slope
        G4double S2TL=lastS2T[blast];           // the low bin of the second mantissa
        theS2=S2TL+shift*(lastS2T[lastL]-S2TL); // the basic value of the second mantissa
        G4double B2TL=lastB2T[blast];           // the low bin of the second slope
        theB2=B2TL+shift*(lastB2T[lastL]-B2TL); // the basic value of the second slope
        G4double S3TL=lastS3T[blast];           // the low bin of the third mantissa
        theS3=S3TL+shift*(lastS3T[lastL]-S3TL); // the basic value of the third mantissa
        G4double B3TL=lastB3T[blast];           // the low bin of the third slope
        theB3=B3TL+shift*(lastB3T[lastL]-B3TL); // the basic value of the third slope
        G4double S4TL=lastS4T[blast];           // the low bin of the 4-th mantissa
        theS4=S4TL+shift*(lastS4T[lastL]-S4TL); // the basic value of the 4-th mantissa
        G4double B4TL=lastB4T[blast];           // the low bin of the 4-th slope
        theB4=B4TL+shift*(lastB4T[lastL]-B4TL); // the basic value of the 4-th slope
      }
    }
  }
  else lastSIG=GetTabValues(lastLP, PDG, tgZ, tgN); // Direct calculation beyond the table
  if(lastSIG<0.) lastSIG = 0.;                   // @@ a Warning print can be added
  return lastSIG;
}

Here is the call graph for this function:

Here is the caller graph for this function:

CrossSectionDescription()

void G4ChipsHyperonElasticXS::CrossSectionDescription ( std::ostream &  outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 143 of file G4ChipsHyperonElasticXS.cc.

{
    outFile << "G4ChipsHyperonElasticXS provides the elastic cross\n"
            << "section for hyperon nucleus scattering as a function of incident\n"
            << "momentum. The cross section is calculated using M. Kossov's\n"
            << "CHIPS parameterization of cross section data.\n";
}

Here is the caller graph for this function:

Default_Name()

static const char* G4ChipsHyperonElasticXS::Default_Name ( )

inlinestatic

Definition at line 54 of file G4ChipsHyperonElasticXS.hh.

{return "ChipsHyperonElasticXS";}

Here is the call graph for this function:

Here is the caller graph for this function:

GetChipsCrossSection()

G4double G4ChipsHyperonElasticXS::GetChipsCrossSection ( G4double  momentum, G4int  Z, G4int  N, G4int  pdg  )

virtual

!The slave functions must provide cross-sections in millibarns (mb) !! (not in IU)

Definition at line 172 of file G4ChipsHyperonElasticXS.cc.

{

  G4bool fCS = false;
  G4double pEn=pMom;

  onlyCS=fCS;

  G4bool in=false;                   // By default the isotope must be found in the AMDB
  lastP   = 0.;                      // New momentum history (nothing to compare with)
  lastN   = tgN;                     // The last N of the calculated nucleus
  lastZ   = tgZ;                     // The last Z of the calculated nucleus
  lastI   = colN.size();             // Size of the Associative Memory DB in the heap
  if(lastI) for(G4int i=0; i<lastI; i++) // Loop over proj/tgZ/tgN lines of DB
  {                                  // The nucleus with projPDG is found in AMDB
    if(colN[i]==tgN && colZ[i]==tgZ) // Isotope is foind in AMDB
    {
      lastI=i;
      lastTH =colTH[i];              // Last THreshold (A-dependent)
      if(pEn<=lastTH)
      {
        return 0.;                   // Energy is below the Threshold value
      }
      lastP  =colP [i];              // Last Momentum  (A-dependent)
      lastCS =colCS[i];              // Last CrossSect (A-dependent)
      //  if(std::fabs(lastP/pMom-1.)<tolerance) //VI (do not use tolerance)
      if(lastP == pMom)              // Do not recalculate
      {
        CalculateCrossSection(fCS,-1,i,pPDG,lastZ,lastN,pMom); // Update param's only
        return lastCS*millibarn;     // Use theLastCS
      }
      in = true;                       // This is the case when the isotop is found in DB
      // Momentum pMom is in IU ! @@ Units
      lastCS=CalculateCrossSection(fCS,-1,i,pPDG,lastZ,lastN,pMom); // read & update
      if(lastCS<=0. && pEn>lastTH)    // Correct the threshold
      {
        lastTH=pEn;
      }
      break;                           // Go out of the LOOP with found lastI
    }
  } // End of attampt to find the nucleus in DB
  if(!in)                            // This nucleus has not been calculated previously
  {
    lastCS=CalculateCrossSection(fCS,0,lastI,pPDG,lastZ,lastN,pMom);//calculate&create
    if(lastCS<=0.)
    {
      lastTH = 0; //ThresholdEnergy(tgZ, tgN); // The Threshold Energy which is now the last
      if(pEn>lastTH)
      {
        lastTH=pEn;
      }
    }
    colN.push_back(tgN);
    colZ.push_back(tgZ);
    colP.push_back(pMom);
    colTH.push_back(lastTH);
    colCS.push_back(lastCS);
    return lastCS*millibarn;
  } // End of creation of the new set of parameters
  else
  {
    colP[lastI]=pMom;
    colCS[lastI]=lastCS;
  }
  return lastCS*millibarn;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetExchangeT()

G4double G4ChipsHyperonElasticXS::GetExchangeT	(	G4int	tZ,
		G4int	tN,
		G4int	pPDG
	)

Definition at line 610 of file G4ChipsHyperonElasticXS.cc.

{
  static const G4double GeVSQ=gigaelectronvolt*gigaelectronvolt;
  static const G4double third=1./3.;
  static const G4double fifth=1./5.;
  static const G4double sevth=1./7.;  
  //AR-04Jun2014  if(PDG==3222 || PDG<3000 || PDG>3334)G4cout<<"*Warning*G4QHyElCS::GET:PDG="<<PDG<<G4endl;
  if(PDG<3000 || PDG>3334)G4cout<<"*Warning*G4QHyElCS::GET:PDG="<<PDG<<G4endl;
  if(onlyCS)G4cout<<"*Warning*G4ChipsHyperonElasticXS::GetExchanT: onlyCS=1"<<G4endl;
  if(lastLP<-4.3) return lastTM*GeVSQ*G4UniformRand();// S-wave for p<14 MeV/c (kinE<.1MeV)
  G4double q2=0.;
  if(tgZ==1 && tgN==0)                // ===> p+p=p+p
  {
    G4double E1=lastTM*theB1;
    G4double R1=(1.-G4Exp(-E1));
    G4double E2=lastTM*theB2;
    G4double R2=(1.-G4Exp(-E2*E2*E2));
    G4double E3=lastTM*theB3;
    G4double R3=(1.-G4Exp(-E3));
    G4double I1=R1*theS1/theB1;
    G4double I2=R2*theS2;
    G4double I3=R3*theS3;
    G4double I12=I1+I2;
    G4double rand=(I12+I3)*G4UniformRand();
    if     (rand<I1 )
    {
      G4double ran=R1*G4UniformRand();
      if(ran>1.) ran=1.;
      q2=-G4Log(1.-ran)/theB1;
    }
    else if(rand<I12)
    {
      G4double ran=R2*G4UniformRand();
      if(ran>1.) ran=1.;
      q2=-G4Log(1.-ran);
      if(q2<0.) q2=0.;
      q2=G4Pow::GetInstance()->powA(q2,third)/theB2;
    }
    else
    {
      G4double ran=R3*G4UniformRand();
      if(ran>1.) ran=1.;
      q2=-G4Log(1.-ran)/theB3;
    }
  }
  else
  {
    G4double a=tgZ+tgN;
    G4double E1=lastTM*(theB1+lastTM*theSS);
    G4double R1=(1.-G4Exp(-E1));
    G4double tss=theSS+theSS; // for future solution of quadratic equation (imediate check)
    G4double tm2=lastTM*lastTM;
    G4double E2=lastTM*tm2*theB2;                   // power 3 for lowA, 5 for HighA (1st)
    if(a>6.5)E2*=tm2;                               // for heavy nuclei
    G4double R2=(1.-G4Exp(-E2));
    G4double E3=lastTM*theB3;
    if(a>6.5)E3*=tm2*tm2*tm2;                       // power 1 for lowA, 7 (2nd) for HighA
    G4double R3=(1.-G4Exp(-E3));
    G4double E4=lastTM*theB4;
    G4double R4=(1.-G4Exp(-E4));
    G4double I1=R1*theS1;
    G4double I2=R2*theS2;
    G4double I3=R3*theS3;
    G4double I4=R4*theS4;
    G4double I12=I1+I2;
    G4double I13=I12+I3;
    G4double rand=(I13+I4)*G4UniformRand();
    if(rand<I1)
    {
      G4double ran=R1*G4UniformRand();
      if(ran>1.) ran=1.;
      q2=-G4Log(1.-ran)/theB1;
      if(std::fabs(tss)>1.e-7) q2=(std::sqrt(theB1*(theB1+(tss+tss)*q2))-theB1)/tss;
    }
    else if(rand<I12)
    {
      G4double ran=R2*G4UniformRand();
      if(ran>1.) ran=1.;
      q2=-G4Log(1.-ran)/theB2;
      if(q2<0.) q2=0.;
      if(a<6.5) q2=G4Pow::GetInstance()->powA(q2,third);
      else      q2=G4Pow::GetInstance()->powA(q2,fifth);
    }
    else if(rand<I13)
    {
      G4double ran=R3*G4UniformRand();
      if(ran>1.) ran=1.;
      q2=-G4Log(1.-ran)/theB3;
      if(q2<0.) q2=0.;
      if(a>6.5) q2=G4Pow::GetInstance()->powA(q2,sevth);
    }
    else
    {
      G4double ran=R4*G4UniformRand();
      if(ran>1.) ran=1.;
      q2=-G4Log(1.-ran)/theB4;
      if(a<6.5) q2=lastTM-q2;                    // u reduced for lightA (starts from 0)
    }
  }
  if(q2<0.) q2=0.;
  if(!(q2>=-1.||q2<=1.))G4cout<<"*NAN*G4QHyElasticCrossSect::GetExchangeT:-t="<<q2<<G4endl;
  if(q2>lastTM)
  {
    q2=lastTM;
  }
  return q2*GeVSQ;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetHMaxT()

G4double G4ChipsHyperonElasticXS::GetHMaxT ( )

private

Definition at line 742 of file G4ChipsHyperonElasticXS.cc.

{
  static const G4double HGeVSQ=gigaelectronvolt*gigaelectronvolt/2.;
  return lastTM*HGeVSQ;
}

Here is the caller graph for this function:

GetIsoCrossSection()

G4double G4ChipsHyperonElasticXS::GetIsoCrossSection ( const G4DynamicParticle *  Pt, G4int  tgZ, G4int  A, const G4Isotope *  iso = 0, const G4Element *  elm = 0, const G4Material *  mat = 0  )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 160 of file G4ChipsHyperonElasticXS.cc.

{
  G4double pMom=Pt->GetTotalMomentum();
  G4int tgN = A - tgZ;
  G4int pdg = Pt->GetDefinition()->GetPDGEncoding();
  
  return GetChipsCrossSection(pMom, tgZ, tgN, pdg);
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetPTables()

G4double G4ChipsHyperonElasticXS::GetPTables ( G4double  lpP, G4double  lPm, G4int  PDG, G4int  tZ, G4int  tN  )

private

Definition at line 366 of file G4ChipsHyperonElasticXS.cc.

{
  // @@ At present all nA==pA ---------> Each neucleus can have not more than 51 parameters
  static const G4double pwd=2727;
  const G4int n_hypel=33;                // #of parameters for pp-elastic (<nPoints=128)
  //                        -0- -1- -2-  -3-  -4- -5-  -6--7--8--9--10--11--12-13--14-
  G4double hyp_el[n_hypel]={1.,.002,.12,.0557,3.5,6.72,99.,2.,3.,5.,74.,3.,3.4,.2,.17,
                            .001,8.,.055,3.64,5.e-5,4000.,1500.,.46,1.2e6,3.5e6,5.e-5,
                            1.e10,8.5e8,1.e10,1.1,3.4e6,6.8e6,0.};
  //                        -15--16- -17- -18- -19-  -20- -21- -22- -23-  -24-   -25-
  //                         -26-  -27- -28-  -29- -30- -31- -32-
  //AR-04Jun2014  if(PDG!=3222 && PDG>3000 && PDG<3335)
  if(PDG>3000 && PDG<3335)
  {
    // -- Total pp elastic cross section cs & s1/b1 (main), s2/b2 (tail1), s3/b3 (tail2) --
    //p2=p*p;p3=p2*p;sp=sqrt(p);p2s=p2*sp;lp=log(p);dl1=lp-(3.=par(3));p4=p2*p2; p=|3-mom|
    //CS=2.865/p2s/(1+.0022/p2s)+(18.9+.6461*dl1*dl1+9./p)/(1.+.425*lp)/(1.+.4276/p4);
    //   par(0)       par(7)     par(1) par(2)      par(4)      par(5)         par(6)
    //dl2=lp-5., s1=(74.+3.*dl2*dl2)/(1+3.4/p4/p)+(.2/p2+17.*p)/(p4+.001*sp),
    //     par(8) par(9) par(10)        par(11)   par(12)par(13)    par(14)
    // b1=8.*p**.055/(1.+3.64/p3); s2=5.e-5+4000./(p4+1500.*p); b2=.46+1.2e6/(p4+3.5e6/sp);
    // par(15) par(16)  par(17)     par(18) par(19)  par(20)   par(21) par(22)  par(23)
    // s3=5.e-5+1.e10/(p4*p4+8.5e8*p2+1.e10); b3=1.1+3.4e6/(p4+6.8e6); ss=0.
    //  par(24) par(25)     par(26)  par(27) par(28) par(29)  par(30)   par(31)
    //
    if(lastPAR[nLast]!=pwd) // A unique flag to avoid the repeatable definition
    {
      if ( tgZ == 1 && tgN == 0 )
      {
        for (G4int ip=0; ip<n_hypel; ip++) lastPAR[ip]=hyp_el[ip]; // Hyperon+P
      }
      else
      {
        G4double a=tgZ+tgN;
        G4double sa=std::sqrt(a);
        G4double ssa=std::sqrt(sa);
        G4double asa=a*sa;
        G4double a2=a*a;
        G4double a3=a2*a;
        G4double a4=a3*a;
        G4double a5=a4*a;
        G4double a6=a4*a2;
        G4double a7=a6*a;
        G4double a8=a7*a;
        G4double a9=a8*a;
        G4double a10=a5*a5;
        G4double a12=a6*a6;
        G4double a14=a7*a7;
        G4double a16=a8*a8;
        G4double a17=a16*a;
        //G4double a20=a16*a4;
        G4double a32=a16*a16;
        // Reaction cross-section parameters (pel=peh_fit.f)
        lastPAR[0]=4./(1.+22/asa);                                           // p1
        lastPAR[1]=2.36*asa/(1.+a*.055/ssa);                                 // p2
        lastPAR[2]=(1.+.00007*a3/ssa)/(1.+.0026*a2);                         // p3
        lastPAR[3]=1.76*a/ssa+.00003*a3;                                     // p4
        lastPAR[4]=(.03+200./a3)/(1.+1.E5/a3/sa);                            // p5
        lastPAR[5]=5.;                                                       // p6
        lastPAR[6]=0.;                                                       // p7 not used
        lastPAR[7]=0.;                                                       // p8 not used
        lastPAR[8]=0.;                                                       // p9 not used
        // @@ the differential cross-section is parameterized separately for A>6 & A<7
        if(a<6.5)
        {
          G4double a28=a16*a12;
          // The main pre-exponent      (pel_sg)
          lastPAR[ 9]=4000*a;                                // p1
          lastPAR[10]=1.2e7*a8+380*a17;                      // p2
          lastPAR[11]=.7/(1.+4.e-12*a16);                    // p3
          lastPAR[12]=2.5/a8/(a4+1.e-16*a32);                // p4
          lastPAR[13]=.28*a;                                 // p5
          lastPAR[14]=1.2*a2+2.3;                            // p6
          lastPAR[15]=3.8/a;                                 // p7
          // The main slope             (pel_sl)
          lastPAR[16]=.01/(1.+.0024*a5);                     // p1
          lastPAR[17]=.2*a;                                  // p2
          lastPAR[18]=9.e-7/(1.+.035*a5);                    // p3
          lastPAR[19]=(42.+2.7e-11*a16)/(1.+.14*a);          // p4
          // The main quadratic         (pel_sh)
          lastPAR[20]=2.25*a3;                               // p1
          lastPAR[21]=18.;                                   // p2
          lastPAR[22]=2.4e-3*a8/(1.+2.6e-4*a7);              // p3
          lastPAR[23]=3.5e-36*a32*a8/(1.+5.e-15*a32/a);      // p4
          // The 1st max pre-exponent   (pel_qq)
          lastPAR[24]=1.e5/(a8+2.5e12/a16);                  // p1
          lastPAR[25]=8.e7/(a12+1.e-27*a28*a28);             // p2 
          lastPAR[26]=.0006*a3;                              // p3
          // The 1st max slope          (pel_qs)
          lastPAR[27]=10.+4.e-8*a12*a;                       // p1
          lastPAR[28]=.114;                                  // p2
          lastPAR[29]=.003;                                  // p3
          lastPAR[30]=2.e-23;                                // p4
          // The effective pre-exponent (pel_ss)
          lastPAR[31]=1./(1.+.0001*a8);                      // p1
          lastPAR[32]=1.5e-4/(1.+5.e-6*a12);                 // p2
          lastPAR[33]=.03;                                   // p3
          // The effective slope        (pel_sb)
          lastPAR[34]=a/2;                                   // p1
          lastPAR[35]=2.e-7*a4;                              // p2
          lastPAR[36]=4.;                                    // p3
          lastPAR[37]=64./a3;                                // p4
          // The gloria pre-exponent    (pel_us)
          lastPAR[38]=1.e8*G4Exp(.32*asa);                // p1
          lastPAR[39]=20.*G4Exp(.45*asa);                 // p2
          lastPAR[40]=7.e3+2.4e6/a5;                         // p3
          lastPAR[41]=2.5e5*G4Exp(.085*a3);               // p4
          lastPAR[42]=2.5*a;                                 // p5
          // The gloria slope           (pel_ub)
          lastPAR[43]=920.+.03*a8*a3;                        // p1
          lastPAR[44]=93.+.0023*a12;                         // p2
        }
        else
        {
          G4double p1a10=2.2e-28*a10;
          G4double r4a16=6.e14/a16;
          G4double s4a16=r4a16*r4a16;
          // a24
          // a36
          // The main pre-exponent      (peh_sg)
          lastPAR[ 9]=4.5*G4Pow::GetInstance()->powA(a,1.15);                  // p1
          lastPAR[10]=.06*G4Pow::GetInstance()->powA(a,.6);                    // p2
          lastPAR[11]=.6*a/(1.+2.e15/a16);                   // p3
          lastPAR[12]=.17/(a+9.e5/a3+1.5e33/a32);            // p4
          lastPAR[13]=(.001+7.e-11*a5)/(1.+4.4e-11*a5);      // p5
          lastPAR[14]=(p1a10*p1a10+2.e-29)/(1.+2.e-22*a12);  // p6
          // The main slope             (peh_sl)
          lastPAR[15]=400./a12+2.e-22*a9;                    // p1
          lastPAR[16]=1.e-32*a12/(1.+5.e22/a14);             // p2
          lastPAR[17]=1000./a2+9.5*sa*ssa;                   // p3
          lastPAR[18]=4.e-6*a*asa+1.e11/a16;                 // p4
          lastPAR[19]=(120./a+.002*a2)/(1.+2.e14/a16);       // p5
          lastPAR[20]=9.+100./a;                             // p6
          // The main quadratic         (peh_sh)
          lastPAR[21]=.002*a3+3.e7/a6;                       // p1
          lastPAR[22]=7.e-15*a4*asa;                         // p2
          lastPAR[23]=9000./a4;                              // p3
          // The 1st max pre-exponent   (peh_qq)
          lastPAR[24]=.0011*asa/(1.+3.e34/a32/a4);           // p1
          lastPAR[25]=1.e-5*a2+2.e14/a16;                    // p2
          lastPAR[26]=1.2e-11*a2/(1.+1.5e19/a12);            // p3
          lastPAR[27]=.016*asa/(1.+5.e16/a16);               // p4
          // The 1st max slope          (peh_qs)
          lastPAR[28]=.002*a4/(1.+7.e7/G4Pow::GetInstance()->powA(a-6.83,14)); // p1
          lastPAR[29]=2.e6/a6+7.2/G4Pow::GetInstance()->powA(a,.11);           // p2
          lastPAR[30]=11.*a3/(1.+7.e23/a16/a8);              // p3
          lastPAR[31]=100./asa;                              // p4
          // The 2nd max pre-exponent   (peh_ss)
          lastPAR[32]=(.1+4.4e-5*a2)/(1.+5.e5/a4);           // p1
          lastPAR[33]=3.5e-4*a2/(1.+1.e8/a8);                // p2
          lastPAR[34]=1.3+3.e5/a4;                           // p3
          lastPAR[35]=500./(a2+50.)+3;                       // p4
          lastPAR[36]=1.e-9/a+s4a16*s4a16;                   // p5
          // The 2nd max slope          (peh_sb)
          lastPAR[37]=.4*asa+3.e-9*a6;                       // p1
          lastPAR[38]=.0005*a5;                              // p2
          lastPAR[39]=.002*a5;                               // p3
          lastPAR[40]=10.;                                   // p4
          // The effective pre-exponent (peh_us)
          lastPAR[41]=.05+.005*a;                            // p1
          lastPAR[42]=7.e-8/sa;                              // p2
          lastPAR[43]=.8*sa;                                 // p3
          lastPAR[44]=.02*sa;                                // p4
          lastPAR[45]=1.e8/a3;                               // p5
          lastPAR[46]=3.e32/(a32+1.e32);                     // p6
          // The effective slope        (peh_ub)
          lastPAR[47]=24.;                                   // p1
          lastPAR[48]=20./sa;                                // p2
          lastPAR[49]=7.e3*a/(sa+1.);                        // p3
          lastPAR[50]=900.*sa/(1.+500./a3);                  // p4
        }
        // Parameter for lowEnergyNeutrons
        lastPAR[51]=1.e15+2.e27/a4/(1.+2.e-18*a16);
      }
      lastPAR[nLast]=pwd;
      // and initialize the zero element of the table
      G4double lp=lPMin;                                      // ln(momentum)
      G4bool memCS=onlyCS;                                    // ??
      onlyCS=false;
      lastCST[0]=GetTabValues(lp, PDG, tgZ, tgN); // Calculate AMDB tables
      onlyCS=memCS;
      lastSST[0]=theSS;
      lastS1T[0]=theS1;
      lastB1T[0]=theB1;
      lastS2T[0]=theS2;
      lastB2T[0]=theB2;
      lastS3T[0]=theS3;
      lastB3T[0]=theB3;
      lastS4T[0]=theS4;
      lastB4T[0]=theB4;
    }
    if(LP>ILP)
    {
      G4int ini = static_cast<int>((ILP-lPMin+.000001)/dlnP)+1; // already inited till this
      if(ini<0) ini=0;
      if(ini<nPoints)
      {
        G4int fin = static_cast<int>((LP-lPMin)/dlnP)+1; // final bin of initialization
        if(fin>=nPoints) fin=nLast;               // Limit of the tabular initialization
        if(fin>=ini)
        {
          G4double lp=0.;
          for(G4int ip=ini; ip<=fin; ip++)        // Calculate tabular CS,S1,B1,S2,B2,S3,B3
          {
            lp=lPMin+ip*dlnP;                     // ln(momentum)
            G4bool memCS=onlyCS;
            onlyCS=false;
            lastCST[ip]=GetTabValues(lp, PDG, tgZ, tgN); // Calculate AMDB tables (ret CS)
            onlyCS=memCS;
            lastSST[ip]=theSS;
            lastS1T[ip]=theS1;
            lastB1T[ip]=theB1;
            lastS2T[ip]=theS2;
            lastB2T[ip]=theB2;
            lastS3T[ip]=theS3;
            lastB3T[ip]=theB3;
            lastS4T[ip]=theS4;
            lastB4T[ip]=theB4;
          }
          return lp;
        }
        else G4cout<<"*Warning*G4ChipsHyperonElasticXS::GetPTables: PDG="<<PDG
                   <<", Z="<<tgZ<<", N="<<tgN<<", i="<<ini<<" > fin="<<fin<<", LP="<<LP
                   <<" > ILP="<<ILP<<" nothing is done!"<<G4endl;
      }
      else G4cout<<"*Warning*G4ChipsHyperonElasticXS::GetPTables: PDG="<<PDG
                 <<", Z="<<tgZ<<", N="<<tgN<<", i="<<ini<<">= max="<<nPoints<<", LP="<<LP
                 <<" > ILP="<<ILP<<", lPMax="<<lPMax<<" nothing is done!"<<G4endl;
    }
  } else {
    // G4cout<<"*Error*G4ChipsHyperonElasticXS::GetPTables: PDG="<<PDG<<", Z="<<tgZ
    //       <<", N="<<tgN<<", while it is defined only for Hyperons"<<G4endl;
    // throw G4QException("G4ChipsHyperonElasticXS::GetPTables:onlyaBA implemented");
    G4ExceptionDescription ed;
    ed << "PDG = " << PDG << ", Z = " << tgZ << ", N = " << tgN
       << ", while it is defined only for Hyperons" << G4endl;
    G4Exception("G4ChipsHyperonElasticXS::GetPTables()", "HAD_CHPS_0000",
                FatalException, ed);
  }
  return ILP;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetQ2max()

G4double G4ChipsHyperonElasticXS::GetQ2max ( G4int  pPDG, G4int  tgZ, G4int  tgN, G4double  pP  )

private

Definition at line 842 of file G4ChipsHyperonElasticXS.cc.

{
  static const G4double mLamb= G4Lambda::Lambda()->GetPDGMass()*.001; // MeV to GeV
  static const G4double mLa2= mLamb*mLamb;
  G4double pP2=pP*pP;                                 // squared momentum of the projectile
  if(tgZ || tgN>-1)                                   // --> Hyperon-A
  {
    G4double mt=G4ParticleTable::GetParticleTable()->GetIonTable()->GetIon(tgZ,tgZ+tgN,0)->GetPDGMass()*.001; // Target mass in GeV

    G4double dmt=mt+mt;
    G4double mds=dmt*std::sqrt(pP2+mLa2)+mLa2+mt*mt;    // Mondelstam mds (@@ other hyperons?)
    return dmt*dmt*pP2/mds;
  }
  else
  {
    // G4cout<<"*Error*G4ChipsHyperonElasticXS::GetQ2ma:PDG="<<PDG<<",Z="<<tgZ<<",N="
    //       <<tgN<<", while it is defined only for p projectiles & Z_target>0"<<G4endl;
    // throw G4QException("G4ChipsHyperonElasticXS::GetQ2max: only HyperA implemented");
    G4ExceptionDescription ed;
    ed << "PDG = " << PDG << ", Z = " << tgZ << ", N = " << tgN
       << ", while it is defined only for p projectiles & Z_target>0" << G4endl;
    G4Exception("G4ChipsHyperonElasticXS::GetQ2max()", "HAD_CHPS_0000",
                FatalException, ed);
    return 0;
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetSlope()

G4double G4ChipsHyperonElasticXS::GetSlope ( G4int  tZ, G4int  tN, G4int  pPDG  )

private

Definition at line 719 of file G4ChipsHyperonElasticXS.cc.

{
  static const G4double GeVSQ=gigaelectronvolt*gigaelectronvolt;
  if(onlyCS)G4cout<<"*Warning*G4ChipsHyperonElasticXS::GetSlope: onlCS=true"<<G4endl;
  if(lastLP<-4.3) return 0.;          // S-wave for p<14 MeV/c (kinE<.1MeV)
  //AR-04Jun2014  if(PDG==3222 || PDG<3000 || PDG>3334)
  if(PDG<3000 || PDG>3334)
  {
    // G4cout<<"*Error*G4ChipsHyperonElasticXS::GetSlope: PDG="<<PDG<<", Z="<<tgZ
    //       <<", N="<<tgN<<", while it is defined only for Hyperons"<<G4endl;
    // throw G4QException("G4ChipsHyperonElasticXS::GetSlope: HypA are implemented");
    G4ExceptionDescription ed;
    ed << "PDG = " << PDG << ", Z = " << tgZ << ", N = " << tgN
       << ", while it is defined only for Hyperons" << G4endl;
    G4Exception("G4ChipsHyperonElasticXS::GetSlope()", "HAD_CHPS_0000",
                FatalException, ed);
  }
  if(theB1<0.) theB1=0.;
  if(!(theB1>=-1.||theB1<=1.)) G4cout<<"*NAN*G4QHyElasticCrossS::Getslope:"<<theB1<<G4endl;
  return theB1/GeVSQ;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetTabValues()

G4double G4ChipsHyperonElasticXS::GetTabValues ( G4double  lp, G4int  pPDG, G4int  tgZ, G4int  tgN  )

private

Definition at line 749 of file G4ChipsHyperonElasticXS.cc.

{
  //AR-04Jun2014  if(PDG==3222 || PDG<3000 || PDG>3334) G4cout<<"*Warning*G4QHypElCS::GTV:P="<<PDG<<G4endl;
  if(PDG<3000 || PDG>3334) G4cout<<"*Warning*G4QHypElCS::GTV:P="<<PDG<<G4endl;
  if(tgZ<0 || tgZ>92)
  {
    G4cout<<"*Warning*G4QHyperonElastCS::GetTabValue:(1-92) NoIsotopesFor Z="<<tgZ<<G4endl;
    return 0.;
  }
  G4int iZ=tgZ-1; // Z index
  if(iZ<0)
  {
    iZ=0;         // conversion of the neutron target to the proton target
    tgZ=1;
    tgN=0;
  }
  G4double p=G4Exp(lp);              // momentum
  G4double sp=std::sqrt(p);             // sqrt(p)
  G4double p2=p*p;            
  G4double p3=p2*p;
  G4double p4=p3*p;
  if ( tgZ == 1 && tgN == 0 ) // Hyperon+P
  {
    G4double dl2=lp-lastPAR[9];
    theSS=lastPAR[32];
    theS1=(lastPAR[10]+lastPAR[11]*dl2*dl2)/(1.+lastPAR[12]/p4/p)+
          (lastPAR[13]/p2+lastPAR[14]*p)/(p4+lastPAR[15]*sp);
    theB1=lastPAR[16]*G4Pow::GetInstance()->powA(p,lastPAR[17])/(1.+lastPAR[18]/p3);
    theS2=lastPAR[19]+lastPAR[20]/(p4+lastPAR[21]*p);
    theB2=lastPAR[22]+lastPAR[23]/(p4+lastPAR[24]/sp); 
    theS3=lastPAR[25]+lastPAR[26]/(p4*p4+lastPAR[27]*p2+lastPAR[28]);
    theB3=lastPAR[29]+lastPAR[30]/(p4+lastPAR[31]); 
    theS4=0.;
    theB4=0.; 
    // Returns the total elastic pim-p cross-section (to avoid spoiling lastSIG)
    G4double dp=lp-lastPAR[4];
    return lastPAR[0]/(lastPAR[1]+p2*(lastPAR[2]+p2))+(lastPAR[3]*dp*dp+lastPAR[5]+
           lastPAR[6]/p2)/(1.+lastPAR[7]/sp+lastPAR[8]/p4);
  }
  else
  {
    G4double p5=p4*p;
    G4double p6=p5*p;
    G4double p8=p6*p2;
    G4double p10=p8*p2;
    G4double p12=p10*p2;
    G4double p16=p8*p8;
    //G4double p24=p16*p8;
    G4double dl=lp-5.;
    G4double a=tgZ+tgN;
    G4double pah=G4Pow::GetInstance()->powA(p,a/2);
    G4double pa=pah*pah;
    G4double pa2=pa*pa;
    if(a<6.5)
    {
      theS1=lastPAR[9]/(1.+lastPAR[10]*p4*pa)+lastPAR[11]/(p4+lastPAR[12]*p4/pa2)+
            (lastPAR[13]*dl*dl+lastPAR[14])/(1.+lastPAR[15]/p2);
      theB1=(lastPAR[16]+lastPAR[17]*p2)/(p4+lastPAR[18]/pah)+lastPAR[19];
      theSS=lastPAR[20]/(1.+lastPAR[21]/p2)+lastPAR[22]/(p6/pa+lastPAR[23]/p16);
      theS2=lastPAR[24]/(pa/p2+lastPAR[25]/p4)+lastPAR[26];
      theB2=lastPAR[27]*G4Pow::GetInstance()->powA(p,lastPAR[28])+lastPAR[29]/(p8+lastPAR[30]/p16);
      theS3=lastPAR[31]/(pa*p+lastPAR[32]/pa)+lastPAR[33];
      theB3=lastPAR[34]/(p3+lastPAR[35]/p6)+lastPAR[36]/(1.+lastPAR[37]/p2);
      theS4=p2*(pah*lastPAR[38]*G4Exp(-pah*lastPAR[39])+
                lastPAR[40]/(1.+lastPAR[41]*G4Pow::GetInstance()->powA(p,lastPAR[42])));
      theB4=lastPAR[43]*pa/p2/(1.+pa*lastPAR[44]);
    }
    else
    {
      theS1=lastPAR[9]/(1.+lastPAR[10]/p4)+lastPAR[11]/(p4+lastPAR[12]/p2)+
            lastPAR[13]/(p5+lastPAR[14]/p16);
      theB1=(lastPAR[15]/p8+lastPAR[19])/(p+lastPAR[16]/G4Pow::GetInstance()->powA(p,lastPAR[20]))+
            lastPAR[17]/(1.+lastPAR[18]/p4);
      theSS=lastPAR[21]/(p4/G4Pow::GetInstance()->powA(p,lastPAR[23])+lastPAR[22]/p4);
      theS2=lastPAR[24]/p4/(G4Pow::GetInstance()->powA(p,lastPAR[25])+lastPAR[26]/p12)+lastPAR[27];
      theB2=lastPAR[28]/G4Pow::GetInstance()->powA(p,lastPAR[29])+lastPAR[30]/G4Pow::GetInstance()->powA(p,lastPAR[31]);
      theS3=lastPAR[32]/G4Pow::GetInstance()->powA(p,lastPAR[35])/(1.+lastPAR[36]/p12)+
            lastPAR[33]/(1.+lastPAR[34]/p6);
      theB3=lastPAR[37]/p8+lastPAR[38]/p2+lastPAR[39]/(1.+lastPAR[40]/p8);
      theS4=(lastPAR[41]/p4+lastPAR[46]/p)/(1.+lastPAR[42]/p10)+
            (lastPAR[43]+lastPAR[44]*dl*dl)/(1.+lastPAR[45]/p12);
      theB4=lastPAR[47]/(1.+lastPAR[48]/p)+lastPAR[49]*p4/(1.+lastPAR[50]*p5);
    }
    // Returns the total elastic (n/p)A cross-section (to avoid spoiling lastSIG)
    G4double dlp=lp-lastPAR[5]; // ax
    //         p1               p2          p3                 p4               p5
    return (lastPAR[0]*dlp*dlp+lastPAR[1])/(1.+lastPAR[2]/p)+lastPAR[3]/(p3+lastPAR[4]);
  }
  return 0.;
} // End of GetTableValues

Here is the call graph for this function:

Here is the caller graph for this function:

IsIsoApplicable()

G4bool G4ChipsHyperonElasticXS::IsIsoApplicable ( const G4DynamicParticle *  Pt, G4int  Z, G4int  A, const G4Element *  elm, const G4Material *  mat  )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 151 of file G4ChipsHyperonElasticXS.cc.

{
  return true;
}

Here is the caller graph for this function:

Member Data Documentation

B1T

std::vector<G4double*> G4ChipsHyperonElasticXS::B1T

private

Definition at line 130 of file G4ChipsHyperonElasticXS.hh.

B2T

std::vector<G4double*> G4ChipsHyperonElasticXS::B2T

private

Definition at line 132 of file G4ChipsHyperonElasticXS.hh.

B3T

std::vector<G4double*> G4ChipsHyperonElasticXS::B3T

private

Definition at line 134 of file G4ChipsHyperonElasticXS.hh.

B4T

std::vector<G4double*> G4ChipsHyperonElasticXS::B4T

private

Definition at line 136 of file G4ChipsHyperonElasticXS.hh.

colCS

std::vector<G4double> G4ChipsHyperonElasticXS::colCS

private

Definition at line 142 of file G4ChipsHyperonElasticXS.hh.

colN

std::vector<G4int> G4ChipsHyperonElasticXS::colN

private

Definition at line 138 of file G4ChipsHyperonElasticXS.hh.

colP

std::vector<G4double> G4ChipsHyperonElasticXS::colP

private

Definition at line 140 of file G4ChipsHyperonElasticXS.hh.

colTH

std::vector<G4double> G4ChipsHyperonElasticXS::colTH

private

Definition at line 141 of file G4ChipsHyperonElasticXS.hh.

colZ

std::vector<G4int> G4ChipsHyperonElasticXS::colZ

private

Definition at line 139 of file G4ChipsHyperonElasticXS.hh.

CST

std::vector<G4double*> G4ChipsHyperonElasticXS::CST

private

Definition at line 127 of file G4ChipsHyperonElasticXS.hh.

dlnP

G4double G4ChipsHyperonElasticXS::dlnP

private

Definition at line 89 of file G4ChipsHyperonElasticXS.hh.

lastB1T

G4double* G4ChipsHyperonElasticXS::lastB1T

private

Definition at line 118 of file G4ChipsHyperonElasticXS.hh.

lastB2T

G4double* G4ChipsHyperonElasticXS::lastB2T

private

Definition at line 120 of file G4ChipsHyperonElasticXS.hh.

lastB3T

G4double* G4ChipsHyperonElasticXS::lastB3T

private

Definition at line 122 of file G4ChipsHyperonElasticXS.hh.

lastB4T

G4double* G4ChipsHyperonElasticXS::lastB4T

private

Definition at line 124 of file G4ChipsHyperonElasticXS.hh.

lastCS

G4double G4ChipsHyperonElasticXS::lastCS

private

Definition at line 99 of file G4ChipsHyperonElasticXS.hh.

lastCST

G4double* G4ChipsHyperonElasticXS::lastCST

private

Definition at line 114 of file G4ChipsHyperonElasticXS.hh.

lastI

G4int G4ChipsHyperonElasticXS::lastI

private

Definition at line 100 of file G4ChipsHyperonElasticXS.hh.

lastLP

G4double G4ChipsHyperonElasticXS::lastLP

private

Definition at line 93 of file G4ChipsHyperonElasticXS.hh.

lastN

G4int G4ChipsHyperonElasticXS::lastN

private

Definition at line 95 of file G4ChipsHyperonElasticXS.hh.

lastP

G4double G4ChipsHyperonElasticXS::lastP

private

Definition at line 97 of file G4ChipsHyperonElasticXS.hh.

lastPAR

G4double* G4ChipsHyperonElasticXS::lastPAR

private

Definition at line 115 of file G4ChipsHyperonElasticXS.hh.

lastPIN

G4double G4ChipsHyperonElasticXS::lastPIN

private

Definition at line 113 of file G4ChipsHyperonElasticXS.hh.

lastS1T

G4double* G4ChipsHyperonElasticXS::lastS1T

private

Definition at line 117 of file G4ChipsHyperonElasticXS.hh.

lastS2T

G4double* G4ChipsHyperonElasticXS::lastS2T

private

Definition at line 119 of file G4ChipsHyperonElasticXS.hh.

lastS3T

G4double* G4ChipsHyperonElasticXS::lastS3T

private

Definition at line 121 of file G4ChipsHyperonElasticXS.hh.

lastS4T

G4double* G4ChipsHyperonElasticXS::lastS4T

private

Definition at line 123 of file G4ChipsHyperonElasticXS.hh.

lastSIG

G4double G4ChipsHyperonElasticXS::lastSIG

private

Definition at line 92 of file G4ChipsHyperonElasticXS.hh.

lastSST

G4double* G4ChipsHyperonElasticXS::lastSST

private

Definition at line 116 of file G4ChipsHyperonElasticXS.hh.

lastTH

G4double G4ChipsHyperonElasticXS::lastTH

private

Definition at line 98 of file G4ChipsHyperonElasticXS.hh.

lastTM

G4double G4ChipsHyperonElasticXS::lastTM

private

Definition at line 94 of file G4ChipsHyperonElasticXS.hh.

lastTN

G4int G4ChipsHyperonElasticXS::lastTN

private

Definition at line 112 of file G4ChipsHyperonElasticXS.hh.

lastTZ

G4int G4ChipsHyperonElasticXS::lastTZ

private

Definition at line 111 of file G4ChipsHyperonElasticXS.hh.

lastZ

G4int G4ChipsHyperonElasticXS::lastZ

private

Definition at line 96 of file G4ChipsHyperonElasticXS.hh.

lPMax

G4double G4ChipsHyperonElasticXS::lPMax

private

Definition at line 88 of file G4ChipsHyperonElasticXS.hh.

lPMin

G4double G4ChipsHyperonElasticXS::lPMin

private

Definition at line 87 of file G4ChipsHyperonElasticXS.hh.

nLast

const G4int G4ChipsHyperonElasticXS::nLast

private

Definition at line 86 of file G4ChipsHyperonElasticXS.hh.

nPoints

const G4int G4ChipsHyperonElasticXS::nPoints

private

Definition at line 85 of file G4ChipsHyperonElasticXS.hh.

onlyCS

G4bool G4ChipsHyperonElasticXS::onlyCS

private

Definition at line 91 of file G4ChipsHyperonElasticXS.hh.

PAR

std::vector<G4double*> G4ChipsHyperonElasticXS::PAR

private

Definition at line 126 of file G4ChipsHyperonElasticXS.hh.

PIN

std::vector<G4double> G4ChipsHyperonElasticXS::PIN

private

Definition at line 144 of file G4ChipsHyperonElasticXS.hh.

S1T

std::vector<G4double*> G4ChipsHyperonElasticXS::S1T

private

Definition at line 129 of file G4ChipsHyperonElasticXS.hh.

S2T

std::vector<G4double*> G4ChipsHyperonElasticXS::S2T

private

Definition at line 131 of file G4ChipsHyperonElasticXS.hh.

S3T

std::vector<G4double*> G4ChipsHyperonElasticXS::S3T

private

Definition at line 133 of file G4ChipsHyperonElasticXS.hh.

S4T

std::vector<G4double*> G4ChipsHyperonElasticXS::S4T

private

Definition at line 135 of file G4ChipsHyperonElasticXS.hh.

SST

std::vector<G4double*> G4ChipsHyperonElasticXS::SST

private

Definition at line 128 of file G4ChipsHyperonElasticXS.hh.

theB1

G4double G4ChipsHyperonElasticXS::theB1

private

Definition at line 103 of file G4ChipsHyperonElasticXS.hh.

theB2

G4double G4ChipsHyperonElasticXS::theB2

private

Definition at line 105 of file G4ChipsHyperonElasticXS.hh.

theB3

G4double G4ChipsHyperonElasticXS::theB3

private

Definition at line 107 of file G4ChipsHyperonElasticXS.hh.

theB4

G4double G4ChipsHyperonElasticXS::theB4

private

Definition at line 109 of file G4ChipsHyperonElasticXS.hh.

theS1

G4double G4ChipsHyperonElasticXS::theS1

private

Definition at line 102 of file G4ChipsHyperonElasticXS.hh.

theS2

G4double G4ChipsHyperonElasticXS::theS2

private

Definition at line 104 of file G4ChipsHyperonElasticXS.hh.

theS3

G4double G4ChipsHyperonElasticXS::theS3

private

Definition at line 106 of file G4ChipsHyperonElasticXS.hh.

theS4

G4double G4ChipsHyperonElasticXS::theS4

private

Definition at line 108 of file G4ChipsHyperonElasticXS.hh.

theSS

G4double G4ChipsHyperonElasticXS::theSS

private

Definition at line 101 of file G4ChipsHyperonElasticXS.hh.

---

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

• G4ChipsHyperonElasticXS.hh
• G4ChipsHyperonElasticXS.cc