69 LEN =
new std::vector<G4double*>;
70 HEN =
new std::vector<G4double*>;
76 for(
G4int i=0; i<lens; ++i)
delete[] (*
LEN)[i];
79 for(
G4int i=0; i<hens; ++i)
delete[] (*
HEN)[i];
145 static G4ThreadLocal std::vector <G4int> *colN_G4MT_TLS_ = 0 ;
if (!colN_G4MT_TLS_) colN_G4MT_TLS_ =
new std::vector <G4int> ; std::vector <G4int> &colN = *colN_G4MT_TLS_;
146 static G4ThreadLocal std::vector <G4int> *colZ_G4MT_TLS_ = 0 ;
if (!colZ_G4MT_TLS_) colZ_G4MT_TLS_ =
new std::vector <G4int> ; std::vector <G4int> &colZ = *colZ_G4MT_TLS_;
147 static G4ThreadLocal std::vector <G4double> *colP_G4MT_TLS_ = 0 ;
if (!colP_G4MT_TLS_) colP_G4MT_TLS_ =
new std::vector <G4double> ; std::vector <G4double> &colP = *colP_G4MT_TLS_;
148 static G4ThreadLocal std::vector <G4double> *colTH_G4MT_TLS_ = 0 ;
if (!colTH_G4MT_TLS_) colTH_G4MT_TLS_ =
new std::vector <G4double> ; std::vector <G4double> &colTH = *colTH_G4MT_TLS_;
149 static G4ThreadLocal std::vector <G4double> *colCS_G4MT_TLS_ = 0 ;
if (!colCS_G4MT_TLS_) colCS_G4MT_TLS_ =
new std::vector <G4double> ; std::vector <G4double> &colCS = *colCS_G4MT_TLS_;
163 if(colN[i]==tgN && colZ[i]==tgZ)
176 if(lastCS<=0. && pMom>
lastTH)
195 colP.push_back(pMom);
224 static const G4double THmiG=THmin*.001;
228 static const G4double Pmin=THmin+(nL-1)*dP;
231 static const G4double milP=std::log(Pmin);
232 static const G4double malP=std::log(Pmax);
233 static const G4double dlP=(malP-milP)/(nH-1);
234 static const G4double milPG=std::log(.001*Pmin);
243 if(sync<=I)
G4cerr<<
"*!*G4QPiMinusNuclCS::CalcCrosSect:Sync="<<sync<<
"<="<<I<<
G4endl;
269 G4cerr<<
"***G4QPiMinusNuclCS::CalcCrossSect: Sinc="<<sync<<
"#"<<I<<
", Z=" <<targZ
270 <<
", N="<<targN<<
", F="<<F<<
G4endl;
278 if (Momentum<
lastTH)
return 0.;
279 else if (Momentum<Pmin)
283 else if (Momentum<Pmax)
293 if(sigma<0.)
return 0.;
323 sigma=(To-El)+.2443*ld2+31.48;
325 else if(tZ==1 && tN==1)
328 sigma=0.6*r*r+67.+90.*std::exp(-lP*.666);
330 else if(tZ<97 && tN<152)
339 G4double c=(170.+3600./a2s)/(1.+65./a2s)+40.*std::pow(a,0.712)/(1.+12.2/
a)/(1.+34./a2);
345 G4cerr<<
"-Warning-G4QAntiBarNuclearCroSect::CSForm:*Bad A* Z="<<tZ<<
", N="<<tN<<
G4endl;
348 if(sigma<0.)
return 0.;
356 G4cerr<<
"***G4ChipsAntiBaryonInelasticXS::EquLinearFit: DX="<<DX<<
", N="<<N<<
G4endl;
362 G4int j=
static_cast<int>(d);
static G4AntiOmegaMinus * AntiOmegaMinus()
G4double CalculateCrossSection(G4int F, G4int I, G4int PDG, G4int Z, G4int N, G4double Momentum)
std::vector< G4double * > * HEN
G4int GetPDGEncoding() const
G4ParticleDefinition * GetDefinition() const
G4double CrossSectionFormula(G4int targZ, G4int targN, G4double P, G4double lP)
static G4AntiSigmaPlus * AntiSigmaPlus()
std::vector< G4double * > * LEN
G4double GetTotalMomentum() const
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)
static G4AntiSigmaMinus * AntiSigmaMinus()
static G4AntiProton * AntiProton()
G4double EquLinearFit(G4double X, G4int N, G4double X0, G4double DX, G4double *Y)
static G4AntiXiMinus * AntiXiMinus()
G4ChipsAntiBaryonInelasticXS()
G4_DECLARE_XS_FACTORY(G4ChipsAntiBaryonInelasticXS)
static const G4double A[nN]
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
static G4AntiLambda * AntiLambda()
G4double CrossSectionLin(G4int targZ, G4int targN, G4double P)
static G4AntiSigmaZero * AntiSigmaZero()
static const double millibarn
static G4AntiXiZero * AntiXiZero()
G4double CrossSectionLog(G4int targZ, G4int targN, G4double lP)
static const G4double THmin
static G4AntiNeutron * AntiNeutron()
~G4ChipsAntiBaryonInelasticXS()
G4GLOB_DLL std::ostream G4cerr