69 LEN =
new std::vector<G4double*>;
70 HEN =
new std::vector<G4double*>;
75 G4int lens=LEN->size();
76 for(
G4int i=0; i<lens; ++i)
delete[] (*LEN)[i];
78 G4int hens=HEN->size();
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_;
153 if(tgN!=lastN || tgZ!=lastZ)
161 if(lastI)
for(
G4int i=0; i<lastI; i++)
163 if(colN[i]==tgN && colZ[i]==tgZ)
175 lastCS=CalculateCrossSection(-1,j,cPDG,lastZ,lastN,pMom);
176 if(lastCS<=0. && pMom>lastTH)
188 lastCS=CalculateCrossSection(0,j,cPDG,lastZ,lastN,pMom);
195 colP.push_back(pMom);
196 colTH.push_back(lastTH);
197 colCS.push_back(lastCS);
207 else if(pMom<=lastTH)
213 lastCS=CalculateCrossSection(1,j,cPDG,lastZ,lastN,pMom);
224 static const G4double THmiG=THmin*.001;
227 static const G4int nL=105;
228 static const G4double Pmin=THmin+(nL-1)*dP;
230 static const G4int nH=224;
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);
242 G4int sync=LEN->size();
243 if(sync<=I)
G4cerr<<
"*!*G4QPiMinusNuclCS::CalcCrosSect:Sync="<<sync<<
"<="<<I<<
G4endl;
253 for(
G4int k=0; k<nL; k++)
255 lastLEN[k] = CrossSectionLin(targZ, targN, P);
261 lastHEN[
n] = CrossSectionLog(targZ, targN, lP);
266 G4int sync=LEN->size();
269 G4cerr<<
"***G4QPiMinusNuclCS::CalcCrossSect: Sinc="<<sync<<
"#"<<I<<
", Z=" <<targZ
270 <<
", N="<<targN<<
", F="<<F<<
G4endl;
273 LEN->push_back(lastLEN);
274 HEN->push_back(lastHEN);
278 if (Momentum<lastTH)
return 0.;
279 else if (Momentum<Pmin)
281 sigma=EquLinearFit(Momentum,nL,THmin,dP,lastLEN);
283 else if (Momentum<Pmax)
286 sigma=EquLinearFit(lP,nH,milP,dlP,lastHEN);
291 sigma=CrossSectionFormula(targZ, targN, P, std::log(P));
293 if(sigma<0.)
return 0.;
301 return CrossSectionFormula(tZ, tN, P, lP);
308 return CrossSectionFormula(tZ, tN, P, lP);
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);
340 G4double r=(170.+0.01*a3)/(1.+a3/28000.);
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()
G4int GetPDGEncoding() const
G4ParticleDefinition * GetDefinition() const
static G4AntiSigmaPlus * AntiSigmaPlus()
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()
static G4AntiXiMinus * AntiXiMinus()
G4ChipsAntiBaryonInelasticXS()
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
#define G4_DECLARE_XS_FACTORY(cross_section)
static G4AntiLambda * AntiLambda()
static G4AntiSigmaZero * AntiSigmaZero()
static G4AntiXiZero * AntiXiZero()
static G4AntiNeutron * AntiNeutron()
~G4ChipsAntiBaryonInelasticXS()
G4GLOB_DLL std::ostream G4cerr