45 #ifndef G4ElasticHadrNucleusHE_h
46 #define G4ElasticHadrNucleusHE_h 1
241 return c1+(p-
p1)*(c2-c1)/(p2-
p1);
255 if(i == n) i = n - 1;
271 if ( numN >= numM && numN <= 240)
return SetBinom[numN][numM];
G4double LineInterpol(G4double p0, G4double p2, G4double c1, G4double c2, G4double p)
static c2_factory< G4double > c2
G4int HadronType[NHADRONS]
G4double SampleT(const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
G4double GetQ2(G4double Ran)
G4NistManager * nistManager
G4double GetBinomCof(G4int n, G4int m)
void DefineHadronValues(G4int Z)
G4double GetDistrFun(G4double Q2)
const G4ParticleDefinition * Hadron()
G4double HadronProtonQ2(const G4ParticleDefinition *aHadron, G4double inLabMom)
G4double LowEdgeEnergy[NENERGY]
void GetKinematics(const G4ParticleDefinition *aHadron, G4double MomentumH)
G4double SetBinom[240][240]
G4ElasticHadrNucleusHE(const G4String &name="hElasticGlauber")
G4double TableCrossSec[NQTABLE]
G4double HadrNucDifferCrSec(G4int Z, G4int Nucleus, G4double Q2)
G4double lowestEnergyLimit
virtual void Description() const
static const G4double A[nN]
G4ElasticHadrNucleusHE & operator=(const G4ElasticHadrNucleusHE &right)
static const G4int NENERGY
virtual ~G4ElasticHadrNucleusHE()
static const G4int NQTABLE
G4ElasticData(const G4ParticleDefinition *h, G4int Z, G4double A, G4double *eGeV)
G4int HadronType1[NHADRONS]
G4double GetLightFq2(G4int Z, G4int A, G4double Q)
G4double GetHeavyFq2(G4int Z, G4int Nucleus, G4double *LineFq2)
G4double GetQ2_2(G4int N, G4double *Q, G4double *F, G4double R)
G4ElasticData & operator=(const G4ElasticData &right)
const G4ParticleDefinition * hadr
G4double HadronNucleusQ2_2(G4ElasticData *pElD, G4int Z, G4double plabGeV, G4double tmax)
void InterpolateHN(G4int n, const G4double EnP[], const G4double C0P[], const G4double C1P[], const G4double B0P[], const G4double B1P[])
G4double GetFt(G4double Q2)
void DefineNucleusParameters(G4double A)
static const G4int NHADRONS
virtual G4double SampleInvariantT(const G4ParticleDefinition *p, G4double plab, G4int Z, G4int A)
G4ElasticData * SetOfElasticData[NHADRONS][93]
G4int HadronCode[NHADRONS]