65 const G4double dlP=(malP-milP)/(nH-1);
66 const G4double milPG=std::log(.001*Pmin);
92 LEN =
new std::vector<G4double*>;
93 HEN =
new std::vector<G4double*>;
99 for(
G4int i=0; i<lens; ++i)
delete[] (*
LEN)[i];
103 for(
G4int i=0; i<hens; ++i)
delete[] (*
HEN)[i];
159 if(lastCS<=0. && pMom>
lastTH)
179 colP.push_back(pMom);
216 if(sync<=I)
G4cerr<<
"*!*G4ChipsKPlusNuclCS::CalcCrosSect:Sync="<<sync<<
"<="<<I<<
G4endl;
242 G4cerr<<
"***G4ChipsKPlusNuclCS::CalcCrossSect: Sinc="<<sync<<
"#"<<I<<
", Z=" <<targZ
243 <<
", N="<<targN<<
", F="<<F<<
G4endl;
252 if (Momentum<
lastTH)
return 0.;
253 else if (Momentum<Pmin)
255 if(A<=1. && Momentum < 600.) sigma=0.;
258 else if (Momentum<Pmax)
268 if(sigma<0.)
return 0.;
276 if(tZ<.99 || tN<0.)
return 0.;
279 if(tZ==1 && tN==0) tM=prM;
280 else dE=tZ/(1.+std::pow(tA,third));
283 return std::sqrt(T*(tpM+T));
313 G4double El=(.0557*ld2+2.23)/(1.-.7/sp+.1/p4);
314 G4double To=(.3*ld2+19.5)/(1.+.46/sp+1.6/p4);
317 else if(tZ<97 && tN<152)
343 c=52.*std::exp(al*.6)*(1.+95./
a2)/(1.+9./a)/(1.+46./
a2);
344 ss=(40.+.14*
a)/(1.+12./a);
346 t=.185*y/(1.+.00012*y);
347 u=(1.+80./asa)/(1.+200./asa);
348 v=(1.+3.E-6*a4*(1.+6.E-7*a3+4.E10/a12))/a3/20000.;
355 sigma=(f*d*d+c)/(1.+r/std::sqrt(P)+1./p4)+(rD+rR)/(1+v/p4/p4);
359 G4cerr<<
"-Warning-G4ChipsKaonPlusNuclearCroSect::CSForm:Bad A, Z="<<tZ<<
", N="<<tN<<
G4endl;
362 if(sigma<0.)
return 0.;
370 G4cerr<<
"***G4ChipsKaonPlusInelasticXS::EquLinearFit: DX="<<DX<<
", N="<<N<<
G4endl;
376 G4int jj=
static_cast<int>(d);
378 else if(jj>N2) jj=N2;
std::vector< G4double > colP
std::vector< G4int > colZ
std::vector< G4int > colN
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
G4double CrossSectionLog(G4int targZ, G4int targN, G4double lP)
G4ChipsKaonPlusInelasticXS()
~G4ChipsKaonPlusInelasticXS()
G4double ThresholdMomentum(G4int targZ, G4int targN)
std::vector< G4double > colCS
virtual G4double GetIsoCrossSection(const G4DynamicParticle *, G4int tgZ, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
std::vector< G4double * > * HEN
#define G4MUTEX_INITIALIZER
std::vector< G4double * > * LEN
G4double GetTotalMomentum() const
G4double CalculateCrossSection(G4int F, G4int I, G4int PDG, G4int Z, G4int N, G4double Momentum)
G4_DECLARE_XS_FACTORY(G4ChipsKaonPlusInelasticXS)
G4double CrossSectionLin(G4int targZ, G4int targN, G4double P)
std::vector< G4double > colTH
static G4Proton * Proton()
static G4PionPlus * PionPlus()
G4double EquLinearFit(G4double X, G4int N, G4double X0, G4double DX, G4double *Y)
static const G4double A[nN]
G4double GetPDGMass() const
G4double CrossSectionFormula(G4int targZ, G4int targN, G4double P, G4double lP)
static const double millibarn
static G4KaonPlus * KaonPlus()
static const G4double THmin
virtual G4bool IsIsoApplicable(const G4DynamicParticle *Pt, G4int Z, G4int A, const G4Element *elm, const G4Material *mat)
G4GLOB_DLL std::ostream G4cerr