67 const G4double kmp_el[
n_kmpel]={5.2,.0557,3.5,2.23,.7,.075,.004,.39,.000156,.15,1.,.0156,5.,
68 74.,3.,3.4,.2,.17,.001,8.,.055,3.64,5.e-5,4000.,1500.,.46,
69 1.2e6,3.5e6,5.e-5,1.e10,8.5e8,1.e10,1.1,3.4e6,6.8e6,0.};
121 std::vector<G4double*>::iterator
pos;
160 outFile <<
"G4ChipsKaonMinusElasticXS provides the elastic cross\n"
161 <<
"section for K- nucleus scattering as a function of incident\n"
162 <<
"momentum. The cross section is calculated using M. Kossov's\n"
163 <<
"CHIPS parameterization of cross section data.\n";
220 if(lastCS<=0. && pEn>
lastTH)
241 colP.push_back(pMom);
324 G4int blast=
static_cast<int>(shift);
343 G4int blast=
static_cast<int>(shift);
398 if ( tgZ == 1 && tgN == 0 )
424 lastPAR[0]=.06*asa/(1.+a*(.01+.1/ssa));
425 lastPAR[1]=.75*asa/(1.+.009*a);
426 lastPAR[2]=.1*a2*ssa/(1.+.0015*a2/ssa);
440 lastPAR[11]=.7/(1.+4.e-12*a16);
441 lastPAR[12]=2.5/a8/(a4+1.e-16*a32);
448 lastPAR[18]=9.e-7/(1.+.035*a5);
449 lastPAR[19]=(42.+2.7e-11*a16)/(1.+.14*a);
453 lastPAR[22]=2.4e-3*a8/(1.+2.6e-4*a7);
454 lastPAR[23]=3.5e-36*a32*a8/(1.+5.e-15*a32/a);
456 lastPAR[24]=1.e5/(a8+2.5e12/a16);
457 lastPAR[25]=8.e7/(a12+1.e-27*a28*a28);
466 lastPAR[32]=1.5e-4/(1.+5.e-6*a12);
474 lastPAR[38]=1.e8*std::exp(.32*asa);
475 lastPAR[39]=20.*std::exp(.45*asa);
477 lastPAR[41]=2.5e5*std::exp(.085*a3);
491 lastPAR[ 9]=4.5*std::pow(a,1.15);
492 lastPAR[10]=.06*std::pow(a,.6);
493 lastPAR[11]=.6*a/(1.+2.e15/a16);
494 lastPAR[12]=.17/(a+9.e5/a3+1.5e33/a32);
495 lastPAR[13]=(.001+7.e-11*a5)/(1.+4.4e-11*a5);
496 lastPAR[14]=(p1a10*p1a10+2.e-29)/(1.+2.e-22*a12);
498 lastPAR[15]=400./a12+2.e-22*a9;
499 lastPAR[16]=1.e-32*a12/(1.+5.e22/a14);
500 lastPAR[17]=1000./a2+9.5*sa*ssa;
501 lastPAR[18]=4.e-6*a*asa+1.e11/a16;
502 lastPAR[19]=(120./a+.002*a2)/(1.+2.e14/a16);
509 lastPAR[24]=.0011*asa/(1.+3.e34/a32/a4);
510 lastPAR[25]=1.e-5*a2+2.e14/a16;
511 lastPAR[26]=1.2e-11*a2/(1.+1.5e19/a12);
512 lastPAR[27]=.016*asa/(1.+5.e16/a16);
514 lastPAR[28]=.002*a4/(1.+7.e7/std::pow(a-6.83,14));
515 lastPAR[29]=2.e6/a6+7.2/std::pow(a,.11);
516 lastPAR[30]=11.*a3/(1.+7.e23/a16/a8);
519 lastPAR[32]=(.1+4.4e-5*a2)/(1.+5.e5/a4);
520 lastPAR[33]=3.5e-4*a2/(1.+1.e8/a8);
523 lastPAR[36]=1.e-9/a+s4a16*s4a16;
540 lastPAR[50]=900.*sa/(1.+500./a3);
543 lastPAR[51]=1.e15+2.e27/a4/(1.+2.e-18*a16);
573 for(
G4int ip=ini; ip<=fin; ip++)
592 else G4cout<<
"*Warning*G4ChipsKaonMinusElasticXS::GetPTables: PDG="<<PDG
593 <<
", Z="<<tgZ<<
", N="<<tgN<<
", i="<<ini<<
" > fin="<<fin<<
", LP="<<LP
594 <<
" > ILP="<<ILP<<
" nothing is done!"<<
G4endl;
596 else G4cout<<
"*Warning*G4ChipsKaonMinusElasticXS::GetPTables: PDG="<<PDG
597 <<
", Z="<<tgZ<<
", N="<<tgN<<
", i="<<ini<<
">= max="<<
nPoints<<
", LP="<<LP
598 <<
" > ILP="<<ILP<<
", lPMax="<<
lPMax<<
" nothing is done!"<<
G4endl;
607 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
608 <<
", while it is defined only for PDG=-321 (K-) " <<
G4endl;
609 G4Exception(
"G4ChipsKaonMinusElasticXS::GetPTables()",
"HAD_CHPS_0000",
618 if(PDG==310 || PDG==130) PDG=-321;
619 if(PDG!=-321)
G4cout<<
"*Warning*G4ChipsKaonMinusElasticXS::GetET:PDG="<<PDG<<
G4endl;
628 G4double R2=(1.-std::exp(-E2*E2*E2));
640 q2=-std::log(1.-ran)/
theB1;
646 q2=-std::log(1.-ran);
654 q2=-std::log(1.-ran)/
theB3;
668 if(a>6.5)E3*=tm2*tm2*tm2;
683 q2=-std::log(1.-ran)/
theB1;
684 if(std::fabs(tss)>1.e-7) q2=(std::sqrt(
theB1*(
theB1+(tss+tss)*q2))-
theB1)/tss;
690 q2=-std::log(1.-ran)/
theB2;
692 if(a<6.5) q2=std::pow(q2,
third);
693 else q2=std::pow(q2,
fifth);
699 q2=-std::log(1.-ran)/
theB3;
701 if(a>6.5) q2=std::pow(q2,
sevth);
707 q2=-std::log(1.-ran)/
theB4;
712 if(!(q2>=-1.||q2<=1.))
G4cout<<
"*NAN*G4QKaonMinusElasticCS::GetExchT: -t="<<q2<<
G4endl;
724 if(
lastLP<-4.3)
return 0.;
731 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
732 <<
", while it is defined only for PDG=-321 (K-)" <<
G4endl;
749 if(PDG!=-321)
G4cout<<
"*Warning*G4ChipsKaonMinusElasticXS::GetTV:PDG="<<PDG<<
G4endl;
752 const G4bool isHeavyElementAllowed =
true;
753 if(tgZ<0 || ( !isHeavyElementAllowed && tgZ>92))
755 G4cout<<
"*Warning*G4QKaonMinusElasticCS::GetTabV:(1-92)NoIsotopes for Z="<<tgZ<<
G4endl;
771 if ( tgZ == 1 && tgN == 0 )
852 return dmt*dmt*pP2/mds;
857 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
858 <<
", while it is defined only for p projectiles & Z_target>0" <<
G4endl;
859 G4Exception(
"G4ChipsKaonMinusElasticXS::GetQ2max()",
"HAD_CHPS_0000",
G4_DECLARE_XS_FACTORY(G4ChipsKaonMinusElasticXS)
static const G4double pos
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
std::ostringstream G4ExceptionDescription
static constexpr double gigaelectronvolt
static constexpr double millibarn
static constexpr double GeV
#define G4MUTEX_INITIALIZER
G4GLOB_DLL std::ostream G4cout
std::vector< G4double > colCS
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
virtual void CrossSectionDescription(std::ostream &) const
std::vector< G4double * > S4T
virtual G4double GetIsoCrossSection(const G4DynamicParticle *, G4int tgZ, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
G4double CalculateCrossSection(G4bool CS, G4int F, G4int I, G4int pPDG, G4int Z, G4int N, G4double pP)
std::vector< G4double * > S3T
std::vector< G4double > colP
G4double GetQ2max(G4int pPDG, G4int tgZ, G4int tgN, G4double pP)
std::vector< G4double * > CST
std::vector< G4double * > B3T
std::vector< G4double * > B2T
G4double GetTabValues(G4double lp, G4int pPDG, G4int tgZ, G4int tgN)
std::vector< G4double * > PAR
~G4ChipsKaonMinusElasticXS()
std::vector< G4double * > B1T
G4double GetSlope(G4int tZ, G4int tN, G4int pPDG)
std::vector< G4int > colZ
G4ChipsKaonMinusElasticXS()
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
G4double GetPTables(G4double lpP, G4double lPm, G4int PDG, G4int tZ, G4int tN)
virtual G4bool IsIsoApplicable(const G4DynamicParticle *Pt, G4int Z, G4int A, const G4Element *elm, const G4Material *mat)
std::vector< G4double > colTH
std::vector< G4double > PIN
std::vector< G4double * > S2T
std::vector< G4int > colN
std::vector< G4double * > SST
std::vector< G4double * > B4T
std::vector< G4double * > S1T
G4double GetTotalMomentum() const
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
static G4KaonMinus * KaonMinus()
G4double GetPDGMass() const
G4IonTable * GetIonTable() const
static G4ParticleTable * GetParticleTable()
const G4double kmp_el[n_kmpel]