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"){
107 if(bindingEnergy==305.07){
110 else if(bindingEnergy==557.94){
114 else if(materialName==
"PY" || materialName==
"PU"
115 || materialName==
"cytosine_PY" || materialName==
"thymine_PY"
116 || materialName==
"adenine_PU" || materialName==
"guanine_PU"
119 if(bindingEnergy==307.52){
122 else if(bindingEnergy==423.44){
126 else if(materialName==
"TMP"|| materialName==
"backbone_TMP"){
127 if(bindingEnergy==209.59 || bindingEnergy==152.4)
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;
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);
void ComputeAugerEffect(std::vector< G4DynamicParticle * > *fvect, const G4String &materialNameIni, G4double bindingEnergy)
ComputeAugerEffect Main method to be called by the ionisation model.
virtual ~G4DNAPTBAugerModel()
~G4DNAPTBAugerModel Destructor
G4double CalculAugerEnergyFor(G4int atomId)
CalculAugerEnergyFor.
static constexpr double twopi
static constexpr double eV
static G4Electron * Electron()
G4GLOB_DLL std::ostream G4cerr
G4int DetermineIonisedAtom(G4int atomId, const G4String &materialName, G4double bindingEnergy)
DetermineIonisedAtom.
const G4String modelName
name of the auger model
G4DNAPTBAugerModel(const G4String &modelName)
G4DNAPTBAugerModel Constructor.
G4GLOB_DLL std::ostream G4cout
G4double minElectronEnergy
virtual void Initialise()
Initialise Set the verbose value.
void GenerateAugerWithRandomDirection(std::vector< G4DynamicParticle * > *fvect, G4double kineticEnergy)
GenerateAugerWithRandomDirection Generates the auger particle.
G4double bindingEnergy(G4int A, G4int Z)
void SetCutForAugerElectrons(G4double cut)
SetCutForAugerElectrons Set the cut for the auger electrons production.