98 std::vector<G4double*>::iterator
pos;
99 for (pos=
CST.begin(); pos<
CST.end(); pos++)
102 for (pos=
PAR.begin(); pos<
PAR.end(); pos++)
105 for (pos=
SST.begin(); pos<
SST.end(); pos++)
108 for (pos=
S1T.begin(); pos<
S1T.end(); pos++)
111 for (pos=
B1T.begin(); pos<
B1T.end(); pos++)
114 for (pos=
S2T.begin(); pos<
S2T.end(); pos++)
117 for (pos=
B2T.begin(); pos<
B2T.end(); pos++)
120 for (pos=
S3T.begin(); pos<
S3T.end(); pos++)
123 for (pos=
B3T.begin(); pos<
B3T.end(); pos++)
126 for (pos=
S4T.begin(); pos<
S4T.end(); pos++)
129 for (pos=
B4T.begin(); pos<
B4T.end(); pos++)
227 if(lastCS<=0. && pEn>
lastTH)
248 colP.push_back(pMom);
331 G4int blast=
static_cast<int>(shift);
350 G4int blast=
static_cast<int>(shift);
391 const G4int n_appel=30;
393 G4double app_el[n_appel]={1.25,3.5,80.,1.,.0557,6.72,5.,74.,3.,3.4,.2,.17,.001,8.,.055,
394 3.64,5.e-5,4000.,1500.,.46,1.2e6,3.5e6,5.e-5,1.e10,8.5e8,
395 1.e10,1.1,3.4e6,6.8e6,0.};
399 if(PDG>-3335 && PDG<-1111)
414 if ( tgZ == 1 && tgN == 0 )
416 for (
G4int ip=0; ip<n_appel; ip++)
lastPAR[ip]=app_el[ip];
441 lastPAR[1]=2.8*asa/(1.+a*(.015+.05/ssa));
443 lastPAR[3]=.013*a2/(1.+a3*(.006+a*.00001));
456 lastPAR[11]=.7/(1.+4.e-12*a16);
457 lastPAR[12]=2.5/a8/(a4+1.e-16*a32);
465 lastPAR[19]=(42.+2.7e-11*a16)/(1.+.14*a);
469 lastPAR[22]=2.4e-3*a8/(1.+2.6e-4*a7);
470 lastPAR[23]=3.5e-36*a32*a8/(1.+5.e-15*a32/
a);
472 lastPAR[24]=1.e5/(a8+2.5e12/a16);
473 lastPAR[25]=8.e7/(a12+1.e-27*a28*a28);
482 lastPAR[32]=1.5e-4/(1.+5.e-6*a12);
490 lastPAR[38]=1.e8*std::exp(.32*asa);
491 lastPAR[39]=20.*std::exp(.45*asa);
493 lastPAR[41]=2.5e5*std::exp(.085*a3);
507 lastPAR[ 9]=4.5*std::pow(a,1.15);
508 lastPAR[10]=.06*std::pow(a,.6);
509 lastPAR[11]=.6*a/(1.+2.e15/a16);
510 lastPAR[12]=.17/(a+9.e5/a3+1.5e33/a32);
511 lastPAR[13]=(.001+7.e-11*
a5)/(1.+4.4e-11*a5);
512 lastPAR[14]=(p1a10*p1a10+2.e-29)/(1.+2.e-22*a12);
514 lastPAR[15]=400./a12+2.e-22*a9;
515 lastPAR[16]=1.e-32*a12/(1.+5.e22/a14);
516 lastPAR[17]=1000./a2+9.5*sa*ssa;
517 lastPAR[18]=4.e-6*a*asa+1.e11/a16;
518 lastPAR[19]=(120./a+.002*
a2)/(1.+2.e14/a16);
526 lastPAR[25]=1.e-5*a2+2.e14/a16;
527 lastPAR[26]=1.2e-11*a2/(1.+1.5e19/a12);
528 lastPAR[27]=.016*asa/(1.+5.e16/a16);
530 lastPAR[28]=.002*a4/(1.+7.e7/std::pow(a-6.83,14));
531 lastPAR[29]=2.e6/a6+7.2/std::pow(a,.11);
532 lastPAR[30]=11.*a3/(1.+7.e23/a16/a8);
536 lastPAR[33]=3.5e-4*a2/(1.+1.e8/a8);
539 lastPAR[36]=1.e-9/a+s4a16*s4a16;
559 lastPAR[51]=1.e15+2.e27/a4/(1.+2.e-18*a16);
589 for(
G4int ip=ini; ip<=fin; ip++)
608 else G4cout<<
"*Warning*G4ChipsAntiBaryonElasticXS::GetPTables: PDG="<<PDG
609 <<
", Z="<<tgZ<<
", N="<<tgN<<
", i="<<ini<<
" > fin="<<fin<<
", LP="<<LP
610 <<
" > ILP="<<ILP<<
" nothing is done!"<<
G4endl;
612 else G4cout<<
"*Warning*G4ChipsAntiBaryonElasticXS::GetPTables: PDG="<<PDG
613 <<
", Z="<<tgZ<<
", N="<<tgN<<
", i="<<ini<<
">= max="<<
nPoints<<
", LP="<<LP
614 <<
" > ILP="<<ILP<<
", lPMax="<<
lPMax<<
" nothing is done!"<<
G4endl;
623 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
624 <<
", while it is defined only for Anti Baryons" <<
G4endl;
625 G4Exception(
"G4ChipsAntiBaryonElasticXS::GetPTables()",
"HAD_CHPS_0000",
639 if(PDG<-3334 || PDG>-1111)
G4cout<<
"*Warning*G4QAntiBaryonElCS::GetExT:PDG="<<PDG<<
G4endl;
648 G4double R2=(1.-std::exp(-E2*E2*E2));
660 q2=-std::log(1.-ran)/
theB1;
666 q2=-std::log(1.-ran);
668 q2=std::pow(q2,third)/
theB2;
674 q2=-std::log(1.-ran)/
theB3;
688 if(a>6.5)E3*=tm2*tm2*tm2;
703 q2=-std::log(1.-ran)/
theB1;
704 if(std::fabs(tss)>1.e-7) q2=(std::sqrt(
theB1*(
theB1+(tss+tss)*q2))-
theB1)/tss;
710 q2=-std::log(1.-ran)/
theB2;
712 if(a<6.5) q2=std::pow(q2,third);
713 else q2=std::pow(q2,fifth);
719 q2=-std::log(1.-ran)/
theB3;
721 if(a>6.5) q2=std::pow(q2,sevth);
727 q2=-std::log(1.-ran)/
theB4;
728 if(a<6.5) q2=lastTM-q2;
732 if(!(q2>=-1.||q2<=1.))
G4cout<<
"*NAN*G4QaBElasticCrossSect::GetExchangeT:-t="<<q2<<
G4endl;
745 if(
lastLP<-4.3)
return 0.;
746 if(PDG<-3334 || PDG>-1111)
752 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
753 <<
", while it is defined only for Anti Baryons" <<
G4endl;
754 G4Exception(
"G4ChipsAntiBaryonElasticXS::GetSlope()",
"HAD_CHPS_0000",
773 if(PDG<-3334 || PDG>-1111)
G4cout<<
"*Warning*G4QAntiBaryElCS::GetTabV:PDG="<<PDG<<
G4endl;
776 G4cout<<
"*Warning*G4QAntiBaryonElCS::GetTabValue:(1-92) NoIsotopesFor Z="<<tgZ<<
G4endl;
791 if ( tgZ == 1 && tgN == 0 )
795 theS1=(lastPAR[7]+lastPAR[8]*dl2*dl2)/(1.+lastPAR[9]/p4/p)+
796 (lastPAR[10]/p2+lastPAR[11]*p)/(p4+lastPAR[12]*sp);
797 theB1=lastPAR[13]*std::pow(p,lastPAR[14])/(1.+lastPAR[15]/p3);
798 theS2=lastPAR[16]+lastPAR[17]/(p4+lastPAR[18]*p);
799 theB2=lastPAR[19]+lastPAR[20]/(p4+lastPAR[21]/
sp);
800 theS3=lastPAR[22]+lastPAR[23]/(p4*p4+lastPAR[24]*p2+lastPAR[25]);
801 theB3=lastPAR[26]+lastPAR[27]/(p4+lastPAR[28]);
805 G4double ye=std::exp(lp*lastPAR[0]);
807 return lastPAR[2]/(ye+lastPAR[3])+lastPAR[4]*dp*dp+lastPAR[5];
834 lastPAR[40]/(1.+lastPAR[41]*std::pow(p,lastPAR[42])));
835 theB4=lastPAR[43]*pa/p2/(1.+pa*lastPAR[44]);
856 return (lastPAR[0]*dlp*dlp+lastPAR[1]+lastPAR[2]/p)/(1.+lastPAR[3]/p);
873 G4double mds=dmt*std::sqrt(pP2+mNuc2)+mNuc2+mt*mt;
874 return dmt*dmt*pP2/mds;
882 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
883 <<
", while it is defined only for p projectiles & Z_target>0" <<
G4endl;
884 G4Exception(
"G4ChipsAntiBaryonElasticXS::GetQ2max()",
"HAD_CHPS_0000",
std::vector< G4double * > B2T
virtual G4bool IsIsoApplicable(const G4DynamicParticle *Pt, G4int Z, G4int A, const G4Element *elm, const G4Material *mat)
std::vector< G4double * > S1T
std::vector< G4double > PIN
std::ostringstream G4ExceptionDescription
static G4AntiOmegaMinus * AntiOmegaMinus()
G4double GetQ2max(G4int pPDG, G4int tgZ, G4int tgN, G4double pP)
std::vector< G4double > colCS
std::vector< G4double * > SST
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
~G4ChipsAntiBaryonElasticXS()
std::vector< G4int > colZ
G4int GetPDGEncoding() const
G4ParticleDefinition * GetDefinition() const
static G4AntiSigmaPlus * AntiSigmaPlus()
std::vector< G4int > colN
G4_DECLARE_XS_FACTORY(G4ChipsAntiBaryonElasticXS)
G4double GetTotalMomentum() const
static G4AntiSigmaMinus * AntiSigmaMinus()
std::vector< G4double * > S4T
G4IonTable * GetIonTable() const
G4GLOB_DLL std::ostream G4cout
static G4AntiProton * AntiProton()
std::vector< G4double * > PAR
G4double GetSlope(G4int tZ, G4int tN, G4int pPDG)
static G4AntiXiMinus * AntiXiMinus()
std::vector< G4double > colTH
static G4Proton * Proton()
std::vector< G4double * > B1T
static G4Neutron * Neutron()
G4double GetTabValues(G4double lp, G4int pPDG, G4int tgZ, G4int tgN)
static const G4double A[nN]
G4ChipsAntiBaryonElasticXS()
std::vector< G4double > colP
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
std::vector< G4double * > B4T
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
std::vector< G4double * > S3T
G4double GetPDGMass() const
static G4ParticleTable * GetParticleTable()
static G4AntiLambda * AntiLambda()
std::vector< G4double * > B3T
static G4AntiSigmaZero * AntiSigmaZero()
static const double gigaelectronvolt
static const double millibarn
G4double GetPTables(G4double lpP, G4double lPm, G4int PDG, G4int tZ, G4int tN)
virtual G4double GetIsoCrossSection(const G4DynamicParticle *, G4int tgZ, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
static G4AntiXiZero * AntiXiZero()
G4double CalculateCrossSection(G4bool CS, G4int F, G4int I, G4int pPDG, G4int Z, G4int N, G4double pP)
std::vector< G4double * > CST
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
static const G4double pos
static G4AntiNeutron * AntiNeutron()
std::vector< G4double * > S2T