#include <G4ChipsHyperonElasticXS.hh>
|
| G4ChipsHyperonElasticXS () |
|
| ~G4ChipsHyperonElasticXS () |
|
virtual void | CrossSectionDescription (std::ostream &) 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) |
|
virtual G4double | GetChipsCrossSection (G4double momentum, G4int Z, G4int N, G4int pdg) |
|
G4double | GetExchangeT (G4int tZ, G4int tN, G4int pPDG) |
|
| G4VCrossSectionDataSet (const G4String &nam="") |
|
virtual | ~G4VCrossSectionDataSet () |
|
virtual G4bool | IsElementApplicable (const G4DynamicParticle *, G4int Z, const G4Material *mat=0) |
|
G4double | GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0) |
|
G4double | ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=0) |
|
virtual G4double | GetElementCrossSection (const G4DynamicParticle *, G4int Z, const G4Material *mat=0) |
|
virtual G4Isotope * | SelectIsotope (const G4Element *, G4double kinEnergy) |
|
virtual void | BuildPhysicsTable (const G4ParticleDefinition &) |
|
virtual void | DumpPhysicsTable (const G4ParticleDefinition &) |
|
virtual G4int | GetVerboseLevel () const |
|
virtual void | SetVerboseLevel (G4int value) |
|
G4double | GetMinKinEnergy () const |
|
void | SetMinKinEnergy (G4double value) |
|
G4double | GetMaxKinEnergy () const |
|
void | SetMaxKinEnergy (G4double value) |
|
const G4String & | GetName () const |
|
Definition at line 46 of file G4ChipsHyperonElasticXS.hh.
G4ChipsHyperonElasticXS::G4ChipsHyperonElasticXS |
( |
| ) |
|
Definition at line 64 of file G4ChipsHyperonElasticXS.cc.
68 dlnP=(lPMax-lPMin)/nLast;
G4VCrossSectionDataSet(const G4String &nam="")
static const char * Default_Name()
G4ChipsHyperonElasticXS::~G4ChipsHyperonElasticXS |
( |
| ) |
|
Definition at line 104 of file G4ChipsHyperonElasticXS.cc.
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++)
static const G4double pos
void G4ChipsHyperonElasticXS::CrossSectionDescription |
( |
std::ostream & |
outFile | ) |
const |
|
virtual |
Reimplemented from G4VCrossSectionDataSet.
Definition at line 143 of file G4ChipsHyperonElasticXS.cc.
145 outFile <<
"G4ChipsHyperonElasticXS provides the elastic cross\n"
146 <<
"section for hyperon nucleus scattering as a function of incident\n"
147 <<
"momentum. The cross section is calculated using M. Kossov's\n"
148 <<
"CHIPS parameterization of cross section data.\n";
static const char* G4ChipsHyperonElasticXS::Default_Name |
( |
| ) |
|
|
inlinestatic |
!The slave functions must provide cross-sections in millibarns (mb) !! (not in IU)
Definition at line 172 of file G4ChipsHyperonElasticXS.cc.
185 if(lastI)
for(
G4int i=0; i<lastI; i++)
187 if(colN[i]==tgN && colZ[i]==tgZ)
200 CalculateCrossSection(fCS,-1,i,pPDG,lastZ,lastN,pMom);
205 lastCS=CalculateCrossSection(fCS,-1,i,pPDG,lastZ,lastN,pMom);
206 if(lastCS<=0. && pEn>lastTH)
216 lastCS=CalculateCrossSection(fCS,0,lastI,pPDG,lastZ,lastN,pMom);
227 colP.push_back(pMom);
228 colTH.push_back(lastTH);
229 colCS.push_back(lastCS);
static constexpr double millibarn
Definition at line 610 of file G4ChipsHyperonElasticXS.cc.
617 if(PDG<3000 || PDG>3334)
G4cout<<
"*Warning*G4QHyElCS::GET:PDG="<<PDG<<
G4endl;
618 if(onlyCS)
G4cout<<
"*Warning*G4ChipsHyperonElasticXS::GetExchanT: onlyCS=1"<<
G4endl;
638 q2=-
G4Log(1.-ran)/theB1;
652 q2=-
G4Log(1.-ran)/theB3;
658 G4double E1=lastTM*(theB1+lastTM*theSS);
666 if(a>6.5)E3*=tm2*tm2*tm2;
681 q2=-
G4Log(1.-ran)/theB1;
682 if(std::fabs(tss)>1.e-7) q2=(std::sqrt(theB1*(theB1+(tss+tss)*q2))-theB1)/tss;
688 q2=-
G4Log(1.-ran)/theB2;
697 q2=-
G4Log(1.-ran)/theB3;
705 q2=-
G4Log(1.-ran)/theB4;
706 if(a<6.5) q2=lastTM-q2;
710 if(!(q2>=-1.||q2<=1.))
G4cout<<
"*NAN*G4QHyElasticCrossSect::GetExchangeT:-t="<<q2<<
G4endl;
static G4Pow * GetInstance()
G4double powA(G4double A, G4double y) const
std::vector< ExP01TrackerHit * > a
G4GLOB_DLL std::ostream G4cout
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
static constexpr double gigaelectronvolt
Reimplemented from G4VCrossSectionDataSet.
Definition at line 160 of file G4ChipsHyperonElasticXS.cc.
G4int GetPDGEncoding() const
G4ParticleDefinition * GetDefinition() const
G4double GetTotalMomentum() const
double A(double temperature)
virtual G4double GetChipsCrossSection(G4double momentum, G4int Z, G4int N, G4int pdg)
The documentation for this class was generated from the following files: