72 G4bool IsScatProjToProjCase,
74 {
if (IsScatProjToProjCase) return ;
103 for (i=0; i<nShells-1; i++){
115 G4double gamma = 1. + electronEnergy/electron_mass_c2;
117 G4double beta = std::sqrt(gamma*gamma-1.)/gamma;
118 G4double b = 0.5*gamma*(gamma-1.)*(gamma-2);
121 if (gamma < 2.) grejsup = gamma*gamma*(1.+b-beta*b);
122 else grejsup = gamma*gamma*(1.+b+beta*b);
125 cos_theta = (rndm+beta)/(rndm*beta+1.);
126 term = 1.-beta*cos_theta;
127 greject = (1.-cos_theta*cos_theta)*(1.+b*term)/(term*term);
136 G4double sin_theta = std::sqrt(1.-cos_theta*cos_theta);
138 G4double dirx = sin_theta*std::cos(Phi),diry = sin_theta*std::sin(Phi),dirz = cos_theta;
140 adjoint_gammaDirection.rotateUz(electronDirection);
182 new_weight*=projectileKinEnergy/adjointPrimKinEnergy;
193 G4bool IsScatProjToProjCase)
197 if (IsScatProjToProjCase)
return 0.;
228 G4bool IsScatProjToProjCase)
240 G4double gammaEnergy = electronEnergy+B0;
243 if (CS >0) adjointCS += CS/gammaEnergy;
245 for (i=1;i<nShells;i++){
250 if (electronEnergy <Bi_-Bi) {
251 gammaEnergy = electronEnergy+Bi;
254 if (CS>0) adjointCS +=CS/gammaEnergy;
259 adjointCS*=electronEnergy;
static G4AdjointGamma * AdjointGamma()
void DefineCurrentMaterialAndElectronEnergy(const G4MaterialCutsCouple *aCouple, G4double eEnergy)
G4int GetNbOfAtomicShells() const
virtual void CorrectPostStepWeight(G4ParticleChange *fParticleChange, G4double old_weight, G4double adjointPrimKinEnergy, G4double projectileKinEnergy, G4bool IsScatProjToProjCase)
std::vector< G4Element * > G4ElementVector
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
G4double GetPostStepWeightCorrection()
const G4DynamicParticle * GetDynamicParticle() const
G4ParticleDefinition * theDirectPrimaryPartDef
~G4AdjointPhotoElectricModel()
void ProposeParentWeight(G4double finalWeight)
G4PEEffectFluoModel * theDirectPEEffectModel
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
static G4AdjointElectron * AdjointElectron()
virtual G4double GetAdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)
G4Material * currentMaterial
G4double totBiasedAdjointCS
G4ParticleDefinition * theAdjEquivOfDirectPrimPartDef
const G4ElementVector * GetElementVector() const
G4double pre_step_AdjointCS
void SetSecondaryWeightByProcess(G4bool)
void SetParentWeightByProcess(G4bool)
const G4double * GetVecNbOfAtomsPerVolume() const
virtual G4double AdjointCrossSection(const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool IsScatProjToProjCase)
G4double post_step_AdjointCS
const G4ThreeVector & GetMomentumDirection() const
static const double twopi
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A, G4double, G4double)
G4double AdjointCrossSectionPerAtom(const G4Element *anElement, G4double electronEnergy)
void SetUseMatrix(G4bool aBool)
G4MaterialCutsCouple * currentCouple
virtual void SampleSecondaries(const G4Track &aTrack, G4bool IsScatProjToProjCase, G4ParticleChange *fParticleChange)
G4double shell_prob[40][40]
void SetCurrentCouple(const G4MaterialCutsCouple *)
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4AdjointPhotoElectricModel()
void AddSecondary(G4Track *aSecondary)
G4double GetWeight() const
G4ParticleDefinition * theAdjEquivOfDirectSecondPartDef
size_t currentCoupleIndex
size_t GetNumberOfElements() const
G4double GetAtomicShell(G4int index) const
void SetApplyCutInRange(G4bool aBool)
void ProposeTrackStatus(G4TrackStatus status)
static G4AdjointCSManager * GetAdjointCSManager()
size_t currentMaterialIndex
const G4Material * GetMaterial() const
G4bool second_part_of_same_type