#include <G4ChipsKaonMinusInelasticXS.hh>

Inheritance diagram for G4ChipsKaonMinusInelasticXS:

Public Member Functions
	G4ChipsKaonMinusInelasticXS ()

	~G4ChipsKaonMinusInelasticXS ()

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)

virtual G4bool	IsIsoApplicable (const G4DynamicParticle , G4int Z, G4int A, const G4Element elm=0, 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 G4double	GetIsoCrossSection (const G4DynamicParticle , G4int Z, G4int A, const G4Isotope iso=0, const G4Element elm=0, const G4Material mat=0)

virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual void	CrossSectionDescription (std::ostream &) const

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 ()

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

Constructor & Destructor Documentation

◆ G4ChipsKaonMinusInelasticXS()

G4ChipsKaonMinusInelasticXS::G4ChipsKaonMinusInelasticXS ( )

Definition at line 53 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 67 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

◆ Default_Name()

static const char * G4ChipsKaonMinusInelasticXS::Default_Name ( )

inlinestatic

Definition at line 56 of file G4ChipsKaonMinusInelasticXS.hh.

56{return "ChipsKaonMinusInelasticXS";}

Referenced by G4FTFBinaryKaonBuilder::Build(), G4QGSBinaryPiKBuilder::Build(), G4QGSPPiKBuilder::Build(), G4HadronInelasticQBBC::ConstructProcess(), HadronPhysicsFTFP_BERT::ConstructProcess(), HadronPhysicsFTFP_BERT_HP::ConstructProcess(), HadronPhysicsFTFP_BERT_TRV::ConstructProcess(), HadronPhysicsQGSP_FTFP_BERT_95::ConstructProcess(), HadronPhysicsShielding::ConstructProcess(), G4ChipsComponentXS::G4ChipsComponentXS(), and G4ChipsKaonZeroInelasticXS::G4ChipsKaonZeroInelasticXS().

◆ 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 101 of file G4ChipsKaonMinusInelasticXS.cc.

{
  static G4int j;                      // A#0f Z/N-records already tested in AMDB
  static std::vector <G4int>    colN;  // Vector of N for calculated nuclei (isotops)
  static std::vector <G4int>    colZ;  // Vector of Z for calculated nuclei (isotops)
  static std::vector <G4double> colP;  // Vector of last momenta for the reaction
  static std::vector <G4double> colTH; // Vector of energy thresholds for the reaction
  static std::vector <G4double> colCS; // Vector of last cross sections for the reaction
  // ***---*** End of the mandatory Static Definitions of the Associative Memory ***---***
 
  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
    {
      //!!The slave functions must provide cross-sections in millibarns (mb) !! (not in IU)
      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;
}

Referenced by G4ChipsKaonZeroInelasticXS::GetChipsCrossSection(), G4ChipsComponentXS::GetInelasticElementCrossSection(), GetIsoCrossSection(), and G4ChipsComponentXS::GetTotalElementCrossSection().

◆ 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 90 of file G4ChipsKaonMinusInelasticXS.cc.

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

◆ IsIsoApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 78 of file G4ChipsKaonMinusInelasticXS.cc.

{
  G4ParticleDefinition* particle = Pt->GetDefinition();
  if (particle ==       G4KaonMinus::KaonMinus()      ) return true;
  return false;
}

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

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4ChipsKaonMinusInelasticXS()

◆ ~G4ChipsKaonMinusInelasticXS()

Member Function Documentation

◆ Default_Name()

◆ GetChipsCrossSection()

◆ GetIsoCrossSection()

◆ IsIsoApplicable()