36 #ifndef G4MuMinusCaptureCascade_h
37 #define G4MuMinusCaptureCascade_h 1
95 if(Xuser <= X[0]) Yuser = Y[0];
96 else if(Xuser >= X[N-1]) Yuser = Y[N-1];
100 if(Xuser <= X[i])
break;
103 if(Xuser == X[i]) Yuser = Y[i];
104 else Yuser = Y[i-1] + (Y[i] - Y[i-1])*(Xuser - X[i-1])/(X[i] - X[i-1]);
117 G4double sint = std::sqrt((1.0 - cost)*(1.0 + cost));
119 randomVect =
G4ThreeVector(sint * std::cos(Phi), sint * std::sin(Phi), cost);
CLHEP::Hep3Vector G4ThreeVector
G4ThreeVector & GetRandomVec()
G4MuMinusCaptureCascade()
G4int DoCascade(const G4double Z, const G4double A, G4GHEKinematicsVector *Cascade)
G4double GetKShellEnergy(G4double Z)
~G4MuMinusCaptureCascade()
void DoBoundMuonMinusDecay(G4double Z, G4int *nCascade, G4GHEKinematicsVector *Cascade)