49const G4int G4QHyperonElasticCrossSection::nPoints=128;
50const G4int G4QHyperonElasticCrossSection::nLast=nPoints-1;
51G4double G4QHyperonElasticCrossSection::lPMin=-8.;
52G4double G4QHyperonElasticCrossSection::lPMax= 8.;
53G4double G4QHyperonElasticCrossSection::dlnP=(lPMax-lPMin)/nLast;
54G4bool G4QHyperonElasticCrossSection::onlyCS=
true;
55G4double G4QHyperonElasticCrossSection::lastSIG=0.;
56G4double G4QHyperonElasticCrossSection::lastLP=-10.;
57G4double G4QHyperonElasticCrossSection::lastTM=0.;
58G4double G4QHyperonElasticCrossSection::theSS=0.;
59G4double G4QHyperonElasticCrossSection::theS1=0.;
60G4double G4QHyperonElasticCrossSection::theB1=0.;
61G4double G4QHyperonElasticCrossSection::theS2=0.;
62G4double G4QHyperonElasticCrossSection::theB2=0.;
63G4double G4QHyperonElasticCrossSection::theS3=0.;
64G4double G4QHyperonElasticCrossSection::theB3=0.;
65G4double G4QHyperonElasticCrossSection::theS4=0.;
66G4double G4QHyperonElasticCrossSection::theB4=0.;
67G4int G4QHyperonElasticCrossSection::lastTZ=0;
68G4int G4QHyperonElasticCrossSection::lastTN=0;
69G4double G4QHyperonElasticCrossSection::lastPIN=0.;
70G4double* G4QHyperonElasticCrossSection::lastCST=0;
71G4double* G4QHyperonElasticCrossSection::lastPAR=0;
72G4double* G4QHyperonElasticCrossSection::lastSST=0;
73G4double* G4QHyperonElasticCrossSection::lastS1T=0;
74G4double* G4QHyperonElasticCrossSection::lastB1T=0;
75G4double* G4QHyperonElasticCrossSection::lastS2T=0;
76G4double* G4QHyperonElasticCrossSection::lastB2T=0;
77G4double* G4QHyperonElasticCrossSection::lastS3T=0;
78G4double* G4QHyperonElasticCrossSection::lastB3T=0;
79G4double* G4QHyperonElasticCrossSection::lastS4T=0;
80G4double* G4QHyperonElasticCrossSection::lastB4T=0;
81G4int G4QHyperonElasticCrossSection::lastN=0;
82G4int G4QHyperonElasticCrossSection::lastZ=0;
83G4double G4QHyperonElasticCrossSection::lastP=0.;
84G4double G4QHyperonElasticCrossSection::lastTH=0.;
85G4double G4QHyperonElasticCrossSection::lastCS=0.;
86G4int G4QHyperonElasticCrossSection::lastI=0;
88std::vector<G4double*> G4QHyperonElasticCrossSection::PAR;
89std::vector<G4double*> G4QHyperonElasticCrossSection::CST;
90std::vector<G4double*> G4QHyperonElasticCrossSection::SST;
91std::vector<G4double*> G4QHyperonElasticCrossSection::S1T;
92std::vector<G4double*> G4QHyperonElasticCrossSection::B1T;
93std::vector<G4double*> G4QHyperonElasticCrossSection::S2T;
94std::vector<G4double*> G4QHyperonElasticCrossSection::B2T;
95std::vector<G4double*> G4QHyperonElasticCrossSection::S3T;
96std::vector<G4double*> G4QHyperonElasticCrossSection::B3T;
97std::vector<G4double*> G4QHyperonElasticCrossSection::S4T;
98std::vector<G4double*> G4QHyperonElasticCrossSection::B4T;
106 std::vector<G4double*>::iterator pos;
107 for (pos=CST.begin(); pos<CST.end(); pos++)
110 for (pos=PAR.begin(); pos<PAR.end(); pos++)
113 for (pos=SST.begin(); pos<SST.end(); pos++)
116 for (pos=S1T.begin(); pos<S1T.end(); pos++)
119 for (pos=B1T.begin(); pos<B1T.end(); pos++)
122 for (pos=S2T.begin(); pos<S2T.end(); pos++)
125 for (pos=B2T.begin(); pos<B2T.end(); pos++)
128 for (pos=S3T.begin(); pos<S3T.end(); pos++)
131 for (pos=B3T.begin(); pos<B3T.end(); pos++)
134 for (pos=S4T.begin(); pos<S4T.end(); pos++)
137 for (pos=B4T.begin(); pos<B4T.end(); pos++)
146 return &theCrossSection;
154 static std::vector <G4int> colN;
155 static std::vector <G4int> colZ;
156 static std::vector <G4double> colP;
157 static std::vector <G4double> colTH;
158 static std::vector <G4double> colCS;
163 G4cout<<
"G4QHyperElCS::GetCS:>>> f="<<fCS<<
", p="<<pMom<<
", Z="<<tgZ<<
"("<<lastZ<<
") ,N="
164 <<tgN<<
"("<<lastN<<
"), T="<<pEn<<
"("<<lastTH<<
")"<<
",Sz="<<colN.size()<<
G4endl;
167 if(pPDG==3222 || pPDG<3000 || pPDG>3334)
169 G4cout<<
"*Warning*G4QHyperonElaCS::GetCS:**> Found pPDG="<<pPDG<<
" =--=> CS=0"<<
G4endl;
178 if(lastI)
for(
G4int i=0; i<lastI; i++)
180 if(colN[i]==tgN && colZ[i]==tgZ)
185 G4cout<<
"G4QHyElCS::GetCS:*Found* P="<<pMom<<
",Threshold="<<lastTH<<
",i="<<i<<
G4endl;
191 G4cout<<
"G4QHyElCS::GetCS:Found T="<<pEn<<
" < Threshold="<<lastTH<<
",CS=0"<<
G4endl;
202 G4cout<<
"G4QHyperonElasticCS::GetCS:P="<<pMom<<
",CS="<<lastCS*millibarn<<
G4endl;
205 return lastCS*millibarn;
210 G4cout<<
"G4QHyElCS::G:UpdateDB P="<<pMom<<
",f="<<fCS<<
",I="<<lastI<<
",i="<<i<<
G4endl;
214 G4cout<<
"G4QHyperElCS::GetCrosSec: *****> New (inDB) Calculated CS="<<lastCS<<
G4endl;
217 if(lastCS<=0. && pEn>lastTH)
220 G4cout<<
"G4QHyperonElCS::GetCS: NewT="<<pEn<<
"(CS=0) > Threshold="<<lastTH<<
G4endl;
227 G4cout<<
"---G4QHyperonElCrossSection::GetCrosSec:pPDG="<<pPDG<<
",i="<<i<<
",N="<<colN[i]
228 <<
",Z["<<i<<
"]="<<colZ[i]<<
G4endl;
235 G4cout<<
"G4QHyElCS::GetCrosSec:CalcNew P="<<pMom<<
",f="<<fCS<<
",lastI="<<lastI<<
G4endl;
243 G4cout<<
"G4QHyperonElasticCrossSection::GetCS:NewThres="<<lastTH<<
",T="<<pEn<<
G4endl;
248 G4cout<<
"G4QHyperonElCS::GetCS:1st T="<<pEn<<
"(CS=0) > Threshold="<<lastTH<<
G4endl;
254 G4cout<<
"G4QHyElCS::GetCrosSec: New CS="<<lastCS<<
",lZ="<<lastN<<
",lN="<<lastZ<<
G4endl;
259 colP.push_back(pMom);
260 colTH.push_back(lastTH);
261 colCS.push_back(lastCS);
263 G4cout<<
"G4QHyElCS::GetCS:1st,P="<<pMom<<
"(MeV),CS="<<lastCS*millibarn<<
"(mb)"<<
G4endl;
266 return lastCS*millibarn;
271 G4cout<<
"G4QHyperonElasticCrossSection::GetCS: Update lastI="<<lastI<<
G4endl;
277 G4cout<<
"G4QHyElCS::GetCSec:End,P="<<pMom<<
"(MeV),CS="<<lastCS*millibarn<<
"(mb)"<<
G4endl;
279 G4cout<<
"G4QHyperonElasticCrossSection::GetCrosSec:***End***,onlyCS="<<onlyCS<<
G4endl;
281 return lastCS*millibarn;
290 static std::vector <G4double> PIN;
295 G4cout<<
"G4QHyperonElasticCS::CalcCS:->onlyCS="<<onlyCS<<
",F="<<F<<
",p="<<pIU<<
G4endl;
297 lastLP=std::log(pMom);
315 G4cout<<
"G4QHyperonElCS::CalcCS: DB's updated for I="<<I<<
",*,PIN4="<<PIN[4]<<
G4endl;
319 G4cout<<
"G4QHyperonElasticCS::CalcCS:*read*, LP="<<lastLP<<
",PIN="<<lastPIN<<
G4endl;
321 if(lastLP>lastPIN && lastLP<lPMax)
323 lastPIN=GetPTables(lastLP,lastPIN,PDG,tgZ,tgN);
325 G4cout<<
"G4QHyElCS::CalcCS:updated(I),LP="<<lastLP<<
"<IN["<<I<<
"]="<<lastPIN<<
G4endl;
345 G4cout<<
"G4QHyperonElasticCrosS::CalcCS:*ini*,lastLP="<<lastLP<<
",min="<<lPMin<<
G4endl;
347 lastPIN = GetPTables(lastLP,lPMin,PDG,tgZ,tgN);
349 G4cout<<
"G4QHypElCS::CCS:i, Z="<<tgZ<<
",N="<<tgN<<
",PDG="<<PDG<<
",LP"<<lastPIN<<
G4endl;
351 PIN.push_back(lastPIN);
352 PAR.push_back(lastPAR);
353 CST.push_back(lastCST);
354 SST.push_back(lastSST);
355 S1T.push_back(lastS1T);
356 B1T.push_back(lastB1T);
357 S2T.push_back(lastS2T);
358 B2T.push_back(lastB2T);
359 S3T.push_back(lastS3T);
360 B3T.push_back(lastB3T);
361 S4T.push_back(lastS4T);
362 B4T.push_back(lastB4T);
366 G4cout<<
"G4QHypElCS::CalCS:?update?,LP="<<lastLP<<
",IN="<<lastPIN<<
",ML="<<lPMax<<
G4endl;
368 if(lastLP>lastPIN && lastLP<lPMax)
370 lastPIN = GetPTables(lastLP,lastPIN,PDG,tgZ,tgN);
372 G4cout<<
"G4QHyperionElCS::CalcCS:*updated(O)*, LP="<<lastLP<<
" < IN="<<lastPIN<<
G4endl;
376 G4cout<<
"G4QHypElCS::CalcCS: lastLP="<<lastLP<<
",lPM="<<lPMin<<
",lPIN="<<lastPIN<<
G4endl;
378 if(!onlyCS) lastTM=GetQ2max(PDG, tgZ, tgN, pMom);
380 G4cout<<
"G4QHyperElCrosSec::CalcCS:oCS="<<onlyCS<<
",-t="<<lastTM<<
", p="<<lastLP<<
G4endl;
382 if(lastLP>lPMin && lastLP<=lastPIN)
386 G4double shift=(lastLP-lPMin)/dlnP+.000001;
387 G4int blast=
static_cast<int>(shift);
388 if(blast<0 || blast>=nLast)
G4cout<<
"G4QHyperElCS::CCS:b="<<blast<<
","<<nLast<<
G4endl;
389 lastSIG = lastCST[blast];
392 theSS = lastSST[blast];
393 theS1 = lastS1T[blast];
394 theB1 = lastB1T[blast];
395 theS2 = lastS2T[blast];
396 theB2 = lastB2T[blast];
397 theS3 = lastS3T[blast];
398 theB3 = lastB3T[blast];
399 theS4 = lastS4T[blast];
400 theB4 = lastB4T[blast];
403 G4cout<<
"G4QHyperonElasticCS::CalculateCS:(E) S1="<<theS1<<
", B1="<<theB1<<
G4endl;
409 G4int blast=
static_cast<int>(shift);
411 if(blast>=nLast) blast=nLast-1;
415 lastSIG= SIGL+shift*(lastCST[lastL]-SIGL);
417 G4cout<<
"G4QHyperonElasticCrossSection::CalcCrossSection:Sig="<<lastSIG<<
",P="
418 <<pMom<<
",Z="<<tgZ<<
",N="<<tgN<<
",PDG="<<PDG<<
",onlyCS="<<onlyCS<<
G4endl;
423 theSS=SSTL+shift*(lastSST[lastL]-SSTL);
425 theS1=S1TL+shift*(lastS1T[lastL]-S1TL);
428 G4cout<<
"G4QHyperonElasticCS::CalcCrossSection:b="<<blast<<
",ls="<<lastL<<
",SL="
429 <<S1TL<<
",SU="<<lastS1T[lastL]<<
",BL="<<B1TL<<
",BU="<<lastB1T[lastL]<<
G4endl;
431 theB1=B1TL+shift*(lastB1T[lastL]-B1TL);
433 theS2=S2TL+shift*(lastS2T[lastL]-S2TL);
435 theB2=B2TL+shift*(lastB2T[lastL]-B2TL);
437 theS3=S3TL+shift*(lastS3T[lastL]-S3TL);
439 G4cout<<
"G4QHyperonElasticCrossSection::CCS:s3l="<<S3TL<<
",sh3="<<shift<<
",s3h="
440 <<lastS3T[lastL]<<
",b="<<blast<<
",l="<<lastL<<
G4endl;
443 theB3=B3TL+shift*(lastB3T[lastL]-B3TL);
445 theS4=S4TL+shift*(lastS4T[lastL]-S4TL);
447 G4cout<<
"G4QHyperonElasticCrossSection::CCS:s4l="<<S4TL<<
",sh4="<<shift<<
",s4h="
448 <<lastS4T[lastL]<<
",b="<<blast<<
",l="<<lastL<<
G4endl;
451 theB4=B4TL+shift*(lastB4T[lastL]-B4TL);
454 G4cout<<
"G4QHyperonElasticCS::CalculateCS:(I) S1="<<theS1<<
", B1="<<theB1<<
G4endl;
458 else lastSIG=GetTabValues(lastLP, PDG, tgZ, tgN);
459 if(lastSIG<0.) lastSIG = 0.;
461 G4cout<<
"G4QHyperonElasticCrossSection::CalculateCS: END, onlyCS="<<onlyCS<<
G4endl;
472 const G4int n_hypel=33;
474 G4double hyp_el[n_hypel]={1.,.002,.12,.0557,3.5,6.72,99.,2.,3.,5.,74.,3.,3.4,.2,.17,
475 .001,8.,.055,3.64,5.e-5,4000.,1500.,.46,1.2e6,3.5e6,5.e-5,
476 1.e10,8.5e8,1.e10,1.1,3.4e6,6.8e6,0.};
479 if(PDG!=3222 && PDG>3000 && PDG<3335)
492 if(lastPAR[nLast]!=pwd)
494 if ( tgZ == 1 && tgN == 0 )
496 for (
G4int ip=0; ip<n_hypel; ip++) lastPAR[ip]=hyp_el[ip];
520 lastPAR[0]=4./(1.+22/asa);
521 lastPAR[1]=2.36*asa/(1.+a*.055/ssa);
522 lastPAR[2]=(1.+.00007*a3/ssa)/(1.+.0026*a2);
523 lastPAR[3]=1.76*a/ssa+.00003*a3;
524 lastPAR[4]=(.03+200./a3)/(1.+1.E5/a3/sa);
535 lastPAR[10]=1.2e7*a8+380*a17;
536 lastPAR[11]=.7/(1.+4.e-12*a16);
537 lastPAR[12]=2.5/a8/(a4+1.e-16*a32);
539 lastPAR[14]=1.2*a2+2.3;
542 lastPAR[16]=.01/(1.+.0024*a5);
544 lastPAR[18]=9.e-7/(1.+.035*a5);
545 lastPAR[19]=(42.+2.7e-11*a16)/(1.+.14*a);
549 lastPAR[22]=2.4e-3*a8/(1.+2.6e-4*a7);
550 lastPAR[23]=3.5e-36*a32*a8/(1.+5.e-15*a32/a);
552 lastPAR[24]=1.e5/(a8+2.5e12/a16);
553 lastPAR[25]=8.e7/(a12+1.e-27*a28*a28);
554 lastPAR[26]=.0006*a3;
556 lastPAR[27]=10.+4.e-8*a12*a;
561 lastPAR[31]=1./(1.+.0001*a8);
562 lastPAR[32]=1.5e-4/(1.+5.e-6*a12);
566 lastPAR[35]=2.e-7*a4;
570 lastPAR[38]=1.e8*std::exp(.32*asa);
571 lastPAR[39]=20.*std::exp(.45*asa);
572 lastPAR[40]=7.e3+2.4e6/a5;
573 lastPAR[41]=2.5e5*std::exp(.085*a3);
576 lastPAR[43]=920.+.03*a8*a3;
577 lastPAR[44]=93.+.0023*a12;
579 G4cout<<
"G4QHyElCS::CalcCS:la "<<lastPAR[38]<<
", "<<lastPAR[39]<<
", "<<lastPAR[40]
580 <<
", "<<lastPAR[42]<<
", "<<lastPAR[43]<<
", "<<lastPAR[44]<<
G4endl;
591 lastPAR[ 9]=4.5*std::pow(a,1.15);
592 lastPAR[10]=.06*std::pow(a,.6);
593 lastPAR[11]=.6*a/(1.+2.e15/a16);
594 lastPAR[12]=.17/(a+9.e5/a3+1.5e33/a32);
595 lastPAR[13]=(.001+7.e-11*a5)/(1.+4.4e-11*a5);
596 lastPAR[14]=(p1a10*p1a10+2.e-29)/(1.+2.e-22*a12);
598 lastPAR[15]=400./a12+2.e-22*a9;
599 lastPAR[16]=1.e-32*a12/(1.+5.e22/a14);
600 lastPAR[17]=1000./a2+9.5*sa*ssa;
601 lastPAR[18]=4.e-6*a*asa+1.e11/a16;
602 lastPAR[19]=(120./a+.002*a2)/(1.+2.e14/a16);
603 lastPAR[20]=9.+100./a;
605 lastPAR[21]=.002*a3+3.e7/a6;
606 lastPAR[22]=7.e-15*a4*asa;
607 lastPAR[23]=9000./a4;
609 lastPAR[24]=.0011*asa/(1.+3.e34/a32/a4);
610 lastPAR[25]=1.e-5*a2+2.e14/a16;
611 lastPAR[26]=1.2e-11*a2/(1.+1.5e19/a12);
612 lastPAR[27]=.016*asa/(1.+5.e16/a16);
614 lastPAR[28]=.002*a4/(1.+7.e7/std::pow(a-6.83,14));
615 lastPAR[29]=2.e6/a6+7.2/std::pow(a,.11);
616 lastPAR[30]=11.*a3/(1.+7.e23/a16/a8);
617 lastPAR[31]=100./asa;
619 lastPAR[32]=(.1+4.4e-5*a2)/(1.+5.e5/a4);
620 lastPAR[33]=3.5e-4*a2/(1.+1.e8/a8);
621 lastPAR[34]=1.3+3.e5/a4;
622 lastPAR[35]=500./(a2+50.)+3;
623 lastPAR[36]=1.e-9/a+s4a16*s4a16;
625 lastPAR[37]=.4*asa+3.e-9*a6;
626 lastPAR[38]=.0005*a5;
630 lastPAR[41]=.05+.005*a;
631 lastPAR[42]=7.e-8/sa;
635 lastPAR[46]=3.e32/(a32+1.e32);
639 lastPAR[49]=7.e3*a/(sa+1.);
640 lastPAR[50]=900.*sa/(1.+500./a3);
642 G4cout<<
"G4QHyElCS::CalcCS:ha "<<lastPAR[41]<<
", "<<lastPAR[42]<<
", "<<lastPAR[43]
643 <<
", "<<lastPAR[44]<<
", "<<lastPAR[45]<<
", "<<lastPAR[46]<<
G4endl;
647 lastPAR[51]=1.e15+2.e27/a4/(1.+2.e-18*a16);
654 lastCST[0]=GetTabValues(lp, PDG, tgZ, tgN);
666 G4cout<<
"G4QHyperonElasticCrossSection::GetPTables:ip=0(init), lp="<<lp<<
",S1="
667 <<theS1<<
",B1="<<theB1<<
",S2="<<theS2<<
",B2="<<theB3<<
",S3="<<theS3
668 <<
",B3="<<theB3<<
",S4="<<theS4<<
",B4="<<theB4<<
G4endl;
673 G4int ini =
static_cast<int>((ILP-lPMin+.000001)/dlnP)+1;
677 G4int fin =
static_cast<int>((LP-lPMin)/dlnP)+1;
678 if(fin>=nPoints) fin=nLast;
682 for(
G4int ip=ini; ip<=fin; ip++)
687 lastCST[ip]=GetTabValues(lp, PDG, tgZ, tgN);
699 G4cout<<
"G4QHyperonElasticCrossSection::GetPTables:ip="<<ip<<
",lp="<<lp
700 <<
",S1="<<theS1<<
",B1="<<theB1<<
",S2="<<theS2<<
",B2="<<theB2<<
",S3="
701 <<theS3<<
",B3="<<theB3<<
",S4="<<theS4<<
",B4="<<theB4<<
G4endl;
706 else G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetPTables: PDG="<<PDG
707 <<
", Z="<<tgZ<<
", N="<<tgN<<
", i="<<ini<<
" > fin="<<fin<<
", LP="<<LP
708 <<
" > ILP="<<ILP<<
" nothing is done!"<<
G4endl;
710 else G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetPTables: PDG="<<PDG
711 <<
", Z="<<tgZ<<
", N="<<tgN<<
", i="<<ini<<
">= max="<<nPoints<<
", LP="<<LP
712 <<
" > ILP="<<ILP<<
", lPMax="<<lPMax<<
" nothing is done!"<<
G4endl;
715 else G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetPTab:PDG="<<PDG<<
",Z="<<tgZ
716 <<
", N="<<tgN<<
", LP="<<LP<<
" <= ILP="<<ILP<<
" nothing is done!"<<
G4endl;
723 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
724 <<
", while it is defined only for Hyperons" <<
G4endl;
725 G4Exception(
"G4QHyperonElasticCrossSection::GetPTables()",
"HAD_CHPS_0000",
734 static const G4double GeVSQ=gigaelectronvolt*gigaelectronvolt;
739 G4cout<<
"G4QHyperElCS::GetExcT:F="<<onlyCS<<
",Z="<<tgZ<<
",N="<<tgN<<
",PDG="<<PDG<<
G4endl;
741 if(PDG==3222 || PDG<3000 || PDG>3334)
G4cout<<
"*Warning*G4QHyElCS::GET:PDG="<<PDG<<
G4endl;
742 if(onlyCS)
G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetExchanT: onlyCS=1"<<
G4endl;
748 G4cout<<
"G4QHyperElCS::GetExchangeT: TM="<<lastTM<<
",S1="<<theS1<<
",B1="<<theB1<<
",S2="
749 <<theS2<<
",B2="<<theB2<<
",S3="<<theS3<<
",B3="<<theB3<<
",GeV2="<<GeVSQ<<
G4endl;
754 G4double ts1=-std::log(1.-R1)/theB1;
755 G4double ds1=std::fabs(ts1-lastTM)/lastTM;
757 G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetExT:1p "<<ts1<<
"#"<<lastTM
758 <<
",d="<<ds1<<
",R1="<<R1<<
",E1="<<E1<<
G4endl;
761 G4double R2=(1.-std::exp(-E2*E2*E2));
763 G4double ts2=std::pow(-std::log(1.-R2),.333333333)/theB2;
764 G4double ds2=std::fabs(ts2-lastTM)/lastTM;
766 G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetExT:2p "<<ts2<<
"#"<<lastTM
767 <<
",d="<<ds2<<
",R2="<<R2<<
",E2="<<E2<<
G4endl;
772 G4double ts3=-std::log(1.-R3)/theB3;
773 G4double ds3=std::fabs(ts3-lastTM)/lastTM;
775 G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetExT:3p "<<ts3<<
"#"<<lastTM
776 <<
",d="<<ds3<<
",R3="<<R1<<
",E3="<<E3<<
G4endl;
787 q2=-std::log(1.-ran)/theB1;
793 q2=-std::log(1.-ran);
795 q2=std::pow(q2,third)/theB2;
801 q2=-std::log(1.-ran)/theB3;
808 G4cout<<
"G4QHyperonElasticCrossSection::GetExT:a="<<a<<
",t="<<lastTM<<
",S1="<<theS1
809 <<
",B1="<<theB1<<
",SS="<<theSS<<
",S2="<<theS2<<
",B2="<<theB2<<
",S3="<<theS3
810 <<
",B3="<<theB3<<
",S4="<<theS4<<
",B4="<<theB4<<
G4endl;
812 G4double E1=lastTM*(theB1+lastTM*theSS);
816 G4double ts1=-std::log(1.-R1)/theB1;
817 if(std::fabs(tss)>1.e-7) ts1=(std::sqrt(theB1*(theB1+(tss+tss)*ts1))-theB1)/tss;
820 G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetExT:1a "<<ts1<<
"#"<<lastTM
821 <<
",d="<<ds1<<
",R1="<<R1<<
",E1="<<E1<<
G4endl;
828 G4double ts2=-std::log(1.-R2)/theB2;
829 if(a<6.5)ts2=std::pow(ts2,third);
830 else ts2=std::pow(ts2,fifth);
831 G4double ds2=std::fabs(ts2-lastTM)/lastTM;
833 G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetExT:2a "<<ts2<<
"#"<<lastTM
834 <<
",d="<<ds2<<
",R2="<<R2<<
",E2="<<E2<<
G4endl;
837 if(a>6.5)E3*=tm2*tm2*tm2;
840 G4double ts3=-std::log(1.-R3)/theB3;
841 if(a>6.5)ts3=std::pow(ts3,sevth);
842 G4double ds3=std::fabs(ts3-lastTM)/lastTM;
844 G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetExT:3a "<<ts3<<
"#"<<lastTM
845 <<
",d="<<ds3<<
",R3="<<R3<<
",E3="<<E3<<
G4endl;
850 G4double ts4=-std::log(1.-R4)/theB4;
851 G4double ds4=std::fabs(ts4-lastTM)/lastTM;
853 G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetExT:4a "<<ts4<<
"#"<<lastTM
854 <<
",d="<<ds4<<
",R4="<<R4<<
",E4="<<E4<<
G4endl;
864 G4cout<<
"G4QHyElCS::GExT:1="<<I1<<
",2="<<I2<<
",3="<<I3<<
",4="<<I4<<
",r="<<rand<<
G4endl;
870 q2=-std::log(1.-ran)/theB1;
871 if(std::fabs(tss)>1.e-7) q2=(std::sqrt(theB1*(theB1+(tss+tss)*q2))-theB1)/tss;
873 G4cout<<
"G4QHyElCS::GExT:Q2="<<q2<<
",ss="<<tss/2<<
",b1="<<theB1<<
",t1="<<ts1<<
G4endl;
880 q2=-std::log(1.-ran)/theB2;
882 if(a<6.5) q2=std::pow(q2,third);
883 else q2=std::pow(q2,fifth);
885 G4cout<<
"G4QHyElCS::GetExT: Q2="<<q2<<
",r2="<<R2<<
",b2="<<theB2<<
",t2="<<ts2<<
G4endl;
892 q2=-std::log(1.-ran)/theB3;
894 if(a>6.5) q2=std::pow(q2,sevth);
896 G4cout<<
"G4QHyElCS::GetExT: Q2="<<q2<<
",r3="<<R2<<
",b3="<<theB2<<
",t3="<<ts2<<
G4endl;
903 q2=-std::log(1.-ran)/theB4;
904 if(a<6.5) q2=lastTM-q2;
906 G4cout<<
"G4QHyElCS::GExT:Q2="<<q2<<
",m="<<lastTM<<
",b4="<<theB3<<
",t4="<<ts3<<
G4endl;
911 if(!(q2>=-1.||q2<=1.))
G4cout<<
"*NAN*G4QHyElasticCrossSect::GetExchangeT:-t="<<q2<<
G4endl;
915 G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GExT: -t="<<q2<<
" > "<<lastTM<<
G4endl;
925 static const G4double GeVSQ=gigaelectronvolt*gigaelectronvolt;
927 G4cout<<
"G4QHyperElCS::GetSlope:"<<onlyCS<<
", Z="<<tgZ<<
",N="<<tgN<<
",PDG="<<PDG<<
G4endl;
929 if(onlyCS)
G4cout<<
"*Warning*G4QHyperonElasticCrossSection::GetSlope: onlCS=true"<<
G4endl;
930 if(lastLP<-4.3)
return 0.;
931 if(PDG==3222 || PDG<3000 || PDG>3334)
937 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
938 <<
", while it is defined only for Hyperons" <<
G4endl;
939 G4Exception(
"G4QHyperonElasticCrossSection::GetSlope()",
"HAD_CHPS_0000",
942 if(theB1<0.) theB1=0.;
943 if(!(theB1>=-1.||theB1<=1.))
G4cout<<
"*NAN*G4QHyElasticCrossS::Getslope:"<<theB1<<
G4endl;
950 static const G4double HGeVSQ=gigaelectronvolt*gigaelectronvolt/2.;
951 return lastTM*HGeVSQ;
958 if(PDG==3222 || PDG<3000 || PDG>3334)
G4cout<<
"*Warning*G4QHypElCS::GTV:P="<<PDG<<
G4endl;
961 G4cout<<
"*Warning*G4QHyperonElastCS::GetTabValue:(1-92) NoIsotopesFor Z="<<tgZ<<
G4endl;
979 G4cout<<
"G4QHyElasticCS::GetTabVal:l="<<lp<<
",Z="<<tgZ<<
",N="<<tgN<<
",PDG="<<PDG<<
G4endl;
986 if ( tgZ == 1 && tgN == 0 )
990 theS1=(lastPAR[10]+lastPAR[11]*dl2*dl2)/(1.+lastPAR[12]/p4/p)+
991 (lastPAR[13]/p2+lastPAR[14]*p)/(p4+lastPAR[15]*sp);
992 theB1=lastPAR[16]*std::pow(p,lastPAR[17])/(1.+lastPAR[18]/p3);
993 theS2=lastPAR[19]+lastPAR[20]/(p4+lastPAR[21]*p);
994 theB2=lastPAR[22]+lastPAR[23]/(p4+lastPAR[24]/
sp);
995 theS3=lastPAR[25]+lastPAR[26]/(p4*p4+lastPAR[27]*p2+lastPAR[28]);
996 theB3=lastPAR[29]+lastPAR[30]/(p4+lastPAR[31]);
1000 G4cout<<
"G4QHyperonElasticCrossSection::GetTableVal:TM="<<lastTM<<
",S1="<<theS1<<
",B1="
1001 <<theB1<<
",S2="<<theS2<<
",B2="<<theB2<<
",S3="<<theS1<<
",B3="<<theB1<<
G4endl;
1005 return lastPAR[0]/(lastPAR[1]+p2*(lastPAR[2]+p2))+(lastPAR[3]*dp*dp+lastPAR[5]+
1006 lastPAR[6]/p2)/(1.+lastPAR[7]/sp+lastPAR[8]/p4);
1024 theS1=lastPAR[9]/(1.+lastPAR[10]*p4*pa)+lastPAR[11]/(p4+lastPAR[12]*p4/pa2)+
1025 (lastPAR[13]*dl*dl+lastPAR[14])/(1.+lastPAR[15]/p2);
1026 theB1=(lastPAR[16]+lastPAR[17]*p2)/(p4+lastPAR[18]/pah)+lastPAR[19];
1027 theSS=lastPAR[20]/(1.+lastPAR[21]/p2)+lastPAR[22]/(p6/pa+lastPAR[23]/p16);
1028 theS2=lastPAR[24]/(pa/p2+lastPAR[25]/p4)+lastPAR[26];
1029 theB2=lastPAR[27]*std::pow(p,lastPAR[28])+lastPAR[29]/(p8+lastPAR[30]/p16);
1030 theS3=lastPAR[31]/(pa*p+lastPAR[32]/pa)+lastPAR[33];
1031 theB3=lastPAR[34]/(p3+lastPAR[35]/p6)+lastPAR[36]/(1.+lastPAR[37]/p2);
1032 theS4=p2*(pah*lastPAR[38]*std::exp(-pah*lastPAR[39])+
1033 lastPAR[40]/(1.+lastPAR[41]*std::pow(p,lastPAR[42])));
1034 theB4=lastPAR[43]*pa/p2/(1.+pa*lastPAR[44]);
1036 G4cout<<
"G4QHyperonElasticCS::GetTabV: lA, p="<<p<<
",S1="<<theS1<<
",B1="<<theB1
1037 <<
",SS="<<theSS<<
",S2="<<theS2<<
",B2="<<theB2<<
",S3="<<theS3<<
",B3="<<theB3
1038 <<
",S4="<<theS4<<
",B4="<<theB4<<
G4endl;
1043 theS1=lastPAR[9]/(1.+lastPAR[10]/p4)+lastPAR[11]/(p4+lastPAR[12]/p2)+
1044 lastPAR[13]/(p5+lastPAR[14]/p16);
1045 theB1=(lastPAR[15]/p8+lastPAR[19])/(p+lastPAR[16]/std::pow(p,lastPAR[20]))+
1046 lastPAR[17]/(1.+lastPAR[18]/p4);
1047 theSS=lastPAR[21]/(p4/std::pow(p,lastPAR[23])+lastPAR[22]/p4);
1048 theS2=lastPAR[24]/p4/(std::pow(p,lastPAR[25])+lastPAR[26]/p12)+lastPAR[27];
1049 theB2=lastPAR[28]/std::pow(p,lastPAR[29])+lastPAR[30]/std::pow(p,lastPAR[31]);
1050 theS3=lastPAR[32]/std::pow(p,lastPAR[35])/(1.+lastPAR[36]/p12)+
1051 lastPAR[33]/(1.+lastPAR[34]/p6);
1052 theB3=lastPAR[37]/p8+lastPAR[38]/p2+lastPAR[39]/(1.+lastPAR[40]/p8);
1053 theS4=(lastPAR[41]/p4+lastPAR[46]/p)/(1.+lastPAR[42]/p10)+
1054 (lastPAR[43]+lastPAR[44]*dl*dl)/(1.+lastPAR[45]/p12);
1055 theB4=lastPAR[47]/(1.+lastPAR[48]/p)+lastPAR[49]*p4/(1.+lastPAR[50]*p5);
1057 G4cout<<
"G4QHyperonElasticCS::GetTabV:hA,p="<<p<<
",S1="<<theS1<<
",B1="<<theB1<<
",SS="
1058 <<theSS<<
",S2="<<theS2<<
",B2="<<theB2<<
",S3="<<theS3<<
",B3="<<theB3<<
",S4="
1059 <<theS4<<
",B4="<<theB4<<
G4endl;
1064 G4cout<<
"G4QHyperonElCS::GetTabV: PDG="<<PDG<<
",P="<<p<<
",N="<<tgN<<
",Z="<<tgZ<<
G4endl;
1068 return (lastPAR[0]*dlp*dlp+lastPAR[1])/(1.+lastPAR[2]/p)+lastPAR[3]/(p3+lastPAR[4]);
1085 static const G4double mLa2= mLamb*mLamb;
1094 G4double s_value=dmt*std::sqrt(pP2+mLa2)+mLa2+mt*mt;
1095 return dmt*dmt*pP2/s_value;
1103 ed <<
"PDG = " << PDG <<
", Z = " << tgZ <<
", N = " << tgN
1104 <<
", while it is defined only for p projectiles & Z_target>0" <<
G4endl;
1105 G4Exception(
"G4QHyperonElasticCrossSection::GetQ2max()",
"HAD_CHPS_0000",
G4DLLIMPORT std::ostream G4cout
G4QHyperonElasticCrossSection()
G4double GetExchangeT(G4int tZ, G4int tN, G4int pPDG)
static G4VQCrossSection * GetPointer()
G4double CalculateCrossSection(G4bool CS, G4int F, G4int I, G4int pPDG, G4int Z, G4int N, G4double pP)
virtual G4double GetCrossSection(G4bool fCS, G4double pMom, G4int tgZ, G4int tgN, G4int pPDG=3122)
G4double GetSlope(G4int tZ, G4int tN, G4int pPDG)
~G4QHyperonElasticCrossSection()
virtual G4double ThresholdEnergy(G4int Z, G4int N, G4int PDG=0)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
std::ostringstream G4ExceptionDescription