65 G4bool IsScatProjToProjCase,
87 IsScatProjToProjCase);
93 G4double gammaE2 = adjointPrimKinEnergy;
94 if (!IsScatProjToProjCase) gammaE2 = gammaE1 - adjointPrimKinEnergy;
104 G4double cos_th = 1.+ electron_mass_c2*(1./gammaE1 -1./gammaE2);
105 if (!IsScatProjToProjCase) {
107 cos_th = (gammaE1 - gammaE2*cos_th)/p_elec;
110 if (std::abs(cos_th)>1){
118 else sin_th = std::sqrt(1.-cos_th*cos_th);
130 gammaMomentum1.rotateUz(dir_parallel);
138 adjointPrimKinEnergy,
140 IsScatProjToProjCase);
142 if (!IsScatProjToProjCase){
157 G4bool IsScatProjToProjCase,
177 if (!IsScatProjToProjCase){
181 if (Emin>=Emax)
return;
183 G4double f2=(Emax-adjointPrimKinEnergy)/Emax/f1;
184 gammaE1=adjointPrimKinEnergy/(1.-f1*std::pow(f2,
G4UniformRand()));;
185 gammaE2=gammaE1-adjointPrimKinEnergy;
186 diffCSUsed= diffCSUsed*(1.+2.*std::log(1.+electron_mass_c2/adjointPrimKinEnergy))*adjointPrimKinEnergy/gammaE1/gammaE2;
192 if (Emin>=Emax)
return;
193 gammaE2 =adjointPrimKinEnergy;
195 diffCSUsed= diffCSUsed/gammaE1;
217 w_corr*=diffCS/diffCSUsed;
229 G4double cos_th = 1.+ electron_mass_c2*(1./gammaE1 -1./gammaE2);
230 if (!IsScatProjToProjCase) {
232 cos_th = (gammaE1 - gammaE2*cos_th)/p_elec;
235 if (std::abs(cos_th)>1){
243 else sin_th = std::sqrt(1.-cos_th*cos_th);
255 gammaMomentum1.rotateUz(dir_parallel);
260 if (!IsScatProjToProjCase){
284 G4double gamEnergy1 = gamEnergy0 - kinEnergyElec;
310 G4double gamEnergy1_max = gamEnergy0;
311 G4double gamEnergy1_min = gamEnergy0/one_plus_two_epsi;
312 if (gamEnergy1 >gamEnergy1_max || gamEnergy1<gamEnergy1_min) {
322 G4double one_plus_two_epsi_2=one_plus_two_epsi*one_plus_two_epsi;
326 CS+=4./epsilon +0.5*(1.-1./one_plus_two_epsi_2);
334 G4double epsilon1 = gamEnergy1 / electron_mass_c2 ;
336 G4double term1 =1.+ 1./epsilon -1/epsilon1;
337 G4double dCS_dE1= 1./v +v + term1*term1 -1.;
338 dCS_dE1 *=1./epsilon/gamEnergy0;
361 {
G4double inv_e_max = 1./PrimAdjEnergy - 2./electron_mass_c2;
370 G4double term=std::sqrt(half_e*(electron_mass_c2+half_e));
378 G4bool IsScatProjToProjCase)
387 if (!IsScatProjToProjCase ){
390 if (Emax_proj>Emin_proj ){
391 Cross= 0.1*std::log((Emax_proj-
float (primEnergy))*Emin_proj/Emax_proj/(Emin_proj-primEnergy))
392 *(1.+2.*std::log(
float(1.+electron_mass_c2/primEnergy)));
398 if (Emax_proj>Emin_proj) {
399 Cross = 0.1*std::log(Emax_proj/Emin_proj);
408 return double(Cross);
414 G4bool IsScatProjToProjCase)
static G4AdjointGamma * AdjointGamma()
G4double GetLambda(G4double &kinEnergy, const G4MaterialCutsCouple *couple)
virtual G4double GetSecondAdjEnergyMaxForProdToProjCase(G4double PrimAdjEnergy)
virtual void SampleSecondaries(const G4Track &aTrack, G4bool IsScatProjToProjCase, G4ParticleChange *fParticleChange)
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
G4double GetPostStepWeightCorrection()
const G4DynamicParticle * GetDynamicParticle() const
G4double SampleAdjSecEnergyFromCSMatrix(size_t MatrixIndex, G4double prim_energy, G4bool IsScatProjToProjCase)
G4ParticleDefinition * theDirectPrimaryPartDef
void ProposeParentWeight(G4double finalWeight)
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
virtual G4double AdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)
void RapidSampleSecondaries(const G4Track &aTrack, G4bool IsScatProjToProjCase, G4ParticleChange *fParticleChange)
static G4AdjointElectron * AdjointElectron()
virtual G4double GetAdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)
G4Material * currentMaterial
G4ParticleDefinition * theAdjEquivOfDirectPrimPartDef
G4double GetTotalMomentum() const
void SetUseMatrixPerElement(G4bool aBool)
void SetSecondaryWeightByProcess(G4bool)
void SetParentWeightByProcess(G4bool)
G4VEmModel * theDirectEMModel
double A(double temperature)
G4double GetElectronDensity() const
const G4ThreeVector & GetMomentumDirection() const
virtual G4double GetSecondAdjEnergyMaxForScatProjToProjCase(G4double PrimAdjEnergy)
void SetUseMatrix(G4bool aBool)
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
G4MaterialCutsCouple * currentCouple
void DefineCurrentMaterial(const G4MaterialCutsCouple *couple)
void SetUseOnlyOneMatrixForAllElements(G4bool aBool)
virtual G4double AdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)
virtual G4double GetSecondAdjEnergyMinForScatProjToProjCase(G4double PrimAdjEnergy, G4double Tcut=0)
void ProposeEnergy(G4double finalEnergy)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
void AddSecondary(G4Track *aSecondary)
G4double GetWeight() const
static const G4double Emin
static const G4double Emax
G4ParticleDefinition * theAdjEquivOfDirectSecondPartDef
G4double currentTcutForDirectSecond
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
void SetApplyCutInRange(G4bool aBool)
void ProposeTrackStatus(G4TrackStatus status)
virtual G4double GetSecondAdjEnergyMinForProdToProjCase(G4double PrimAdjEnergy)
virtual G4double DiffCrossSectionPerAtomPrimToSecond(G4double kinEnergyProj, G4double kinEnergyProd, G4double Z, G4double A=0.)
static G4AdjointCSManager * GetAdjointCSManager()
double epsilon(double density, double temperature)
virtual void CorrectPostStepWeight(G4ParticleChange *fParticleChange, G4double old_weight, G4double adjointPrimKinEnergy, G4double projectileKinEnergy, G4bool IsScatProjToProjCase)
virtual G4double DiffCrossSectionPerAtomPrimToScatPrim(G4double kinEnergyProj, G4double kinEnergyScatProj, G4double Z, G4double A=0.)
G4VEmProcess * theDirectEMProcess
G4bool second_part_of_same_type