63 LEN =
new std::vector<G4double*>;
64 HEN =
new std::vector<G4double*>;
70 for(
G4int i=0; i<lens; ++i)
delete[] (*
LEN)[i];
73 for(
G4int i=0; i<hens; ++i)
delete[] (*
HEN)[i];
125 if(lastCS<=0. && pMom>
lastTH)
144 colP.push_back(pMom);
173 static const G4double THmiG=THmin*.001;
177 static const G4double Pmin=THmin+(nL-1)*dP;
180 static const G4double milP=std::log(Pmin);
181 static const G4double malP=std::log(Pmax);
182 static const G4double dlP=(malP-milP)/(nH-1);
183 static const G4double milPG=std::log(.001*Pmin);
192 if(sync<=I)
G4cerr<<
"*!*G4ChipsPiMinusNuclCS::CalcCrosSect:Sync="<<sync<<
"<="<<I<<
G4endl;
218 G4cerr<<
"***G4ChipsPiMinusNuclCS::CalcCrossSect: Sinc="<<sync<<
"#"<<I<<
", Z=" <<targZ
219 <<
", N="<<targN<<
", F="<<F<<
G4endl;
227 if (Momentum<
lastTH)
return 0.;
228 else if (Momentum<Pmin)
232 else if (Momentum<Pmax)
242 if(sigma<0.)
return 0.;
279 sigma=(To-El)+.4/md+.01/hd;
282 else if(tZ==1 && tN==1)
288 sigma=(.55*d*d+38.+23./std::sqrt(P))/(1.+.3/p2/p2)+18./(f*f+.1089)+.02/(gg*gg+.0025);
290 else if(tZ<97 && tN<152)
300 G4double c=41.*std::exp(al*.68)*(1.+44./
a2)/(1.+8./a)/(1.+200./a2/
a2);
305 sigma=(c+d*d)/(1.+.17/p4)+f/(u*u+h*h);
309 G4cerr<<
"-Warning-G4ChipsPiMinusNuclearCroSect::CSForm:*Bad A* Z="<<tZ<<
", N="<<tN<<
G4endl;
312 if(sigma<0.)
return 0.;
320 G4cerr<<
"***G4ChipsPionMinusInelasticXS::EquLinearFit: DX="<<DX<<
", N="<<N<<
G4endl;
326 G4int jj=
static_cast<int>(d);
328 else if(jj>N2) jj=N2;
G4_DECLARE_XS_FACTORY(G4ChipsPionMinusInelasticXS)
virtual G4bool IsIsoApplicable(const G4DynamicParticle *Pt, G4int Z, G4int A, const G4Element *elm, const G4Material *mat)
G4double EquLinearFit(G4double X, G4int N, G4double X0, G4double DX, G4double *Y)
virtual G4double GetIsoCrossSection(const G4DynamicParticle *, G4int tgZ, G4int A, const G4Isotope *iso=0, const G4Element *elm=0, const G4Material *mat=0)
G4double CalculateCrossSection(G4int F, G4int I, G4int PDG, G4int Z, G4int N, G4double Momentum)
~G4ChipsPionMinusInelasticXS()
G4ChipsPionMinusInelasticXS()
std::vector< G4double * > * LEN
G4double GetTotalMomentum() const
std::vector< G4double > colP
std::vector< G4double * > * HEN
G4double CrossSectionLin(G4int targZ, G4int targN, G4double P)
std::vector< G4double > colCS
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
static const G4double A[nN]
G4double CrossSectionLog(G4int targZ, G4int targN, G4double lP)
G4double CrossSectionFormula(G4int targZ, G4int targN, G4double P, G4double lP)
static const double millibarn
std::vector< G4double > colTH
std::vector< G4int > colN
std::vector< G4int > colZ
static const G4double THmin
G4GLOB_DLL std::ostream G4cerr