69 G4String materialName = materialNameIni;
70 if(materialName.find(
"_MODIFIED")){
71 materialName = materialName.substr(0,materialName.size()-9);
86 G4cerr<<
"FatalError. Auger kineticEnergy: "<<kineticEnergy<<
G4endl;
90 if(atomId==1 || atomId==2 || atomId==3)
106 if(materialName==
"THF" || materialName==
"backbone_THF"){
114 else if(materialName==
"PY" || materialName==
"PU"
115 || materialName==
"cytosine_PY" || materialName==
"thymine_PY"
116 || materialName==
"adenine_PU" || materialName==
"guanine_PU"
126 else if(materialName==
"TMP"|| materialName==
"backbone_TMP"){
142 kineticEnergy = 495*
eV;
157 else if(atomId == 4){
173 kineticEnergy = pow(10, f1*pow( abs( log10(
Y) ) , g1) + f2*pow( abs( log10(
Y) ) , g2) + f3 )*
eV;
176 return kineticEnergy;
192 G4double newsinTh = std::sqrt(1.-newcosTh*newcosTh);
195 G4double xDir = newsinTh*std::sin(newPhi);
196 G4double yDir = newsinTh*std::cos(newPhi);
203 fvect->push_back(dp);
G4double Y(G4double density)
static constexpr double twopi
static constexpr double eV
G4GLOB_DLL std::ostream G4cerr
G4GLOB_DLL std::ostream G4cout
const G4String modelName
name of the auger model
virtual ~G4DNAPTBAugerModel()
~G4DNAPTBAugerModel Destructor
void GenerateAugerWithRandomDirection(std::vector< G4DynamicParticle * > *fvect, G4double kineticEnergy)
GenerateAugerWithRandomDirection Generates the auger particle.
G4int DetermineIonisedAtom(G4int atomId, const G4String &materialName, G4double bindingEnergy)
DetermineIonisedAtom.
void ComputeAugerEffect(std::vector< G4DynamicParticle * > *fvect, const G4String &materialNameIni, G4double bindingEnergy)
ComputeAugerEffect Main method to be called by the ionisation model.
void SetCutForAugerElectrons(G4double cut)
SetCutForAugerElectrons Set the cut for the auger electrons production.
G4double CalculAugerEnergyFor(G4int atomId)
CalculAugerEnergyFor.
virtual void Initialise()
Initialise Set the verbose value.
G4DNAPTBAugerModel(const G4String &modelName)
G4DNAPTBAugerModel Constructor.
G4double minElectronEnergy
static G4Electron * Electron()
G4double bindingEnergy(G4int A, G4int Z)