#include <G4ChipsAntiBaryonInelasticXS.hh>

Inheritance diagram for G4ChipsAntiBaryonInelasticXS:

Collaboration diagram for G4ChipsAntiBaryonInelasticXS:

Public Member Functions
	G4ChipsAntiBaryonInelasticXS ()

	~G4ChipsAntiBaryonInelasticXS ()

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 44 of file G4ChipsAntiBaryonInelasticXS.hh.

Constructor & Destructor Documentation

◆ G4ChipsAntiBaryonInelasticXS()

G4ChipsAntiBaryonInelasticXS::G4ChipsAntiBaryonInelasticXS ( )

Definition at line 62 of file G4ChipsAntiBaryonInelasticXS.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 Cross Section 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*>;
 }

◆ ~G4ChipsAntiBaryonInelasticXS()

G4ChipsAntiBaryonInelasticXS::~G4ChipsAntiBaryonInelasticXS ( )

Definition at line 76 of file G4ChipsAntiBaryonInelasticXS.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 G4ChipsAntiBaryonInelasticXS::CalculateCrossSection	(	G4int	F,
		G4int	I,
		G4int	PDG,
		G4int	Z,
		G4int	N,
		G4double	Momentum
	)

private

Definition at line 227 of file G4ChipsAntiBaryonInelasticXS.cc.

 {
   static const G4double THmin=27.;     // default minimum Momentum (MeV/c) Threshold
   static const G4double THmiG=THmin*.001; // minimum Momentum (GeV/c) Threshold
   static const G4double dP=10.;        // step for the LEN (Low ENergy) table MeV/c
   static const G4double dPG=dP*.001;   // step for the LEN (Low ENergy) table GeV/c
   static const G4int    nL=105;        // A#of LEN points in E (step 10 MeV/c)
   static const G4double Pmin=THmin+(nL-1)*dP; // minP for the HighE part with safety
   static const G4double Pmax=227000.;  // maxP for the HEN (High ENergy) part 227 GeV
   static const G4int    nH=224;        // A#of HEN points in lnE
   static const G4double milP=G4Log(Pmin);// Low logarithm energy for the HEN part
   static const G4double malP=G4Log(Pmax);// High logarithm energy (each 2.75 percent)
   static const G4double dlP=(malP-milP)/(nH-1); // Step in log energy in the HEN part
   static const G4double milPG=G4Log(.001*Pmin);// Low logarithmEnergy for HEN part GeV/c
   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<<"***G4QPiMinusNuclCS::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=G4Log(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, G4Log(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 G4ChipsAntiBaryonInelasticXS::CrossSectionDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 86 of file G4ChipsAntiBaryonInelasticXS.cc.

 {
   outFile << "G4ChipsAntiBaryonInelasticXS provides the inelastic cross\n"
           << "section for anti-baryon 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 G4ChipsAntiBaryonInelasticXS::CrossSectionFormula	(	G4int	targZ,
		G4int	targN,
		G4double	P,
		G4double	lP
	)

private

Definition at line 318 of file G4ChipsAntiBaryonInelasticXS.cc.

 {
   G4double sigma=0.;
   if(tZ==1 && !tN)                        // AntiBar-Prot interaction from G4QuasiElRatios
   {
     G4double ld=lP-3.5;
     G4double ld2=ld*ld;
     G4double ye=G4Exp(lP*1.25);
     G4double yt=G4Exp(lP*0.35);
     G4double El=80./(ye+1.);
     G4double To=(80./yt+.3)/yt;
     sigma=(To-El)+.2443*ld2+31.48;
   }
   else if(tZ==1 && tN==1)
   {
     G4double r=lP-3.7;
     sigma=0.6*r*r+67.+90.*G4Exp(-lP*.666);
   }
   else if(tZ<97 && tN<152)                // General solution
   {
     G4double d=lP-4.2;
     G4double sp=std::sqrt(P);
     G4double a=tN+tZ;                      // A of the target
     G4double sa=std::sqrt(a);
     G4double a2=a*a;
     G4double a3=a2*a;
     G4double a2s=a2*sa;
     G4double c=(170.+3600./a2s)/(1.+65./a2s)+40.*G4Pow::GetInstance()->powA(a,0.712)/(1.+12.2/a)/(1.+34./a2);
     G4double r=(170.+0.01*a3)/(1.+a3/28000.);
     sigma=c+d*d+r/sp;
   }
   else
   {
     G4cerr<<"-Warning-G4QAntiBarNuclearCroSect::CSForm:*Bad A* Z="<<tZ<<", N="<<tN<<G4endl;
     sigma=0.;
   }
   if(sigma<0.) return 0.;
   return sigma;  
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ CrossSectionLin()

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

private

Definition at line 305 of file G4ChipsAntiBaryonInelasticXS.cc.

 {
   G4double lP=G4Log(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 G4ChipsAntiBaryonInelasticXS::CrossSectionLog	(	G4int	targZ,
		G4int	targN,
		G4double	lP
	)

private

Definition at line 312 of file G4ChipsAntiBaryonInelasticXS.cc.

 {
   G4double P=G4Exp(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* G4ChipsAntiBaryonInelasticXS::Default_Name ( )

inlinestatic

Definition at line 52 of file G4ChipsAntiBaryonInelasticXS.hh.

52 {return "ChipsAntiBaryonInelasticXS";}

Here is the call graph for this function:

Here is the caller graph for this function:

◆ EquLinearFit()

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

private

Definition at line 359 of file G4ChipsAntiBaryonInelasticXS.cc.

 {
   if(DX<=0. || N<2)
     {
       G4cerr<<"***G4ChipsAntiBaryonInelasticXS::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 G4ChipsAntiBaryonInelasticXS::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 156 of file G4ChipsAntiBaryonInelasticXS.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,cPDG,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,cPDG,lastZ,lastN,pMom); //calculate & create
       //if(lastCS>0.)                   // It means that the AMBD was initialized
       //{
 
       lastTH = 0; //ThresholdEnergy(tgZ, tgN); // The Threshold Energy which is now the last
         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,cPDG,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 G4ChipsAntiBaryonInelasticXS::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 144 of file G4ChipsAntiBaryonInelasticXS.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:

◆ IsIsoApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 95 of file G4ChipsAntiBaryonInelasticXS.cc.

 {
   /*
   const G4ParticleDefinition* particle = Pt->GetDefinition();
 
   if(particle == G4AntiNeutron::AntiNeutron())
   {
     return true;
   }
   else if(particle == G4AntiProton::AntiProton())
   {
     return true;
   }
   else if(particle == G4AntiLambda::AntiLambda())
   {
     return true;
   }
   else if(particle == G4AntiSigmaPlus::AntiSigmaPlus())
   {
     return true;
   }
   else if(particle == G4AntiSigmaMinus::AntiSigmaMinus())
   {
     return true;
   }
   else if(particle == G4AntiSigmaZero::AntiSigmaZero())
   {
     return true;
   }
   else if(particle == G4AntiXiMinus::AntiXiMinus())
   {
     return true;
   }
   else if(particle == G4AntiXiZero::AntiXiZero())
   {
     return true;
   }
   else if(particle == G4AntiOmegaMinus::AntiOmegaMinus())
   {
     return true;
   }
   */
   return true;
 }