G4ChipsKaonMinusInelasticXS Class Reference

#include <G4ChipsKaonMinusInelasticXS.hh>

Inheritance diagram for G4ChipsKaonMinusInelasticXS:

Collaboration diagram for G4ChipsKaonMinusInelasticXS:

Public Member Functions
	G4ChipsKaonMinusInelasticXS ()
	~G4ChipsKaonMinusInelasticXS ()
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)
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 (G4int F, G4int I, G4int PDG, G4int Z, G4int N, G4double Momentum)
G4double	CrossSectionLin (G4int targZ, G4int targN, G4double P)
G4double	CrossSectionLog (G4int targZ, G4int targN, G4double lP)
G4double	CrossSectionFormula (G4int targZ, G4int targN, G4double P, G4double lP)
G4double	EquLinearFit (G4double X, G4int N, G4double X0, G4double DX, G4double *Y)

Private Attributes
G4double *	lastLEN
G4double *	lastHEN
G4int	lastN
G4int	lastZ
G4double	lastP
G4double	lastTH
G4double	lastCS
G4int	lastI
std::vector< G4double * > *	LEN
std::vector< G4double * > *	HEN
G4int	j
std::vector< G4int >	colN
std::vector< G4int >	colZ
std::vector< G4double >	colP
std::vector< G4double >	colTH
std::vector< G4double >	colCS

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 45 of file G4ChipsKaonMinusInelasticXS.hh.

Constructor & Destructor Documentation

G4ChipsKaonMinusInelasticXS()

G4ChipsKaonMinusInelasticXS::G4ChipsKaonMinusInelasticXS

(

)

Definition at line 66 of file G4ChipsKaonMinusInelasticXS.cc.

                                                        :G4VCrossSectionDataSet(Default_Name())
{
  lastLEN=0; // Pointer to lastArray of LowEn CS
  lastHEN=0; // Pointer to lastArray of HighEn CS
  lastN=0;   // The last N of calculated nucleus
  lastZ=0;   // The last Z of calculated nucleus
  lastP=0.;  // Last used in CrossSection Momentum
  lastTH=0.; // Last threshold momentum
  lastCS=0.; // Last value of the Cross Section
  lastI=0;   // The last position in the DAMDB
  LEN = new std::vector<G4double*>;
  HEN = new std::vector<G4double*>;
}

~G4ChipsKaonMinusInelasticXS()

G4ChipsKaonMinusInelasticXS::~G4ChipsKaonMinusInelasticXS

(

)

Definition at line 80 of file G4ChipsKaonMinusInelasticXS.cc.

{  
  G4int lens=LEN->size();
  for(G4int i=0; i<lens; ++i) delete[] (*LEN)[i];
  delete LEN;
 
  G4int hens=HEN->size();
  for(G4int i=0; i<hens; ++i) delete[] (*HEN)[i];
  delete HEN; 
}

Member Function Documentation

CalculateCrossSection()

G4double G4ChipsKaonMinusInelasticXS::CalculateCrossSection ( G4int  F, G4int  I, G4int  PDG, G4int  Z, G4int  N, G4double  Momentum  )

private

Definition at line 193 of file G4ChipsKaonMinusInelasticXS.cc.

{
  G4double sigma=0.;
  if(F&&I) sigma=0.;                   // @@ *!* Fake line *!* to use F & I !!!Temporary!!!
  //G4double A=targN+targZ;              // A of the target
  if(F<=0)                             // This isotope was not the last used isotop
  {
    if(F<0)                            // This isotope was found in DAMDB =-----=> RETRIEVE
    {
      G4int sync=LEN->size();
      if(sync<=I) G4cerr<<"*!*G4QPiMinusNuclCS::CalcCrosSect:Sync="<<sync<<"<="<<I<<G4endl;
      lastLEN=(*LEN)[I];                // Pointer to prepared LowEnergy cross sections
      lastHEN=(*HEN)[I];                // Pointer to prepared High Energy cross sections
    }
    else                               // This isotope wasn't calculated before => CREATE
    {
      lastLEN = new G4double[nL];      // Allocate memory for the new LEN cross sections
      lastHEN = new G4double[nH];      // Allocate memory for the new HEN cross sections
      // --- Instead of making a separate function ---
      G4double P=THmiG;                // Table threshold in GeV/c
      for(G4int k=0; k<nL; k++)
      {
        lastLEN[k] = CrossSectionLin(targZ, targN, P);
        P+=dPG;
      }
      G4double lP=milPG;
      for(G4int n=0; n<nH; n++)
      {
        lastHEN[n] = CrossSectionLog(targZ, targN, lP);
        lP+=dlP;
      }
      // --- End of possible separate function
      // *** The synchronization check ***
      G4int sync=LEN->size();
      if(sync!=I)
      {
        G4cerr<<"***G4ChipsKaonMinusCS::CalcCrossSect: Sinc="<<sync<<"#"<<I<<", Z=" <<targZ
              <<", N="<<targN<<", F="<<F<<G4endl;
        //G4Exception("G4PiMinusNuclearCS::CalculateCS:","39",FatalException,"DBoverflow");
      }
      LEN->push_back(lastLEN);         // remember the Low Energy Table
      HEN->push_back(lastHEN);         // remember the High Energy Table
    } // End of creation of the new set of parameters
  } // End of parameters udate
  // =------------------= NOW the Magic Formula =--------------------------=
  if (Momentum<lastTH) return 0.;      // It must be already checked in the interface class
  else if (Momentum<Pmin)              // High Energy region
  {
    sigma=EquLinearFit(Momentum,nL,THmin,dP,lastLEN);
  }
  else if (Momentum<Pmax)              // High Energy region
  {
    G4double lP=std::log(Momentum);
    sigma=EquLinearFit(lP,nH,milP,dlP,lastHEN);
  }
  else                                 // UHE region (calculation, not frequent)
  {
    G4double P=0.001*Momentum;         // Approximation formula is for P in GeV/c
    sigma=CrossSectionFormula(targZ, targN, P, std::log(P));
  }
  if(sigma<0.) return 0.;
  return sigma;
}

Here is the call graph for this function:

Here is the caller graph for this function:

CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 92 of file G4ChipsKaonMinusInelasticXS.cc.

{
    outFile << "G4ChipsKaonMinusInelasticXS provides the inelastic cross\n"
            << "section for K- 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:

CrossSectionFormula()

G4double G4ChipsKaonMinusInelasticXS::CrossSectionFormula ( G4int  targZ, G4int  targN, G4double  P, G4double  lP  )

private

Definition at line 272 of file G4ChipsKaonMinusInelasticXS.cc.

{
  G4double sigma=0.;
  if(tZ==1 && !tN)                        // PiMin-Proton interaction from G4QuasiElRatios
  {
    G4double ld=lP-3.5;
    G4double ld2=ld*ld;
    G4double p2=P*P;
    G4double p4=p2*p2;
    G4double sp=std::sqrt(P);
    G4double psp=P*sp;
    G4double lm=P-.39;
    G4double md=lm*lm+.000156;
    G4double lh=P-1.;
    G4double hd=lh*lh+.0156;
    G4double El=(.0557*ld2+2.23)/(1.-.7/sp+.075/p4);
    G4double To=(.3*ld2+19.5)/(1.-.21/sp+.52/p4);
    sigma=8.8/psp+(To-El)+.002/md+.15/hd;
  }
  else if(tZ==1 && tN==1)                  // kmp_tot
  {
    G4double p2=P*P;
    G4double dX=lP-3.7;
    G4double dR=P-.94;
    G4double sp=std::sqrt(P);
    sigma=(.6*dX*dX+36.)/(1.-.11/sp+.52/p2/p2)+.7/(dR*dR+.0256)+18./P/sp;
  }
  else if(tZ<97 && tN<152)                // General solution
  {
    G4double d=lP-4.2;
    G4double sp=std::sqrt(P);
    G4double p2=P*P;
    G4double a=tN+tZ;                       // A of the target
    G4double sa=std::sqrt(a);
    G4double al=std::log(a);
    G4double a2=a*a;
    G4double c=52.*std::exp(al*0.6)*(1.+97./a2)/(1.+9.8/a)/(1.+47./a2);
    G4double gg=-.2-.003*a;
    G4double h=.5+.07*a;
    G4double v=P-1.;
    G4double f=.6*a*sa/(1.+.00002*a2);
    G4double u=.125+.127*al;
    sigma=(c+d*d)/(1.+gg/sp+h/p2/p2)+f/(v*v+u*u)+20.*sa/P/sp;
  }
  else
  {
    G4cerr<<"-Warning-G4ChipsKMinusNuclearCroSect::CSForm:*Bad A* Z="<<tZ<<", N="<<tN<<G4endl;
    sigma=0.;
  }
  if(sigma<0.) return 0.;
  return sigma;  
}

Here is the caller graph for this function:

CrossSectionLin()

G4double G4ChipsKaonMinusInelasticXS::CrossSectionLin ( G4int  targZ, G4int  targN, G4double  P  )

private

Definition at line 259 of file G4ChipsKaonMinusInelasticXS.cc.

{
  G4double lP=std::log(P);
  return CrossSectionFormula(tZ, tN, P, lP);
}

Here is the call graph for this function:

Here is the caller graph for this function:

CrossSectionLog()

G4double G4ChipsKaonMinusInelasticXS::CrossSectionLog ( G4int  targZ, G4int  targN, G4double  lP  )

private

Definition at line 266 of file G4ChipsKaonMinusInelasticXS.cc.

{
  G4double P=std::exp(lP);
  return CrossSectionFormula(tZ, tN, P, lP);
}

Here is the call graph for this function:

Here is the caller graph for this function:

Default_Name()

static const char* G4ChipsKaonMinusInelasticXS::Default_Name ( )

inlinestatic

Definition at line 54 of file G4ChipsKaonMinusInelasticXS.hh.

{return "ChipsKaonMinusInelasticXS";}

Here is the call graph for this function:

Here is the caller graph for this function:

EquLinearFit()

G4double G4ChipsKaonMinusInelasticXS::EquLinearFit ( G4double  X, G4int  N, G4double  X0, G4double  DX, G4double *  Y  )

private

Definition at line 326 of file G4ChipsKaonMinusInelasticXS.cc.

{
  if(DX<=0. || N<2)
    {
      G4cerr<<"***G4ChipsKaonMinusInelasticXS::EquLinearFit: DX="<<DX<<", N="<<N<<G4endl;
      return Y[0];
    }
  
  G4int    N2=N-2;
  G4double d=(X-X0)/DX;
  G4int         jj=static_cast<int>(d);
  if     (jj<0)  jj=0;
  else if(jj>N2) jj=N2;
  d-=jj; // excess
  G4double yi=Y[jj];
  G4double sigma=yi+(Y[jj+1]-yi)*d;
  
  return sigma;
}

Here is the caller graph for this function:

GetChipsCrossSection()

G4double G4ChipsKaonMinusInelasticXS::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 121 of file G4ChipsKaonMinusInelasticXS.cc.

{
  G4bool in=false;                     // By default the isotope must be found in the AMDB
  if(tgN!=lastN || tgZ!=lastZ)         // The nucleus was not the last used isotope
  {
    in = false;                        // By default the isotope haven't be found in 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
    j  = 0;                            // A#0f records found in DB for this projectile
    if(lastI) for(G4int i=0; i<lastI; i++) // AMDB exists, try to find the (Z,N) isotope
    {
      if(colN[i]==tgN && colZ[i]==tgZ) // Try the record "i" in the AMDB
      {
        lastI=i;                       // Remember the index for future fast/last use
        lastTH =colTH[i];              // The last THreshold (A-dependent)
        if(pMom<=lastTH)
        {
          return 0.;                   // Energy is below the Threshold value
        }
        lastP  =colP [i];              // Last Momentum  (A-dependent)
        lastCS =colCS[i];              // Last CrossSect (A-dependent)
        in = true;                     // This is the case when the isotop is found in DB
        // Momentum pMom is in IU ! @@ Units
        lastCS=CalculateCrossSection(-1,j,-321,lastZ,lastN,pMom); // read & update
        if(lastCS<=0. && pMom>lastTH)  // Correct the threshold (@@ No intermediate Zeros)
        {
          lastCS=0.;
          lastTH=pMom;
        }
        break;                         // Go out of the LOOP
      }
      j++;                             // Increment a#0f records found in DB
    }
    if(!in)                            // This isotope has not been calculated previously
    {
      lastCS=CalculateCrossSection(0,j,-321,lastZ,lastN,pMom); //calculate & create
      //if(lastCS>0.)                   // It means that the AMBD was initialized
      //{

      //        lastTH = ThresholdEnergy(tgZ, tgN); // The Threshold Energy which is now the last

      lastTH = 0; // WP - to be checked!!!
        colN.push_back(tgN);
        colZ.push_back(tgZ);
        colP.push_back(pMom);
        colTH.push_back(lastTH);
        colCS.push_back(lastCS);
      //} // M.K. Presence of H1 with high threshold breaks the syncronization
      return lastCS*millibarn;
    } // End of creation of the new set of parameters
    else
    {
      colP[lastI]=pMom;
      colCS[lastI]=lastCS;
    }
  } // End of parameters udate
  else if(pMom<=lastTH)
  {
    return 0.;                         // Momentum is below the Threshold Value -> CS=0
  }
  else                                 // It is the last used -> use the current tables
  {
    lastCS=CalculateCrossSection(1,j,-321,lastZ,lastN,pMom); // Only read and UpdateDB
    lastP=pMom;
  }
  return lastCS*millibarn;
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetIsoCrossSection()

G4double G4ChipsKaonMinusInelasticXS::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 110 of file G4ChipsKaonMinusInelasticXS.cc.

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

Here is the call graph for this function:

Here is the caller graph for this function:

IsIsoApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 100 of file G4ChipsKaonMinusInelasticXS.cc.

{
  return true;
}

Here is the caller graph for this function:

Member Data Documentation

colCS

std::vector<G4double> G4ChipsKaonMinusInelasticXS::colCS

private

Definition at line 97 of file G4ChipsKaonMinusInelasticXS.hh.

colN

std::vector<G4int> G4ChipsKaonMinusInelasticXS::colN

private

Definition at line 93 of file G4ChipsKaonMinusInelasticXS.hh.

colP

std::vector<G4double> G4ChipsKaonMinusInelasticXS::colP

private

Definition at line 95 of file G4ChipsKaonMinusInelasticXS.hh.

colTH

std::vector<G4double> G4ChipsKaonMinusInelasticXS::colTH

private

Definition at line 96 of file G4ChipsKaonMinusInelasticXS.hh.

colZ

std::vector<G4int> G4ChipsKaonMinusInelasticXS::colZ

private

Definition at line 94 of file G4ChipsKaonMinusInelasticXS.hh.

HEN

std::vector<G4double*>* G4ChipsKaonMinusInelasticXS::HEN

private

Definition at line 90 of file G4ChipsKaonMinusInelasticXS.hh.

j

G4int G4ChipsKaonMinusInelasticXS::j

private

Definition at line 92 of file G4ChipsKaonMinusInelasticXS.hh.

lastCS

G4double G4ChipsKaonMinusInelasticXS::lastCS

private

Definition at line 87 of file G4ChipsKaonMinusInelasticXS.hh.

lastHEN

G4double* G4ChipsKaonMinusInelasticXS::lastHEN

private

Definition at line 82 of file G4ChipsKaonMinusInelasticXS.hh.

lastI

G4int G4ChipsKaonMinusInelasticXS::lastI

private

Definition at line 88 of file G4ChipsKaonMinusInelasticXS.hh.

lastLEN

G4double* G4ChipsKaonMinusInelasticXS::lastLEN

private

Definition at line 81 of file G4ChipsKaonMinusInelasticXS.hh.

lastN

G4int G4ChipsKaonMinusInelasticXS::lastN

private

Definition at line 83 of file G4ChipsKaonMinusInelasticXS.hh.

lastP

G4double G4ChipsKaonMinusInelasticXS::lastP

private

Definition at line 85 of file G4ChipsKaonMinusInelasticXS.hh.

lastTH

G4double G4ChipsKaonMinusInelasticXS::lastTH

private

Definition at line 86 of file G4ChipsKaonMinusInelasticXS.hh.

lastZ

G4int G4ChipsKaonMinusInelasticXS::lastZ

private

Definition at line 84 of file G4ChipsKaonMinusInelasticXS.hh.

LEN

std::vector<G4double*>* G4ChipsKaonMinusInelasticXS::LEN

private

Definition at line 89 of file G4ChipsKaonMinusInelasticXS.hh.

---

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

• G4ChipsKaonMinusInelasticXS.hh
• G4ChipsKaonMinusInelasticXS.cc