65 leptonKE(0.0), photonEnergy(0.0), photonQ2(0.0)
106 outFile <<
"G4ElectroVDNuclearModel handles the inelastic scattering\n"
107 <<
"of e- and e+ from nuclei using the equivalent photon\n"
108 <<
"approximation in which the incoming lepton generates a\n"
109 <<
"virtual photon at the electromagnetic vertex, and the\n"
110 <<
"virtual photon is converted to a real photon. At low\n"
111 <<
"energies, the photon interacts directly with the nucleus\n"
112 <<
"using the Bertini cascade. At high energies the photon\n"
113 <<
"is converted to a pi0 which interacts using the FTFP\n"
114 <<
"model. The electro- and gamma-nuclear cross sections of\n"
115 <<
"M. Kossov are used to generate the virtual photon spectrum\n";
184 G4double iniP = std::sqrt(iniE*iniE-mProj2);
185 G4double finP = std::sqrt(finE*finE-mProj2);
187 if (cost > 1.) cost= 1.;
188 if (cost < -1.) cost=-1.;
189 G4double sint = std::sqrt(1.-cost*cost);
203 photonEnergy, photonMomentum);
215 if (gammaE < 10*
GeV) {
221 G4double piMom = std::sqrt(gammaE*gammaE - piMass*piMass);
CLHEP::Hep3Vector G4ThreeVector
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
void SetMinEnergy(G4double anEnergy)
void SetMaxEnergy(const G4double anEnergy)
void CalculateHadronicVertex(G4DynamicParticle *incident, G4Nucleus &target)
void SetKineticEnergy(G4double aEnergy)
void SetMomentumChange(const G4ThreeVector &aV)
G4LundStringFragmentation * theFragmentation
const G4ThreeVector & GetMomentumDirection() const
G4double GetEquivalentPhotonEnergy()
static G4CrossSectionDataSetRegistry * Instance()
static G4Proton * Proton()
static constexpr double PeV
G4double GetPDGMass() const
void SetEnergyChange(G4double anEnergy)
void SetDeExcitation(G4VPreCompoundModel *ptr)
G4HadFinalState * ApplyYourself(const G4HadProjectile &thePrimary, G4Nucleus &theNucleus)
void SetTransport(G4VIntraNuclearTransportModel *const value)
static G4HadronicInteractionRegistry * Instance()
G4DynamicParticle * CalculateEMVertex(const G4HadProjectile &aTrack, G4Nucleus &aTargetNucleus)
static const char * Default_Name()
static constexpr double twopi
virtual G4double GetElementCrossSection(const G4DynamicParticle *, G4int Z, const G4Material *mat)
G4ExcitedStringDecay * theStringDecay
G4double GetEquivalentPhotonQ2(G4double nu)
virtual G4double GetElementCrossSection(const G4DynamicParticle *, G4int Z, const G4Material *)
~G4ElectroVDNuclearModel()
const G4LorentzVector & Get4Momentum() const
Hep3Vector cross(const Hep3Vector &) const
G4double GetKineticEnergy() const
G4ElectroVDNuclearModel()
const G4ParticleDefinition * GetDefinition() const
G4ElectroNuclearCrossSection * electroXS
G4CascadeInterface * bert
virtual void ModelDescription(std::ostream &outFile) const
void SetHighEnergyGenerator(G4VHighEnergyGenerator *const value)
static const char * Default_Name()
static G4Neutron * Neutron()
G4double GetVirtualFactor(G4double nu, G4double Q2)
G4PhotoNuclearCrossSection * gammaXS
static G4PionZero * PionZero()
void SetFragmentationModel(G4VStringFragmentation *aModel)
Hep3Vector orthogonal() const
G4double GetTotalEnergy() const
G4HadronicInteraction * FindModel(const G4String &name)
void AddSecondaries(const std::vector< G4HadSecondary > &addSecs)
static constexpr double GeV
G4HadFinalState theParticleChange
G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &theNucleus)
void SetStatusChange(G4HadFinalStateStatus aS)