60 #ifndef G4MuPairProductionModel_h
61 #define G4MuPairProductionModel_h 1
133 void MakeSamplingTables();
135 inline G4double InterpolatedIntegralCrossSection(
172 static G4double zdat[5], adat[5], tdat[8];
173 G4double ya[1001], proba[5][8][1001];
179 G4bool samplingTablesAreFilled;
180 std::vector<G4double> partialSum;
216 inline G4double G4MuPairProductionModel::InterpolatedIntegralCrossSection(
223 (fac*proba[
iz][it-1][iy]-fac1*proba[
iz-1][it-1][iy])*dz;
225 (fac*proba[
iz][it ][iy]-fac1*proba[
iz-1][it ][iy])*dz;
226 return (f0 + (
f1-f0)*dt)*
z*(
z+1.);