106 std::vector<G4double*>::iterator
pos;
107 for (pos=
CST.begin(); pos<
CST.end(); pos++)
110 for (pos=
PAR.begin(); pos<
PAR.end(); pos++)
113 for (pos=
SST.begin(); pos<
SST.end(); pos++)
116 for (pos=
S1T.begin(); pos<
S1T.end(); pos++)
119 for (pos=
B1T.begin(); pos<
B1T.end(); pos++)
122 for (pos=
S2T.begin(); pos<
S2T.end(); pos++)
125 for (pos=
B2T.begin(); pos<
B2T.end(); pos++)
128 for (pos=
S3T.begin(); pos<
S3T.end(); pos++)
131 for (pos=
B3T.begin(); pos<
B3T.end(); pos++)
134 for (pos=
S4T.begin(); pos<
S4T.end(); pos++)
137 for (pos=
B4T.begin(); pos<
B4T.end(); pos++)
145 outFile <<
"G4ChipsHyperonElasticXS provides the elastic cross\n"
146 <<
"section for hyperon nucleus scattering as a function of incident\n"
147 <<
"momentum. The cross section is calculated using M. Kossov's\n"
148 <<
"CHIPS parameterization of cross section data.\n";
206 if(lastCS<=0. && pEn>
lastTH)
227 colP.push_back(pMom);
311 G4int blast=
static_cast<int>(shift);
330 G4int blast=
static_cast<int>(shift);
371 const G4int n_hypel=33;
373 G4double hyp_el[n_hypel]={1.,.002,.12,.0557,3.5,6.72,99.,2.,3.,5.,74.,3.,3.4,.2,.17,
374 .001,8.,.055,3.64,5.e-5,4000.,1500.,.46,1.2e6,3.5e6,5.e-5,
375 1.e10,8.5e8,1.e10,1.1,3.4e6,6.8e6,0.};
379 if(PDG>3000 && PDG<3335)
394 if ( tgZ == 1 && tgN == 0 )
396 for (
G4int ip=0; ip<n_hypel; ip++)
lastPAR[ip]=hyp_el[ip];
421 lastPAR[1]=2.36*asa/(1.+a*.055/ssa);
422 lastPAR[2]=(1.+.00007*a3/ssa)/(1.+.0026*a2);
423 lastPAR[3]=1.76*a/ssa+.00003*a3;
424 lastPAR[4]=(.03+200./a3)/(1.+1.E5/a3/sa);
436 lastPAR[11]=.7/(1.+4.e-12*a16);
437 lastPAR[12]=2.5/a8/(a4+1.e-16*a32);
444 lastPAR[18]=9.e-7/(1.+.035*a5);
445 lastPAR[19]=(42.+2.7e-11*a16)/(1.+.14*a);
449 lastPAR[22]=2.4e-3*a8/(1.+2.6e-4*a7);
450 lastPAR[23]=3.5e-36*a32*a8/(1.+5.e-15*a32/
a);
452 lastPAR[24]=1.e5/(a8+2.5e12/a16);
453 lastPAR[25]=8.e7/(a12+1.e-27*a28*a28);
462 lastPAR[32]=1.5e-4/(1.+5.e-6*a12);
489 lastPAR[11]=.6*a/(1.+2.e15/a16);
490 lastPAR[12]=.17/(a+9.e5/a3+1.5e33/a32);
491 lastPAR[13]=(.001+7.e-11*a5)/(1.+4.4
e-11*a5);
492 lastPAR[14]=(p1a10*p1a10+2.e-29)/(1.+2.
e-22*a12);
494 lastPAR[15]=400./a12+2.e-22*a9;
495 lastPAR[16]=1.e-32*a12/(1.+5.e22/a14);
496 lastPAR[17]=1000./a2+9.5*sa*ssa;
497 lastPAR[18]=4.e-6*a*asa+1.e11/a16;
498 lastPAR[19]=(120./a+.002*a2)/(1.+2.e14/a16);
505 lastPAR[24]=.0011*asa/(1.+3.e34/a32/a4);
506 lastPAR[25]=1.e-5*a2+2.e14/a16;
507 lastPAR[26]=1.2e-11*a2/(1.+1.5e19/a12);
508 lastPAR[27]=.016*asa/(1.+5.e16/a16);
512 lastPAR[30]=11.*a3/(1.+7.e23/a16/a8);
515 lastPAR[32]=(.1+4.4e-5*a2)/(1.+5.e5/a4);
516 lastPAR[33]=3.5e-4*a2/(1.+1.e8/a8);
519 lastPAR[36]=1.e-9/a+s4a16*s4a16;
536 lastPAR[50]=900.*sa/(1.+500./a3);
539 lastPAR[51]=1.e15+2.e27/a4/(1.+2.e-18*a16);
588 else G4cout<<
"*Warning*G4ChipsHyperonElasticXS::GetPTables: PDG="<<PDG
589 <<
", Z="<<tgZ<<
", N="<<tgN<<
", i="<<ini<<
" > fin="<<fin<<
", LP="<<LP
590 <<
" > ILP="<<ILP<<
" nothing is done!"<<
G4endl;
592 else G4cout<<
"*Warning*G4ChipsHyperonElasticXS::GetPTables: PDG="<<PDG
593 <<
", Z="<<tgZ<<
", N="<<tgN<<
", i="<<ini<<
">= max="<<
nPoints<<
", LP="<<LP
594 <<
" > ILP="<<ILP<<
", lPMax="<<
lPMax<<
" nothing is done!"<<
G4endl;
601 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
602 <<
", while it is defined only for Hyperons" <<
G4endl;
603 G4Exception(
"G4ChipsHyperonElasticXS::GetPTables()",
"HAD_CHPS_0000",
617 if(PDG<3000 || PDG>3334)
G4cout<<
"*Warning*G4QHyElCS::GET:PDG="<<PDG<<
G4endl;
666 if(a>6.5)E3*=tm2*tm2*tm2;
682 if(std::fabs(tss)>1.
e-7) q2=(std::sqrt(
theB1*(
theB1+(tss+tss)*q2))-
theB1)/tss;
706 if(a<6.5) q2=lastTM-q2;
710 if(!(q2>=-1.||q2<=1.))
G4cout<<
"*NAN*G4QHyElasticCrossSect::GetExchangeT:-t="<<q2<<
G4endl;
723 if(
lastLP<-4.3)
return 0.;
725 if(PDG<3000 || PDG>3334)
731 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
732 <<
", while it is defined only for Hyperons" <<
G4endl;
733 G4Exception(
"G4ChipsHyperonElasticXS::GetSlope()",
"HAD_CHPS_0000",
753 if(PDG<3000 || PDG>3334)
G4cout<<
"*Warning*G4QHypElCS::GTV:P="<<PDG<<
G4endl;
757 const G4bool isHeavyElementAllowed =
true;
758 if(tgZ<0 || ( !isHeavyElementAllowed && tgZ>92))
760 G4cout<<
"*Warning*G4QHyperonElastCS::GetTabValue:(1-92) NoIsotopesFor Z="<<tgZ<<
G4endl;
775 if ( tgZ == 1 && tgN == 0 )
779 theS1=(lastPAR[10]+lastPAR[11]*dl2*dl2)/(1.+lastPAR[12]/p4/p)+
780 (lastPAR[13]/p2+lastPAR[14]*
p)/(p4+lastPAR[15]*sp);
782 theS2=lastPAR[19]+lastPAR[20]/(p4+lastPAR[21]*
p);
783 theB2=lastPAR[22]+lastPAR[23]/(p4+lastPAR[24]/
sp);
784 theS3=lastPAR[25]+lastPAR[26]/(p4*p4+lastPAR[27]*p2+lastPAR[28]);
785 theB3=lastPAR[29]+lastPAR[30]/(p4+lastPAR[31]);
790 return lastPAR[0]/(lastPAR[1]+p2*(lastPAR[2]+p2))+(lastPAR[3]*dp*dp+lastPAR[5]+
791 lastPAR[6]/p2)/(1.+lastPAR[7]/sp+lastPAR[8]/p4);
819 theB4=lastPAR[43]*pa/p2/(1.+pa*lastPAR[44]);
840 return (lastPAR[0]*dlp*dlp+lastPAR[1])/(1.+lastPAR[2]/
p)+lastPAR[3]/(p3+lastPAR[4]);
850 static const G4double mLa2= mLamb*mLamb;
857 G4double mds=dmt*std::sqrt(pP2+mLa2)+mLa2+mt*mt;
858 return dmt*dmt*pP2/mds;
866 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
867 <<
", while it is defined only for p projectiles & Z_target>0" <<
G4endl;
868 G4Exception(
"G4ChipsHyperonElasticXS::GetQ2max()",
"HAD_CHPS_0000",
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
static G4Lambda * Lambda()
std::ostringstream G4ExceptionDescription
std::vector< ExP01TrackerHit * > a
static const G4double pos
std::vector< G4int > colN
static G4ParticleTable * GetParticleTable()
std::vector< G4double * > SST
std::vector< G4double * > S4T
std::vector< G4double * > S1T
std::vector< G4double * > CST
std::vector< G4int > colZ
G4IonTable * GetIonTable() const
std::vector< G4double * > B4T
G4double GetPDGMass() const
G4double G4Log(G4double x)
std::vector< G4double * > B1T
static constexpr double gigaelectronvolt
std::vector< G4double * > S3T
std::vector< G4double * > S2T
G4double GetPTables(G4double lpP, G4double lPm, G4int PDG, G4int tZ, G4int tN)
G4ChipsHyperonElasticXS()
static constexpr double millibarn
static G4Pow * GetInstance()
G4ParticleDefinition * GetDefinition() const
std::vector< G4double > PIN
virtual void CrossSectionDescription(std::ostream &) const
std::vector< G4double > colCS
G4double powA(G4double A, G4double y) const
G4double GetTabValues(G4double lp, G4int pPDG, G4int tgZ, G4int tgN)
std::vector< G4double * > B2T
~G4ChipsHyperonElasticXS()
double A(double temperature)
G4double GetQ2max(G4int pPDG, G4int tgZ, G4int tgN, G4double pP)
#define G4_DECLARE_XS_FACTORY(cross_section)
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)
std::vector< G4double * > PAR
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
G4double CalculateCrossSection(G4bool CS, G4int F, G4int I, G4int pPDG, G4int Z, G4int N, G4double pP)
std::vector< G4double > colTH
std::vector< G4double * > B3T
G4int GetPDGEncoding() const
G4GLOB_DLL std::ostream G4cout
G4double GetSlope(G4int tZ, G4int tN, G4int pPDG)
G4double GetTotalMomentum() const
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
static constexpr double GeV
std::vector< G4double > colP
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)