68 outFile <<
"G4NeutrinoElectronNcModel is a neutrino-electron (neutral current) elastic scattering\n"
69 <<
"model which uses the standard model \n"
70 <<
"transfer parameterization. The model is fully relativistic\n";
86 if( ( pName ==
"nu_e" || pName ==
"anti_nu_e" ||
87 pName ==
"nu_mu" || pName ==
"anti_nu_nu" ||
88 pName ==
"nu_tau" || pName ==
"anti_nu_tau" ) &&
125 if( cost2 > 1. ) cost2 = 1.;
126 if( cost2 < 0. ) cost2 = 0.;
129 G4double sint = std::sqrt( (1.0 - cost)*(1.0 + cost) );
132 G4ThreeVector eP( sint*std::cos(phi), sint*std::sin(phi), cost );
147 else if( eTkin > 0.0 )
172 G4double result = 0., xi, cofL, cofR, cofL2, cofR2, cofLR;
175 if(
energy == 0.)
return result;
184 else if( pName ==
"anti_nu_e")
189 else if( pName ==
"nu_mu")
194 else if( pName ==
"anti_nu_mu")
199 else if( pName ==
"nu_tau")
204 else if( pName ==
"anti_nu_tau")
229 G4double d = -( a*xMax3 + b*xMax2 + c*xMax );
240 q += 2*b*b*b/a/a/a/27.;
G4double B(G4double temperature)
G4double D(G4double temp)
static constexpr double eV
static constexpr double GeV
std::complex< G4double > G4complex
static G4Electron * Electron()
void AddSecondary(G4DynamicParticle *aP, G4int mod=-1)
void SetEnergyChange(G4double anEnergy)
void SetMomentumChange(const G4ThreeVector &aV)
void SetLocalEnergyDeposit(G4double aE)
const G4ParticleDefinition * GetDefinition() const
G4double GetKineticEnergy() const
const G4LorentzVector & Get4Momentum() const
G4double GetTotalEnergy() const
G4double LowestEnergyLimit() const
void SetLowestEnergyLimit(G4double value)
G4HadFinalState theParticleChange
void SetMinEnergy(G4double anEnergy)
void SetMaxEnergy(const G4double anEnergy)
G4double GetMaxEnergy() const
static G4HadronicParameters * Instance()
G4NeutrinoElectronNcModel(const G4String &name="nu-e-elastic")
virtual void ModelDescription(std::ostream &) const
G4ParticleDefinition * theElectron
G4double SampleElectronTkin(const G4HadProjectile *aParticle)
virtual G4bool IsApplicable(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
virtual ~G4NeutrinoElectronNcModel()
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
G4double GetPDGMass() const
const G4String & GetParticleName() const
static G4int GetModelID(const G4int modelIndex)
static constexpr double twopi
G4double energy(const ThreeVector &p, const G4double m)
const char * name(G4int ptype)